An interferon-inducible signature of airway disease from blood gene expression profiling

BACKGROUND

The molecular basis of airway remodelling in chronic obstructive pulmonary disease (COPD) remains poorly understood. We identified gene expression signatures associated with chest computed tomography (CT) scan airway measures to understand molecular pathways associated with airway disease.

METHODS

In 2396 subjects in the COPDGene Study, we examined the relationship between quantitative CT airway phenotypes and blood transcriptomes to identify airway disease-specific genes and to define an airway wall thickness (AWT) gene set score. Multivariable regression analyses were performed to identify associations of the AWT score with clinical phenotypes, bronchial gene expression and genetic variants.

RESULTS

Type 1 interferon (IFN)-induced genes were consistently associated with AWT, square root wall area of a hypothetical airway with 10 mm internal perimeter (Pi10) and wall area percentage, with the strongest enrichment in AWT. A score derived from 18 genes whose expression was associated with AWT was associated with COPD-related phenotypes including reduced lung function (forced expiratory volume in 1 s percentage predicted β= -3.4; p<0.05) and increased exacerbations (incidence rate ratio 1.7; p<0.05). The AWT score was reproducibly associated with AWT in bronchial samples from 23 subjects (β=3.22; p<0.05). The blood AWT score was associated with genetic variant rs876039, an expression quantitative trait locus for IKZF1, a gene that regulates IFN signalling and is associated with inflammatory diseases.

CONCLUSIONS

A gene expression signature with IFN-stimulated genes from peripheral blood and bronchial brushings is associated with CT AWT, lung function and exacerbations. Shared genes and genetic associations suggest viral responses and/or autoimmune dysregulation as potential underlying mechanisms of airway disease in COPD.

Lung tissue shows divergent gene expression between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis

BACKGROUND

Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF.

METHODS

We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets.

RESULTS

We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts.

CONCLUSIONS

We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

Blood RNA sequencing shows overlapping gene expression across COPD phenotype domains

RATIONALE

COPD can be assessed using multidimensional grading systems with components from three domains: pulmonary function tests, symptoms and systemic features. Clinically, measures may be used interchangeably, though it is not known if they share similar pathobiology.

OBJECTIVE

To use RNA sequencing (RNA-seq) to determine if there is an overlap in the underlying biological mechanisms and consequences driving different components of the multidimensional grading systems.

METHODS

Whole blood was collected for RNA-seq from current and former smokers in the Genetic Epidemiology of COPD study. We tested the overlap in gene expression and biological pathways associated with case-control status and quantitative COPD phenotypes within and between the three domains.

RESULTS

In 2647 subjects, there were 3030 genes differentially expressed in any of the three domains or case-control status. There were five genes that overlapped between the three domains and case-control status, including G protein-coupled receptor 15(GPR15), sestrin 1 (SESN1) and interferon-induced guanylate-binding protein 1 (GBP1), which were associated with longitudinal decline in FEV1. The overlap between the three domains was enriched for pathways related to cellular components.

CONCLUSIONS

We identified gene sets and pathways that overlap between 12 COPD-related phenotypes and case-control status. There were no pathways represented in the overlap between the three domains and case-control status, but we identified multiple genes that demonstrated a consistent pattern of expression across several of the phenotypes. Patterns of gene expression correlation were generally similar to the correlation of clinical phenotypes in the PFT and symptom domains but not the systemic features.

Zonisamide safety in pregnancy: Data from the UK and Ireland epilepsy and pregnancy register

BACKGROUND

Animal data suggest teratogenic effects with zonisamide use and risk of pregnancy losses. Human data following zonisamide exposure are presently limited, but suggest low risk of malformation with elevated risk of low birth weight.

OBJECTIVE

To calculate the major congenital malformation (MCM) rate of zonisamide in human pregnancy and assess for a signal of any specific malformation pattern and associations with birth weight.

METHODS AND MATERIALS

Data were obtained from the UK and Ireland Epilepsy and Pregnancy register (UKIEPR) which is an observational, registration, and follow up study from December 1996 to July 2020. Eligibility criteria were use of zonisamide and to have been referred to the UKIEPR before the outcome of the pregnancy was known. Primary outcome was evidence of MCM.

