Genome-Wide Interaction Analysis of Air Pollution Exposure and Childhood Asthma with Functional Follow-up.

RATIONALE: The evidence supporting an association between traffic-related air pollution exposure and incident childhood asthma is inconsistent and may depend on genetic factors. OBJECTIVES: To identify gene-environment interaction effects on childhood asthma using genome-wide single-nucleotide polymorphism (SNP) data and air pollution exposure. Identified loci were further analyzed at epigenetic and transcriptomic levels. METHODS: We used land use regression models to estimate individual air pollution exposure (represented by outdoor NO2 levels) at the birth address and performed a genome-wide interaction study for doctors' diagnoses of asthma up to 8 years in three European birth cohorts (n = 1,534) with look-up for interaction in two separate North American cohorts, CHS (Children's Health Study) and CAPPS/SAGE (Canadian Asthma Primary Prevention Study/Study of Asthma, Genetics and Environment) (n = 1,602 and 186 subjects, respectively). We assessed expression quantitative trait locus effects in human lung specimens and blood, as well as associations among air pollution exposure, methylation, and transcriptomic patterns. MEASUREMENTS AND MAIN RESULTS: In the European cohorts, 186 SNPs had an interaction P < 1 × 10-4 and a look-up evaluation of these disclosed 8 SNPs in 4 loci, with an interaction P < 0.05 in the large CHS study, but not in CAPPS/SAGE. Three SNPs within adenylate cyclase 2 (ADCY2) showed the same direction of the interaction effect and were found to influence ADCY2 gene expression in peripheral blood (P = 4.50 × 10-4). One other SNP with P < 0.05 for interaction in CHS, rs686237, strongly influenced UDP-Gal:betaGlcNAc ß-1,4-galactosyltransferase, polypeptide 5 (B4GALT5) expression in lung tissue (P = 1.18 × 10-17). Air pollution exposure was associated with differential discs, large homolog 2 (DLG2) methylation and expression. CONCLUSIONS: Our results indicated that gene-environment interactions are important for asthma development and provided supportive evidence for interaction with air pollution for ADCY2, B4GALT5, and DLG2.

Meta-analysis identifies seven susceptibility loci involved in the atopic march.

Eczema often precedes the development of asthma in a disease course called the 'atopic march'. To unravel the genes underlying this characteristic pattern of allergic disease, we conduct a multi-stage genome-wide association study on infantile eczema followed by childhood asthma in 12 populations including 2,428 cases and 17,034 controls. Here we report two novel loci specific for the combined eczema plus asthma phenotype, which are associated with allergic disease for the first time; rs9357733 located in EFHC1 on chromosome 6p12.3 (OR 1.27; P=2.1 × 10(-8)) and rs993226 between TMTC2 and SLC6A15 on chromosome 12q21.3 (OR 1.58; P=5.3 × 10(-9)). Additional susceptibility loci identified at genome-wide significance are FLG (1q21.3), IL4/KIF3A (5q31.1), AP5B1/OVOL1 (11q13.1), C11orf30/LRRC32 (11q13.5) and IKZF3 (17q21). We show that predominantly eczema loci increase the risk for the atopic march. Our findings suggest that eczema may play an important role in the development of asthma after eczema.

BMP2 response pattern in human lung fibroblasts predicts outcome in lung adenocarcinomas.

BACKGROUND: Bone morphogenetic proteins play important roles in development, morphogenesis and cancer. With this study we aimed to characterize the response of lung stromal fibroblasts to BMPs and their antagonists on a genome wide level and investigate its potential role in human lung adenocarcinomas. METHODS: We used an ex vivo culture model and measured gene expression changes in human lung fibroblasts after stimulation with BMPs and their antagonists using HEEBO microarrays. The in vitro data were correlated with in vivo observations in published expression datasets of human lung adenocarcinomas. RESULTS: We have systematically analyzed the response to BMP2, BMP4, BMP7 and their antagonists, Gremlin and Noggin, to define common and specific gene expression patterns. A BMP2 induced gene expression signature was defined, which is specific for stromal fibroblasts. Gene expression profiles from lung adenocarcinoma biopsies were analyzed to determine the prognostic significance of the "Fibroblast specific BMP2 induced gene list". This gene list successfully segregated patients with different prognostic outcome in 3 datasets. In a small dataset (Garber et al.) there was a strong trend for a worse prognosis of patients with adenocarcinomas of all stages over-expressing the "Fibroblast specific BMP2 induced gene list". In two larger datasets with stage I adenocarcinomas we observed a significantly worse disease-free (p = 0.002, Lee et al. and p = 0.002, Bhattacharjee et al.) and overall survival (p = 0.0002). CONCLUSIONS: The effects of BMPs and their antagonists are heterogeneous in different cell types. The gene expression pattern induced by BMP2 in primary lung fibroblasts may predict outcomes of patients with lung adenocarcinomas.