RESULTS

From December 1996 through July 2020 there were 112 cases of first trimester exposure to zonisamide, including 26 monotherapy cases. There were 3 MCM for monotherapy cases (MCM rate 13.0% (95% confidence interval 4.5-32.1)), and 5 MCM for polytherapy cases (MCM rate 6.9% (95% confidence interval 3.0-15.2)). While the median birth weight was on 71st and 44th centile for monotherapy and polytherapy cases respectively, there was a high rate of infants born small for gestational age (21% for both).

CONCLUSION

These data raise concerns about a signal for potential teratogenicity with zonisamide in human pregnancy. Given the low numbers reported, further data will be required to adequately counsel women who use zonisamide in pregnancy.

Sex-specific associations with DNA methylation in lung tissue demonstrate smoking interactions

Cigarette smoking impacts DNA methylation, but the investigation of sex-specific features of lung tissue DNA methylation in smokers has been limited. Women appear more susceptible to cigarette smoke, and often develop more severe lung disease at an earlier age with less smoke exposure. We aimed to analyse whether there are sex differences in DNA methylation in lung tissue and whether these DNA methylation marks interact with smoking. We collected lung tissue samples from former smokers who underwent lung tissue resection. One hundred thirty samples from white subjects were included for this analysis. Regression models for sex as a predictor of methylation were adjusted for age, presence of COPD, smoking variables and technical batch variables revealed 710 associated sites. 294 sites demonstrated robust sex-specific methylation associations in foetal lung tissue. Pathway analysis identified 6 nominally significant pathways including the mitophagy pathway. Three CpG sites demonstrated a suggested interaction between sex and pack-years of smoking: GPR132, ANKRD44 and C19orf60. All of them were nominally significant in both male- and female-specific models, and the effect estimates were in opposite directions for male and female; GPR132 demonstrated significant association between DNA methylation and gene expression in lung tissue (P < 0.05). Sex-specific associations with DNA methylation in lung tissue are wide-spread and may reveal genes and pathways relevant to sex differences for lung damaging effects of cigarette smoking.

ScribeMD: Mixed-methods analysis of a pen pal program in undergraduate medical education

Introduction: A growing emphasis on humanistic medical care has led to the development of programs to imbue more humanistic values into training physicians. The ScribeMD intra-class pen pal program is a narrative medicine tool designed to meet this need with a focus on the journey through medical education. This study intends to evaluate the efficacy of this novel program on the participants' professional identity formation (PIF), emotional intelligence development (EID), and narrative competency improvement (NCI).Methods: This program evaluation involved quantitative and qualitative survey solicitations from first and second-year medical student participants and controls. An efficacy-triangulation model was developed as a quantitative outcome-measuring tool linking objective, writing prompts, and survey data.Results: The quantitative results showed statistically significant improvements in line with the Efficacy Triangulation Model in the participant group. Qualitatively, reflective and insight-driven gains were elicited, as were notable themes of personal and community improvements.Conclusion: ScribeMD poses a unique avenue for the development of more humanistic physicians during their medical education. We recommend additional development of the Efficacy Triangulation Model as well as the program design itself. We also recommend the piloting of this program in other healthcare-education settings both with medical education and other health professions.

Integrated transcriptomic correlation network analysis identifies COPD molecular determinants

Chronic obstructive pulmonary disease (COPD) is a complex and heterogeneous syndrome. Network-based analysis implemented by SWIM software can be exploited to identify key molecular switches - called "switch genes" - for the disease. Genes contributing to common biological processes or defining given cell types are usually co-regulated and co-expressed, forming expression network modules. Consistently, we found that the COPD correlation network built by SWIM consists of three well-characterized modules: one populated by switch genes, all up-regulated in COPD cases and related to the regulation of immune response, inflammatory response, and hypoxia (like TIMP1, HIF1A, SYK, LY96, BLNK and PRDX4); one populated by well-recognized immune signature genes, all up-regulated in COPD cases; one where the GWAS genes AGER and CAVIN1 are the most representative module genes, both down-regulated in COPD cases. Interestingly, 70% of AGER negative interactors are switch genes including PRDX4, whose activation strongly correlates with the activation of known COPD GWAS interactors SERPINE2, CD79A, and POUF2AF1. These results suggest that SWIM analysis can identify key network modules related to complex diseases like COPD.

Whole exome sequencing analysis in severe chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD), one of the leading causes of death worldwide, is substantially influenced by genetic factors. Alpha-1 antitrypsin deficiency demonstrates that rare coding variants of large effect can influence COPD susceptibility. To identify additional rare coding variants in patients with severe COPD, we conducted whole exome sequencing analysis in 2543 subjects from two family-based studies (Boston Early-Onset COPD Study and International COPD Genetics Network) and one case-control study (COPDGene). Applying a gene-based segregation test in the family-based data, we identified significant segregation of rare loss of function variants in TBC1D10A and RFPL1 (P-value < 2x10-6), but were unable to find similar variants in the case-control study. In single-variant, gene-based and pathway association analyses, we were unable to find significant findings that replicated or were significant in meta-analysis. However, we found that the top results in the two datasets were in proximity to each other in the protein-protein interaction network (P-value = 0.014), suggesting enrichment of these results for similar biological processes. A network of these association results and their neighbors was significantly enriched in the transforming growth factor beta-receptor binding and cilia-related pathways. Finally, in a more detailed examination of candidate genes, we identified individuals with putative high-risk variants, including patients harboring homozygous mutations in genes associated with cutis laxa and Niemann-Pick Disease Type C. Our results likely reflect heterogeneity of genetic risk for COPD along with limitations of statistical power and functional annotation, and highlight the potential of network analysis to gain insight into genetic association studies.

Transcriptomic Analysis of Lung Tissue from Cigarette Smoke-Induced Emphysema Murine Models and Human Chronic Obstructive Pulmonary Disease Show Shared and Distinct Pathways

Although cigarette smoke (CS) is the primary risk factor for chronic obstructive pulmonary disease (COPD), the underlying molecular mechanisms for the significant variability in developing COPD in response to CS are incompletely understood. We performed lung gene expression profiling of two different wild-type murine strains (C57BL/6 and NZW/LacJ) and two genetic models with mutations in COPD genome-wide association study genes (HHIP and FAM13A) after 6 months of chronic CS exposure and compared the results to human COPD lung tissues. We identified gene expression patterns that correlate with severity of emphysema in murine and human lungs. Xenobiotic metabolism and nuclear erythroid 2-related factor 2-mediated oxidative stress response were commonly regulated molecular response patterns in C57BL/6, Hhip^+/-, and Fam13a^-/- murine strains exposed chronically to CS. The CS-resistant Fam13a^-/- mouse and NZW/LacJ strain revealed gene expression response pattern differences. The Fam13a^-/- strain diverged in gene expression compared with C57BL/6 control only after CS exposure. However, the NZW/LacJ strain had a unique baseline expression pattern, enriched for nuclear erythroid 2-related factor 2-mediated oxidative stress response and xenobiotic metabolism, and converged to a gene expression pattern similar to the more susceptible wild-type C57BL/6 after CS exposure. These results suggest that distinct molecular pathways may account for resistance to emphysema. Surprisingly, there were few genes commonly modulated in mice and humans. Our study suggests that gene expression responses to CS may be largely species and model dependent, yet shared pathways could provide biologically significant insights underlying individual susceptibility to CS.

Exome Sequencing Analysis in Severe, Early-Onset Chronic Obstructive Pulmonary Disease

RATIONALE

Genomic regions identified by genome-wide association studies explain only a small fraction of heritability for chronic obstructive pulmonary disease (COPD). Alpha-1 antitrypsin deficiency shows that rare coding variants of large effect also influence COPD susceptibility. We hypothesized that exome sequencing in families identified through a proband with severe, early-onset COPD would identify additional rare genetic determinants of large effect.

OBJECTIVES

To identify rare genetic determinants of severe COPD.

METHODS

We applied filtering approaches to identify potential causal variants for COPD in whole exomes from 347 subjects in 49 extended pedigrees from the Boston Early-Onset COPD Study. We assessed the power of this approach under different levels of genetic heterogeneity using simulations. We tested genes identified in these families using gene-based association tests in exomes of 204 cases with severe COPD and 195 resistant smokers from the COPDGene study. In addition, we examined previously described loci associated with COPD using these datasets.