IgA measurements in over 12 000 Swedish twins reveal sex differential heritability and regulatory locus near CD30L.

In a broad attempt to improve the understanding of the genetic regulation of serum IgA levels, the heritability was estimated in over 12 000 Swedish twins, and a genome-wide association study was conducted in a subsample of 9617. Using the classical twin model the heritability was found to be significantly larger among females (61%) compared with males (21%), while contribution from shared environment (20%) was only seen for males. By modeling the genetic relationship matrix with IgA levels, we estimate that a substantial proportion (31%) of variance in IgA levels can ultimately be explained by the investigated SNPs. The genome-wide association study revealed significant association to two loci: (i) rs6928791 located on chromosome 6, 22 kb upstream of the gene SAM and SH3 domain containing 1 (SASH1) and (ii) rs13300483 on chromosome 9, situated 12 kb downstream the CD30 ligand (CD30L) encoding gene. The association to rs13300483 was replicated in two additional independent Swedish materials. The heritability of IgA levels is moderate and can partly be attributable to common variation in the CD30L locus.

DNA methylation in the Neuropeptide S Receptor 1 (NPSR1) promoter in relation to asthma and environmental factors.

Asthma and allergy are complex disorders influenced by both inheritance and environment, a relationship that might be further clarified by epigenetics. Neuropeptide S Receptor 1 (NPSR1) has been associated with asthma and allergy and a study suggested modulation of the genetic risk by environmental factors. We aimed to study DNA methylation in the promoter region of NPSR1 in relation to asthma and environmental exposures. Electrophoretic Mobility Shift Assay (EMSA) was used to investigate potential functional roles of both genotypes and methylation status in the NPSR1 promoter. DNA methylation was analysed using EpiTYPER in blood samples from two well-characterized cohorts; the BIOAIR study of severe asthma in adults and the Swedish birth cohort BAMSE. We observed that DNA methylation and genetic variants in the promoter influenced the binding of nuclear proteins to DNA, suggesting functional relevance. Significant, although small, differences in methylation were related to both adult severe asthma (p = 0.0001) and childhood allergic asthma (p = 0.01). Furthermore, DNA methylation was associated with exposures such as current smoking in adults for two CpG sites (p = 0.005 and 0.04), parental smoking during infancy in the children (p = 0.02) and in which month the sample was taken (p = 0.01). In summary, DNA methylation levels in the promoter of NPSR1 showed small but significant associations with asthma, both in adults and in children, and to related traits such as allergy and certain environmental exposures. Both genetic variation and the methylated state of CpG sites seem to have an effect on the binding of nuclear proteins in the regulatory region of NPSR1 suggesting complex regulation of this gene in asthma and allergy.

Characterization of EGFR and ErbB2 expression in atopic dermatitis patients.

Atopic dermatitis (AD) is one of the most common chronic inflammatory skin diseases in industrialized countries. To identify candidate genes involved in the pathogenesis of AD, we previously undertook a genome-wide approach using DNA microarrays. A transcript encoding the epidermal growth factor receptor (EGFR) was found to be among the down-regulated transcripts in AD skin. Here, we further investigated the expression pattern of two EGFR family members (EGFR and ErbB2) in AD skin on a protein level. Immunohistochemical (IHC) analysis of EGFR and ErbB2 showed decreased expression of EGFR and ErbB2 proteins in AD lesional skin as compared to skin from healthy individuals. Interestingly, we found that EGFR and ErbB2 were reciprocally expressed in an in vitro model of keratinocyte proliferation and differentiation, paralleling the expression patterns found in epidermis of healthy skin. The highest levels of EGFR transcripts were found in proliferating cells, while ErbB2 was found in differentiated cells. We show that blocking EGFR activity combined with co-stimulation of the Th2-cytokine IL4 in keratinocytes leads to induction of the inflammatory chemokine CCL26/eotaxin-3 in vitro. Accordingly, increased CCL26 transcriptional levels were observed in AD lesional skin. Taken together, suppression of EGFR may contribute to the pathogenesis of AD via the regulation of inflammatory chemokines.