MEASUREMENTS AND MAIN RESULTS

We identified 69 genes with predicted deleterious nonsynonymous, stop, or splice variants that segregated with severe COPD in at least two pedigrees. Four genes (DNAH8, ALCAM, RARS, and GBF1) also demonstrated an increase in rare nonsynonymous, stop, and/or splice mutations in cases compared with resistant smokers from the COPDGene study; however, these results were not statistically significant. We demonstrate the limitations of the power of this approach under genetic heterogeneity through simulation.

CONCLUSIONS

Rare deleterious coding variants may increase risk for COPD, but multiple genes likely contribute to COPD susceptibility.

Detecting disease-associated genomic outcomes using constrained mixture of Bayesian hierarchical models for paired data

Detecting disease-associated genomic outcomes is one of the key steps in precision medicine research. Cutting-edge high-throughput technologies enable researchers to unbiasedly test if genomic outcomes are associated with disease of interest. However, these technologies also include the challenges associated with the analysis of genome-wide data. Two big challenges are (1) how to reduce the effects of technical noise; and (2) how to handle the curse of dimensionality (i.e., number of variables are way larger than the number of samples). To tackle these challenges, we propose a constrained mixture of Bayesian hierarchical models (MBHM) for detecting disease-associated genomic outcomes for data obtained from paired/matched designs. Paired/matched designs can effectively reduce effects of confounding factors. MBHM does not involve multiple testing, hence does not have the problem of the curse of dimensionality. It also could borrow information across genes so that it can be used for whole genome data with small sample sizes.

Gene-based segregation method for identifying rare variants in family-based sequencing studies

Whole-exome sequencing using family data has identified rare coding variants in Mendelian diseases or complex diseases with Mendelian subtypes, using filters based on variant novelty, functionality, and segregation with the phenotype within families. However, formal statistical approaches are limited. We propose a gene-based segregation test (GESE) that quantifies the uncertainty of the filtering approach. It is constructed using the probability of segregation events under the null hypothesis of Mendelian transmission. This test takes into account different degrees of relatedness in families, the number of functional rare variants in the gene, and their minor allele frequencies in the corresponding population. In addition, a weighted version of this test allows incorporating additional subject phenotypes to improve statistical power. We show via simulations that the GESE and weighted GESE tests maintain appropriate type I error rate, and have greater power than several commonly used region-based methods. We apply our method to whole-exome sequencing data from 49 extended pedigrees with severe, early-onset chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD study (BEOCOPD) and identify several promising candidate genes. Our proposed methods show great potential for identifying rare coding variants of large effect and high penetrance for family-based sequencing data. The proposed tests are implemented in an R package that is available on CRAN (https://cran.r-project.org/web/packages/GESE/).

Differential DNA methylation marks and gene comethylation of COPD in African-Americans with COPD exacerbations

Background

Chronic obstructive pulmonary disease (COPD) is the third-leading cause of death worldwide. Identifying COPD-associated DNA methylation marks in African-Americans may contribute to our understanding of racial disparities in COPD susceptibility. We determined differentially methylated genes and co-methylation network modules associated with COPD in African-Americans recruited during exacerbations of COPD and smoking controls from the Pennsylvania Study of Chronic Obstructive Pulmonary Exacerbations (PA-SCOPE) cohort.

Methods

We assessed DNA methylation from whole blood samples in 362 African-American smokers in the PA-SCOPE cohort using the Illumina Infinium HumanMethylation27 BeadChip Array. Final analysis included 19302 CpG probes annotated to the nearest gene transcript after quality control. We tested methylation associations with COPD case-control status using mixed linear models. Weighted gene comethylation networks were constructed using weighted gene coexpression network analysis (WGCNA) and network modules were analyzed for association with COPD.

Results

There were five differentially methylated CpG probes significantly associated with COPD among African-Americans at an FDR less than 5 %, and seven additional probes that approached significance at an FDR less than 10 %. The top ranked gene association was MAML1, which has been shown to affect NOTCH-dependent angiogenesis in murine lung. Network modeling yielded the “yellow” and “blue” comethylation modules which were significantly associated with COPD (p-value 4 × 10^-10 and 4 × 10^-9, respectively). The yellow module was enriched for gene sets related to inflammatory pathways known to be relevant to COPD. The blue module contained the top ranked genes in the concurrent differential methylation analysis (FXYD1/LGI4, gene significance p-value 1.2 × 10^-26; MAML1, p-value 2.0 × 10^-26; CD72, p-value 2.1 × 10^-25; and LPO, p-value 7.2 × 10^-25), and was significantly associated with lung development processes in Gene Ontology gene-set enrichment analysis.

Conclusion

We identified 12 differentially methylated CpG sites associated with COPD that mapped to biologically plausible genes. Network module comethylation patterns have identified candidate genes that may be contributing to racial differences in COPD susceptibility and severity. COPD-associated comethylation modules contained genes previously associated with lung disease and inflammation and recapitulated known COPD-associated genes. The genes implicated by differential methylation and WGCNA analysis may provide mechanistic targets contributing to COPD susceptibility, exacerbations, and outcomes among African-Americans.

Trial registration

Trial Registration: NCT00774176, Registry: ClinicalTrials.gov, URL: www.clinicaltrials.gov, Date of Enrollment of First Participant: June 2004, Date Registered: 04 January 2008 (retrospectively registered).

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-016-0459-8) contains supplementary material, which is available to authorized users.

The impact of genetic variation and cigarette smoke on DNA methylation in current and former smokers from the COPDGene study

DNA methylation can be affected by systemic exposures, such as cigarette smoking and genetic sequence variation; however, the relative impact of each on the epigenome is unknown. We aimed to assess if cigarette smoking and genetic variation are associated with overlapping or distinct sets of DNA methylation marks and pathways. We selected 85 Caucasian current and former smokers with genome-wide single nucleotide polymorphism (SNP) genotyping available from the COPDGene study. Genome-wide methylation was obtained on DNA from whole blood using the Illumina HumanMethylation27 platform. To determine the impact of local sequence variation on DNA methylation (mQTL), we examined the association between methylation and SNPs within 50 kb of each CpG site. To examine the impact of cigarette smoking on DNA methylation, we examined the differences in methylation by current cigarette smoking status. We detected 770 CpG sites annotated to 708 genes associated at an FDR < 0.05 in the cis-mQTL analysis and 1,287 CpG sites annotated to 1,242 genes, which were nominally associated in the smoking-CpG association analysis (P(unadjusted) < 0.05). Forty-three CpG sites annotated to 40 genes were associated with both SNP variation and current smoking; this overlap was not greater than that expected by chance. Our results suggest that cigarette smoking and genetic variants impact distinct sets of DNA methylation marks, the further elucidation of which may partially explain the variable susceptibility to the health effects of cigarette smoking. Ascertaining how genetic variation and systemic exposures differentially impact the human epigenome has relevance for both biomarker identification and therapeutic target development for smoking-related diseases.

The Role of Vitamin D in the Transcriptional Program of Human Pregnancy

Background

Patterns of gene expression of human pregnancy are poorly understood. In a trial of vitamin D supplementation in pregnant women, peripheral blood transcriptomes were measured longitudinally on 30 women and used to characterize gene co-expression networks.

Objective

Studies suggest that increased maternal Vitamin D levels may reduce the risk of asthma in early life, yet the underlying mechanisms have not been examined. In this study, we used a network-based approach to examine changes in gene expression profiles during the course of normal pregnancy and evaluated their association with maternal Vitamin D levels.

Design

The VDAART study is a randomized clinical trial of vitamin D supplementation in pregnancy for reduction of pediatric asthma risk. The trial enrolled 881 women at 10–18 weeks of gestation. Longitudinal gene expression measures were obtained on thirty pregnant women, using RNA isolated from peripheral blood samples obtained in the first and third trimesters. Differentially expressed genes were identified using significance of analysis of microarrays (SAM), and clustered using a weighted gene co-expression network analysis (WGCNA). Gene-set enrichment was performed to identify major biological pathways.