Birt-Hogg-Dubé renal tumors are genetically distinct from other renal neoplasias and are associated with up-regulation of mitochondrial gene expression.

BACKGROUND: Germline mutations in the folliculin (FLCN) gene are associated with the development of Birt-Hogg-Dubé syndrome (BHDS), a disease characterized by papular skin lesions, a high occurrence of spontaneous pneumothorax, and the development of renal neoplasias. The majority of renal tumors that arise in BHDS-affected individuals are histologically similar to sporadic chromophobe renal cell carcinoma (RCC) and sporadic renal oncocytoma. However, most sporadic tumors lack FLCN mutations and the extent to which the BHDS-derived renal tumors share genetic defects associated with the sporadic tumors has not been well studied. METHODS: BHDS individuals were identified symptomatically and FLCN mutations were confirmed by DNA sequencing. Comparative gene expression profiling analyses were carried out on renal tumors isolated from individuals afflicted with BHDS and a panel of sporadic renal tumors of different subtypes using discriminate and clustering approaches. qRT-PCR was used to confirm selected results of the gene expression analyses. We further analyzed differentially expressed genes using gene set enrichment analysis and pathway analysis approaches. Pathway analysis results were confirmed by generation of independent pathway signatures and application to additional datasets. RESULTS: Renal tumors isolated from individuals with BHDS showed distinct gene expression and cytogenetic characteristics from sporadic renal oncocytoma and chromophobe RCC. The most prominent molecular feature of BHDS-derived kidney tumors was high expression of mitochondria-and oxidative phosphorylation (OXPHOS)-associated genes. This mitochondria expression phenotype was associated with deregulation of the PGC-1α-TFAM signaling axis. Loss of FLCN expression across various tumor types is also associated with increased nuclear mitochondrial gene expression. CONCLUSIONS: Our results support a genetic distinction between BHDS-associated tumors and other renal neoplasias. In addition, deregulation of the PGC-1α-TFAM signaling axis is most pronounced in renal tumors that harbor FLCN mutations and in tumors from other organs that have relatively low expression of FLCN. These results are consistent with the recently discovered interaction between FLCN and AMPK and support a model in which FLCN is a regulator of mitochondrial function.

Characterisation of hairy cell leukaemia by tiling resolution array-based comparative genome hybridisation: a series of 13 cases and review of the literature.

Little is known about the cytogenetic features and molecular mechanisms behind hairy cell leukaemia (HCL), despite the advances in diagnosis and treatment. Therefore, we used high-resolution genome-wide array-based comparative genomic hybridisation (array-CGH) and multiplex ligation-dependent probe amplification (MLPA) to characterise copy number alterations (CNAs) in DNA from 13 cases of HCL. We also summarise CNAs and cytogenetic features in 109 HCL cases comprising our 13 cases and 96 cases from the literature. Genomic array-CGH revealed imbalances in two out of 13 cases in addition to previously described copy number variants (CNVs) found in healthy individuals. In one case, a 700 kb deletion of 20q11.22 was detected encompassing ten characterised genes, among them the TP53INP2, DNCL2A and ITCH genes. In the second case, trisomy 5, and a deletion of 5p15.2 encompassing a non-characterised gene AY328033 was found. Altogether only 20/81 (25%) of all cases studied by CGH or gene dose array revealed CNAs. The most common recurrent deletions and breakpoints were 14q22-32 (33%), 6q25 (16%), 2p12 (10%), 22q11 (10%), 17p11-13 (10%), 7q32-36 (9%), 18q11-13 (7%), 1q32-44 (6%), 8p22-23 (6%) and 7q11 (6%). Trisomy 5 occurred in 15%. In addition, several other recurrent breakpoints were identified. Although a number of genomic imbalances were identified in the HCL samples, the genome appeared remarkably stable.