Results

Comparison of transcriptional profiles between first and third trimesters of pregnancy identified 5839 significantly differentially expressed genes (FDR<0.05). Weighted gene co-expression network analysis clustered these transcripts into 14 co-expression modules of which two showed significant correlation with maternal vitamin D levels. Pathway analysis of these two modules revealed genes enriched in immune defense pathways and extracellular matrix reorganization as well as genes enriched in notch signaling and transcription factor networks.

Conclusion

Our data show that gene expression profiles of healthy pregnant women change during the course of pregnancy and suggest that maternal Vitamin D levels influence transcriptional profiles. These alterations of the maternal transcriptome may contribute to fetal immune imprinting and reduce allergic sensitization in early life.

Trial Registration

clinicaltrials.gov NCT00920621

Exome Array Analysis Identifies a Common Variant in IL27 Associated with Chronic Obstructive Pulmonary Disease

RATIONALE

Chronic obstructive pulmonary disease (COPD) susceptibility is in part related to genetic variants. Most genetic studies have been focused on genome-wide common variants without a specific focus on coding variants, but common and rare coding variants may also affect COPD susceptibility.

OBJECTIVES

To identify coding variants associated with COPD.

METHODS

We tested nonsynonymous, splice, and stop variants derived from the Illumina HumanExome array for association with COPD in five study populations enriched for COPD. We evaluated single variants with a minor allele frequency greater than 0.5% using logistic regression. Results were combined using a fixed effects meta-analysis. We replicated novel single-variant associations in three additional COPD cohorts.

MEASUREMENTS AND MAIN RESULTS

We included 6,004 control subjects and 6,161 COPD cases across five cohorts for analysis. Our top result was rs16969968 (P = 1.7 × 10(-14)) in CHRNA5, a locus previously associated with COPD susceptibility and nicotine dependence. Additional top results were found in AGER, MMP3, and SERPINA1. A nonsynonymous variant, rs181206, in IL27 (P = 4.7 × 10(-6)) was just below the level of exome-wide significance but attained exome-wide significance (P = 5.7 × 10(-8)) when combined with results from other cohorts. Gene expression datasets revealed an association of rs181206 and the surrounding locus with expression of multiple genes; several were differentially expressed in COPD lung tissue, including TUFM.

CONCLUSIONS

In an exome array analysis of COPD, we identified nonsynonymous variants at previously described loci and a novel exome-wide significant variant in IL27. This variant is at a locus previously described in genome-wide associations with diabetes, inflammatory bowel disease, and obesity and appears to affect genes potentially related to COPD pathogenesis.

A Comparative Study of Tests for Homogeneity of Variances with Application to DNA Methylation Data

Variable DNA methylation has been associated with cancers and complex diseases. Researchers have identified many DNA methylation markers that have different mean methylation levels between diseased subjects and normal subjects. Recently, researchers found that DNA methylation markers with different variabilities between subject groups could also have biological meaning. In this article, we aimed to help researchers choose the right test of equal variance in DNA methylation data analysis. We performed systematic simulation studies and a real data analysis to compare the performances of 7 equal-variance tests, including 2 tests recently proposed in the DNA methylation analysis literature. Our results showed that the Brown-Forsythe test and trimmed-mean-based Levene's test had good performance in testing for equality of variance in our simulation studies and real data analyses. Our results also showed that outlier profiles could be biologically very important.

Smoking-Associated Site-Specific Differential Methylation in Buccal Mucosa in the COPDGene Study

DNA methylation is a complex, tissue-specific phenomenon that can reflect both endogenous factors and exogenous exposures. Buccal brushings represent an easily accessible source of DNA, which may be an appropriate surrogate tissue in the study of environmental exposures and chronic respiratory diseases. Buccal brushings were obtained from a subset of current and former smokers from the COPDGene study. Genome-wide DNA methylation data were obtained in the discovery cohort (n = 82) using the Illumina HumanMethylation450K array. Empirical Bayes methods were used to test for differential methylation by current smoking status at 468,219 autosomal CpG sites using linear models adjusted for age, sex, and race. Pyrosequencing was performed in a nonoverlapping replication cohort (n = 130). Current smokers were significantly younger than former smokers in both the discovery and replication cohorts. Seven CpG sites were associated with current smoking at a false discovery rate less than 0.05 in the discovery cohort. Six of the seven significant sites were pyrosequenced in the replication cohort; five CpG sites, including sites annotated to CYP1B1 and PARVA, were replicated. Correlations between cumulative smoke exposure and time since smoking cessation were observed in a subset of the significantly associated CpG sites. A significant correlation between reduced lung function and increased radiographic emphysema with methylation at cg02162897 (CYP1B1) was observed among female subjects. Site-specific methylation of DNA isolated from buccal mucosa is associated with exposure to cigarette smoke, and may provide insights into the mechanisms underlying differential susceptibility toward the development of smoking-related chronic respiratory diseases.

Fetal lung and placental methylation is associated with in utero nicotine exposure

In utero smoke exposure has been shown to have detrimental effects on lung function and to be associated with persistent wheezing and asthma in children. One potential mechanism of IUS effects could be alterations in DNA methylation, which may have life-long implications. The goal of this study was to examine the association between DNA methylation and nicotine exposure in fetal lung and placental tissue in early development; nicotine exposure in this analysis represents a likely surrogate for in-utero smoke. We performed an epigenome-wide analysis of DNA methylation in fetal lung tissue (n = 85, 41 smoke exposed (48%), 44 controls) and the corresponding placental tissue samples (n = 80, 39 smoke exposed (49%), 41 controls) using the Illumina HumanMethylation450 BeadChip array. Differential methylation analyses were conducted to evaluate the variation associated with nicotine exposure. The most significant CpG sites in the fetal lung analysis mapped to the PKP3 (P = 2.94 × 10(-03)), ANKRD33B (P = 3.12 × 10(-03)), CNTD2 (P = 4.9 × 10(-03)) and DPP10 (P = 5.43 × 10(-03)) genes. In the placental methylome, the most significant CpG sites mapped to the GTF2H2C and GTF2H2D genes (P = 2.87 × 10(-06) - 3.48 × 10(-05)). One hundred and one unique CpG sites with P-values < 0.05 were concordant between lung and placental tissue analyses. Gene Set Enrichment Analysis demonstrated enrichment of specific disorders, such as asthma and immune disorders. Our findings demonstrate an association between in utero nicotine exposure and variable DNA methylation in fetal lung and placental tissues, suggesting a role for DNA methylation variation in the fetal origins of chronic diseases.

Genome-wide site-specific differential methylation in the blood of individuals with Klinefelter syndrome

Klinefelter syndrome (KS) (47 XXY) is a common sex-chromosome aneuploidy with an estimated prevalence of one in every 660 male births. Investigations into the associations between DNA methylation and the highly variable clinical manifestations of KS have largely focused on the supernumerary X chromosome; systematic investigations of the epigenome have been limited. We obtained genome-wide DNA methylation data from peripheral blood using the Illumina HumanMethylation450K platform in 5 KS (47 XXY) versus 102 male (46 XY) and 113 female (46 XX) control subjects participating in the COPDGene Study. Empirical Bayes-mediated models were used to test for differential methylation by KS status. CpG sites with a false-discovery rate < 0.05 in the discovery cohort which were available on the first-generation HumanMethylation 27 K platform were further examined in an independent replication cohort of 2 KS subjects, 590 male, and 495 female controls drawn from the International COPD Genetics Network (ICGN). Differential methylation at sites throughout the genome were identified, including 86 CpG sites that were differentially methylated in KS subjects relative to both male and female controls. CpG sites annotated to the HEN1 methyltransferase homolog 1 (HENMT1), calcyclin-binding protein (CACYBP), and GTPase-activating protein (SH3 domain)-binding protein 1 (G3BP1) genes were among the "KS-specific" loci that were replicated in ICGN. We conclude that site-specific differential methylation exists throughout the genome in KS. The functional impact and clinical relevance of these differentially methylated loci should be explored in future studies.

Malformation risks of antiepileptic drug monotherapies in pregnancy: updated results from the UK and Ireland Epilepsy and Pregnancy Registers

OBJECTIVES