Global expression profiling in atopic eczema reveals reciprocal expression of inflammatory and lipid genes.

BACKGROUND: Atopic eczema (AE) is a common chronic inflammatory skin disorder. In order to dissect the genetic background several linkage and genetic association studies have been performed. Yet very little is known about specific genes involved in this complex skin disease, and the underlying molecular mechanisms are not fully understood. METHODOLOGY/FINDINGS: We used human DNA microarrays to identify a molecular picture of the programmed responses of the human genome to AE. The transcriptional program was analyzed in skin biopsy samples from lesional and patch-tested skin from AE patients sensitized to Malassezia sympodialis (M. sympodialis), and corresponding biopsies from healthy individuals. The most notable feature of the global gene-expression pattern observed in AE skin was a reciprocal expression of induced inflammatory genes and repressed lipid metabolism genes. The overall transcriptional response in M. sympodialis patch-tested AE skin was similar to the gene-expression signature identified in lesional AE skin. In the constellation of genes differentially expressed in AE skin compared to healthy control skin, we have identified several potential susceptibility genes that may play a critical role in the pathological condition of AE. Many of these genes, including genes with a role in immune responses, lipid homeostasis, and epidermal differentiation, are localized on chromosomal regions previously linked to AE. CONCLUSIONS/SIGNIFICANCE: Through genome-wide expression profiling, we were able to discover a distinct reciprocal expression pattern of induced inflammatory genes and repressed lipid metabolism genes in skin from AE patients. We found a significant enrichment of differentially expressed genes in AE with cytobands associated to the disease, and furthermore new chromosomal regions were found that could potentially guide future region-specific linkage mapping in AE. The full data set is available at http://microarray-pubs.stanford.edu/eczema.

Parallels between global transcriptional programs of polarizing Caco-2 intestinal epithelial cells in vitro and gene expression programs in normal colon and colon cancer.

Posttranslational mechanisms are implicated in the development of epithelial cell polarity, but little is known about the patterns of gene expression and transcriptional regulation during this process. We characterized temporal patterns of gene expression during cell-cell adhesion-initiated polarization of cultured human Caco-2 cells, which develop structural and functional polarity resembling enterocytes in vivo. A distinctive switch in gene expression patterns occurred upon formation of cell-cell contacts. Comparison to gene expression patterns in normal human colon and colon tumors revealed that the pattern in proliferating, nonpolarized Caco-2 cells paralleled patterns seen in human colon cancer in vivo, including expression of genes involved in cell proliferation. The pattern switched in polarized Caco-2 cells to one more closely resembling that in normal colon tissue, indicating that regulation of transcription underlying Caco-2 cell polarization is similar to that during enterocyte differentiation in vivo. Surprisingly, the temporal program of gene expression in polarizing Caco-2 cells involved changes in signaling pathways (e.g., Wnt, Hh, BMP, FGF) in patterns similar to those during migration and differentiation of intestinal epithelial cells in vivo, despite the absence of morphogen gradients and interactions with stromal cells characteristic of enterocyte differentiation in situ. The full data set is available at http://microarray-pubs.stanford.edu/CACO2.

Transcriptional modulation of genes encoding structural characteristics of differentiating enterocytes during development of a polarized epithelium in vitro.

Although there is considerable evidence implicating posttranslational mechanisms in the development of epithelial cell polarity, little is known about the patterns of gene expression and transcriptional regulation during this process. We characterized the temporal program of gene expression during cell-cell adhesion-initiated polarization of human Caco-2 cells in tissue culture, which develop structural and functional polarity similar to that of enterocytes in vivo. A distinctive switch in gene expression patterns occurred upon formation of cell-cell contacts between neighboring cells. Expression of genes involved in cell proliferation was down-regulated concomitant with induction of genes necessary for functional specialization of polarized epithelial cells. Transcriptional up-regulation of these latter genes correlated with formation of important structural and functional features in enterocyte differentiation and establishment of structural and functional cell polarity; components of the apical microvilli were induced as the brush border formed during polarization; as barrier function was established, expression of tight junction transmembrane proteins peaked; transcripts encoding components of the apical, but not the basal-lateral trafficking machinery were increased during polarization. Coordinated expression of genes encoding components of functional cell structures were often observed indicating temporal control of expression and assembly of multiprotein complexes.