Antiepileptic drug (AED) exposure during pregnancy increases the risk of major congenital malformations (MCMs). The magnitude of this risk varies by AED exposure. Here we provide updated results from the UK Epilepsy and Pregnancy Register of the risk of MCMs after monotherapy exposure to valproate, carbamazepine and lamotrigine.

METHODS

Fifteen-year prospective observational study from 1996 until 2012. The main outcome measure is the MCM rate.

RESULTS

Informative outcomes were available for 5206 cases. 1290 women were exposed to valproate monotherapy, 1718 to carbamazepine monotherapy and 2198 to lamotrigine monotherapy. The MCM risk with valproate monotherapy exposure in utero was 6.7% (95% CI 5.5% to 8.3%) compared with 2.6% with carbamazepine (95% CI 1.9% to 3.5%) and 2.3% with lamotrigine (95% CI 1.8% to 3.1%). A significant dose effect was seen with valproate (p=0.0006) and carbamazepine (p=0.03) exposed pregnancies. A non-significant trend towards higher MCM rate with increasing dose was found with lamotrigine. MCM rate for high-dose lamotrigine (>400 mg daily) was lower than the MCM rate for pregnancies exposed to <600 mg daily of valproate, but this was not significant (3.4% vs 5.0%, p=0.31).

CONCLUSIONS

In utero exposure to valproate carries a significantly higher MCM risk than lamotrigine (p=0.0001) and carbamazepine (p=0.0001) monotherapy. In contrast to prior findings, high-dose lamotrigine was associated with fewer MCMs than all doses of valproate. While lamotrigine has a favourable profile compared with valproate for adverse pregnancy outcomes, the requirements for seizure control should not be overlooked.

The effect of socioeconomic status on treatment and pregnancy outcomes in women with epilepsy in Scotland

Compared to the background population, people with epilepsy tend to have lower rates of education and employment, lower rates of marriage and childbearing, and lower overall socioeconomic status (SES). Disparities in epilepsy care based on sociodemographic factors have been observed in the literature, but it is not known whether any such disparities exist in the UK. The UK Epilepsy and Pregnancy Register is a prospective, observational, registration and follow-up study that was set up to determine the relative safety of all AEDs taken in pregnancy. Here, we report outcomes of registered pregnancies to women with epilepsy living in Scotland from December 1996 to June 2012, based on the degree of socioeconomic deprivation of their postcode area. The Scottish Index of Multiple Deprivation (SIMD) quintile scores from 2006 were used to determine degree of socioeconomic deprivation, and group 1 (most deprived) and group 5 (least deprived) were compared. There were 1526 pregnancies with complete outcome data to women living in Scotland. Of these, 1453 (95.1%) resulted in a live birth and 68 (4.7%) had a major congenital malformation (MCM). Postcodes could not be reliably identified or verified for an additional three women, who have been excluded from SIMD group analysis. Of all women included in this study, 32.4% were in group 1 and 13.2% in group 5. No difference in MCM rate was observed between the two groups (4.4% in group 1 compared to 4.7% in group 5, p=0.84). Women in group 5 were more likely to take preconceptual folic acid (56.8% compared to 14.0%, relative risk: 4.1; 95% CI: 3.1-5.2) and less likely to have generalized tonic-clonic seizures in pregnancy (13.0% compared to 29.2%, relative risk: 0.4; 95% CI: 0.3-0.7) than those in group 1. Women in group 5 were more likely to be on monotherapy regimens (80.2% compared to 65.9%, relative risk: 1.2; 95% CI: 1.1-1.3), less likely to be on valproate (19.5% compared to 28.0%, p=0.05), and more likely to be on lower doses of the drug (825.9mg/day compared to 1012.0mg/day, p=0.05) compared to those in group 1. Although no change in MCM rate was seen based on SES, differences in treatment between socioeconomic groups do exist, particularly for preconceptual folic acid consumption, AED regimen, and seizure frequency. Greater emphasis on the importance of preconceptual counseling, both to discuss AED choice and folic acid intake, would be of benefit, particularly to those living in areas of high socioeconomic deprivation, to improve equity of healthcare delivery for women with epilepsy in Scotland.