Author Correction: Progression of whole-blood transcriptional signatures from interferon-induced to neutrophil-associated patterns in severe influenza.

In the version of this article initially published, a source of funding was not included in the Acknowledgements section. That section should include the following: P.J.M.O. was supported by EU FP7 PREPARE project 602525. The error has been corrected in the HTML and PDF version of the article.

Progression of whole-blood transcriptional signatures from interferon-induced to neutrophil-associated patterns in severe influenza.

Transcriptional profiles and host-response biomarkers are used increasingly to investigate the severity, subtype and pathogenesis of disease. We now describe whole-blood mRNA signatures and concentrations of local and systemic immunological mediators in 131 adults hospitalized with influenza, from whom extensive clinical and investigational data were obtained by MOSAIC investigators. Signatures reflective of interferon-related antiviral pathways were common up to day 4 of symptoms in patients who did not require mechanical ventilator support; in those who needed mechanical ventilation, an inflammatory, activated-neutrophil and cell-stress or death ('bacterial') pattern was seen, even early in disease. Identifiable bacterial co-infection was not necessary for this 'bacterial' signature but was able to enhance its development while attenuating the early 'viral' signature. Our findings emphasize the importance of timing and severity in the interpretation of host responses to acute viral infection and identify specific patterns of immune-system activation that might enable the development of novel diagnostic and therapeutic tools for severe influenza.

M ycobacterium avium complex genomics and transmission in a London hospital.

BACKGROUND: Non-tuberculous mycobacteria (NTM) are environmental microorganisms and opportunistic pathogens in individuals with pre-existing lung conditions such as cystic fibrosis (CF) and non-CF bronchiectasis. While recent studies of Mycobacterium abscessus have identified transmission within single CF centres as well as nationally and globally, transmission of other NTM species is less well studied. METHODS: To investigate the potential for transmission of the Mycobacterium avium complex (MAC) we sequenced 996 isolates from 354 CF and non-CF patients at the Royal Brompton Hospital (London, UK; collected 2013-2016) and analysed them in a global context. Epidemiological links were identified from patient records. Previously published genomes were used to characterise global population structures. RESULTS: We identified putative transmission clusters in three MAC species, although few epidemiological links could be identified. For M. avium, lineages were largely limited to single countries, while for Mycobacterium chimaera, global transmission clusters previously associated with heater-cooler units (HCUs) were found. However, the immediate ancestor of the lineage causing the major HCU-associated outbreak was a lineage already circulating in patients. CONCLUSIONS: CF and non-CF patients shared transmission chains, although the lack of epidemiological links suggested that most transmission is indirect and may involve environmental intermediates or asymptomatic carriage in the wider population.

The Genetics and Genomics of Asthma.

Asthma is a common, clinically heterogeneous disease with strong evidence of heritability. Progress in defining the genetic underpinnings of asthma, however, has been slow and hampered by issues of inconsistency. Recent advances in the tools available for analysis-assaying transcription, sequence variation, and epigenetic marks on a genome-wide scale-have substantially altered this landscape. Applications of such approaches are consistent with heterogeneity at the level of causation and specify patterns of commonality with a wide range of alternative disease traits. Looking beyond the individual as the unit of study, advances in technology have also fostered comprehensive analysis of the human microbiome and its varied roles in health and disease. In this article, we consider the implications of these technological advances for our current understanding of the genetics and genomics of asthma.

A large-scale genome-wide association analysis of lung function in the Chinese population identifies novel loci and highlights shared genetic aetiology with obesity.

BACKGROUND: Lung function is a heritable complex phenotype with obesity being one of its important risk factors. However, knowledge of their shared genetic basis is limited. Most genome-wide association studies (GWASs) for lung function have been based on European populations, limiting the generalisability across populations. Large-scale lung function GWASs in other populations are lacking. METHODS: We included 100 285 subjects from the China Kadoorie Biobank (CKB). To identify novel loci for lung function, single-trait GWAS analyses were performed on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC in the CKB. We then performed genome-wide cross-trait analysis between lung function and obesity traits (body mass index (BMI), BMI-adjusted waist-to-hip ratio and BMI-adjusted waist circumference) to investigate the shared genetic effects in the CKB. Finally, polygenic risk scores (PRSs) of lung function were developed in the CKB and their interaction with BMI's association on lung function were examined. We also conducted cross-trait analysis in parallel with the CKB using up to 457 756 subjects from the UK Biobank (UKB) for replication and investigation of ancestry-specific effects. RESULTS: We identified nine genome-wide significant novel loci for FEV1, six for FVC and three for FEV1/FVC in the CKB. FEV1 and FVC showed significant negative genetic correlation with obesity traits in both the CKB and UKB. Genetic loci shared between lung function and obesity traits highlighted important biological pathways, including cell proliferation, embryo, skeletal and tissue development, and regulation of gene expression. Mendelian randomisation analysis suggested significant negative causal effects of BMI on FEV1 and on FVC in both the CKB and UKB. Lung function PRSs significantly modified the effect of change in BMI on change in lung function during an average follow-up of 8 years. CONCLUSION: This large-scale GWAS of lung function identified novel loci and shared genetic aetiology between lung function and obesity. Change in BMI might affect change in lung function differently according to a subject's polygenic background. These findings may open new avenues for the development of molecular-targeted therapies for obesity and lung function improvement.

An epigenome-wide association study of total serum immunoglobulin E concentration.

Immunoglobulin E (IgE) is a central mediator of allergic (atopic) inflammation. Therapies directed against IgE can alleviate hay fever and allergic asthma. Genetic association studies have not yet identified novel therapeutic targets or pathways underlying IgE regulation. We therefore surveyed epigenetic associations between serum IgE concentrations and methylation at loci concentrated in CpG islands genome wide in 95 nuclear pedigrees, using DNA from peripheral blood leukocytes. We validated positive results in additional families and in subjects from the general population. Here we show replicated associations--with a meta-analysis false discovery rate less than 10(-4)--between IgE and low methylation at 36 loci. Genes annotated to these loci encode known eosinophil products, and also implicate phospholipid inflammatory mediators, specific transcription factors and mitochondrial proteins. We confirmed that methylation at these loci differed significantly in isolated eosinophils from subjects with and without asthma and high IgE levels. The top three loci accounted for 13% of IgE variation in the primary subject panel, explaining the tenfold higher variance found compared with that derived from large single-nucleotide polymorphism genome-wide association studies. This study identifies novel therapeutic targets and biomarkers for patient stratification for allergic diseases.

ORMDL3 regulates poly I:C induced inflammatory responses in airway epithelial cells.

BACKGROUND: Oroscomucoid 3 (ORMDL3) has been linked to susceptibility of childhood asthma and respiratory viral infection. Polyinosinic-polycytidylic acid (poly I:C) is a synthetic analog of viral double-stranded RNA, a toll-like receptor 3 (TLR3) ligand and mimic of viral infection. METHODS: To investigate the functional role of ORMDL3 in the poly I:C-induced inflammatory response in airway epithelial cells, ORMDL3 knockdown and over-expression models were established in human A549 epithelial cells and primary normal human bronchial epithelial (NHBE) cells. The cells were stimulated with poly I:C or the Th17 cytokine IL-17A. IL-6 and IL-8 levels in supernatants,  mRNA levels of genes in the TLR3 pathway and inflammatory response from cell pellets were measured. ORMDL3 knockdown models in A549 and BEAS-2B epithelial cells were then infected with live human rhinovirus (HRV16) followed by IL-6 and IL-8 measurement. RESULTS: ORMDL3 knockdown and over-expression had little influence on the transcript levels of TLR3 in airway epithelial cells. Time course studies showed that ORMDL3-deficient A549 and NHBE cells had an attenuated IL-6 and IL-8 response to poly I:C stimulation. A549 and NHBE cells over-expressing ORMDL3 released relatively more IL-6 and IL-8 following poly I:C stimulation. IL-17A exhibited a similar inflammatory response in ORMDL3 knockdown and over-expressing cells, but co-stimulation of poly I:C and IL-17A did not significantly enhance the IL-6 and IL-8 response. Transcript abundance of IFNB following poly I:C stimulation was not significantly altered by ORMDL3 knockdown or over-expression. Dampening of the IL-6 response by ORMDL3 knockdown was confirmed in HRV16 infected BEAS-2B and A549 cells. CONCLUSIONS: ORMDL3 regulates the viral inflammatory response in airway epithelial cells via mechanisms independent of the TLR3 pathway.

Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank.

BACKGROUND: Clinical and epidemiologic studies have shown that obesity is associated with asthma and that these associations differ by asthma subtype. Little is known about the shared genetic components between obesity and asthma. OBJECTIVE: We sought to identify shared genetic associations between obesity-related traits and asthma subtypes in adults. METHODS: A cross-trait genome-wide association study (GWAS) was performed using 457,822 subjects of European ancestry from the UK Biobank. Experimental evidence to support the role of genes significantly associated with both obesity-related traits and asthma through a GWAS was sought by using results from obese versus lean mouse RNA sequencing and RT-PCR experiments. RESULTS: We found a substantial positive genetic correlation between body mass index and later-onset asthma defined by asthma age of onset at 16 years or greater (Rg = 0.25, P = 9.56 × 10-22). Mendelian randomization analysis provided strong evidence in support of body mass index causally increasing asthma risk. Cross-trait meta-analysis identified 34 shared loci among 3 obesity-related traits and 2 asthma subtypes. GWAS functional analyses identified potential causal relationships between the shared loci and Genotype-Tissue Expression (GTEx) quantitative trait loci and shared immune- and cell differentiation-related pathways between obesity and asthma. Finally, RNA sequencing data from lungs of obese versus control mice found that 2 genes (acyl-coenzyme A oxidase-like [ACOXL] and myosin light chain 6 [MYL6]) from the cross-trait meta-analysis were differentially expressed, and these findings were validated by using RT-PCR in an independent set of mice. CONCLUSIONS: Our work identified shared genetic components between obesity-related traits and specific asthma subtypes, reinforcing the hypothesis that obesity causally increases the risk of asthma and identifying molecular pathways that might underlie both obesity and asthma.

Presence of pleomorphic features but not growth patterns improves prognostic stratification of epithelioid malignant pleural mesothelioma by 2-tier nuclear grade.

AIMS: Nuclear grade has been recently validated as a powerful prognostic tool in epithelioid malignant pleural mesothelioma (E-MPM). In other studies histological parameters including pleomorphic features and growth patterns were also shown to exert prognostic impact. The primary aims of our study are (i) externally validate the prognostic role of pleomorphic features in E-MPM and (ii) investigate if evaluating growth pattern in addition to 2-tier nuclear grade improves prognostication. METHODS AND RESULTS: 614 consecutive cases of E-MPM from our institution over a period of 15 years were retrospectively reviewed, of which 51 showed pleomorphic features. E-MPM with pleomorphic features showed significantly worse overall survival compared to those without (5.4 versus 14.7 months). Tumours with predominantly micropapillary pattern showed the worst survival (6.2 months) followed by solid (10.5 months), microcystic (15.3 months), discohesive (16.1 months), trabecular (17.6 months) and tubulo-papillary (18.6 months). Sub-classification of growth patterns into high grade (solid, micropapillary) and low grade (all others) led to good separation of overall survival (10.5 versus 18.0 months) but did not predict survival independent of 2-tier nuclear grade. A composite score comprised of growth pattern and 2-tier nuclear grade did not improve prognostication compared with nuclear grade alone. Intra-tumoural heterogeneity in growth patterns is ubiquitous. CONCLUSIONS: Our findings support the incorporation of E-MPM with pleomorphic features in the epithelioid subtype as a highly aggressive variant distinct from 2-tier nuclear grade. E-MPM demonstrates extensive heterogeneity in growth pattern but its evaluation does not offer additional prognostic utility to 2-tier nuclear grade.

The ORMDL3 Asthma Gene Regulates ICAM1 and Has Multiple Effects on Cellular Inflammation.

RATIONALE: Polymorphisms on chromosome 17q21 confer the major genetic susceptibility to childhood-onset asthma. Risk alleles positively correlate with ORMDL3 (orosomucoid-like 3) expression. The locus influences disease severity and the frequency of human rhinovirus (HRV)-initiated exacerbations. ORMDL3 is known to regulate sphingolipid synthesis by binding serine palmitoyltransferase, but its role in inflammation is incompletely understood. OBJECTIVES: To investigate the role of ORMDL3 in cellular inflammation. METHODS: We modeled a time series of IL1B-induced inflammation in A549 cells, using cytokine production as outputs and testing effects of ORMDL3 siRNA knockdown, ORMDL3 overexpression, and the serine palmitoyltransferase inhibitor myriocin. We replicated selected findings in normal human bronchial epithelial cells. Cytokine and metabolite levels were analyzed by analysis of variance. Transcript abundances were analyzed by group means parameterization, controlling the false discovery rate below 0.05. MEASUREMENTS AND MAIN RESULTS: Silencing ORMDL3 led to steroid-independent reduction of IL6 and IL8 release and reduced endoplasmic reticulum stress after IL1B stimulation. Overexpression and myriocin conversely augmented cytokine release. Knockdown reduced expression of genes regulating host-pathogen interactions, stress responses, and ubiquitination: in particular, ORMDL3 knockdown strongly reduced expression of the HRV receptor ICAM1. Silencing led to changes in levels of transcripts and metabolites integral to glycolysis. Increased levels of ceramides and the immune mediator sphingosine-1-phosphate were also observed. CONCLUSIONS: The results show ORMDL3 has pleiotropic effects during cellular inflammation, consistent with its substantial genetic influence on childhood asthma. Actions on ICAM1 provide a mechanism for the locus to confer susceptibility to HRV-induced asthma.

Genome-wide interaction study of early-life smoking exposure on time-to-asthma onset in childhood.

BACKGROUND: Asthma, a heterogeneous disease with variable age of onset, results from the interplay between genetic and environmental factors. Early-life tobacco smoke (ELTS) exposure is a major asthma risk factor. Only a few genetic loci have been reported to interact with ELTS exposure in asthma. OBJECTIVE: Our aim was to identify new loci interacting with ELTS exposure on time-to-asthma onset (TAO) in childhood. METHODS: We conducted genome-wide interaction analyses of ELTS exposure on time-to-asthma onset in childhood in five European-ancestry studies (totalling 8273 subjects) using Cox proportional-hazard model. The results of all five genome-wide analyses were meta-analysed. RESULTS: The 13q21 locus showed genome-wide significant interaction with ELTS exposure (P = 4.3 × 10-8 for rs7334050 within KLHL1 with consistent results across the five studies). Suggestive interactions (P < 5 × 10-6 ) were found at three other loci: 20p12 (rs13037508 within MACROD2; P = 4.9 × 10-7 ), 14q22 (rs7493885 near NIN; P = 2.9 × 10-6 ) and 2p22 (rs232542 near CYP1B1; P = 4.1 × 10-6 ). Functional annotations and the literature showed that the lead SNPs at these four loci influence DNA methylation in the blood and are located nearby CpG sites reported to be associated with exposure to tobacco smoke components, which strongly support our findings. CONCLUSIONS AND CLINICAL RELEVANCE: We identified novel candidate genes interacting with ELTS exposure on time-to-asthma onset in childhood. These genes have plausible biological relevance related to tobacco smoke exposure. Further epigenetic and functional studies are needed to confirm these findings and to shed light on the underlying mechanisms.

A novel role for ciliary function in atopy: ADGRV1 and DNAH5 interactions.

BACKGROUND: Atopy, an endotype underlying allergic diseases, has a substantial genetic component. OBJECTIVE: Our goal was to identify novel genes associated with atopy in asthma-ascertained families. METHODS: We implemented a 3-step analysis strategy in 3 data sets: the Epidemiological Study on the Genetics and Environment of Asthma (EGEA) data set (1660 subjects), the Saguenay-Lac-Saint-Jean study data set (1138 subjects), and the Medical Research Council (MRC) data set (446 subjects). This strategy included a single nucleotide polymorphism (SNP) genome-wide association study (GWAS), the selection of related gene pairs based on statistical filtering of GWAS results, and text-mining filtering using Gene Relationships Across Implicated Loci and SNP-SNP interaction analysis of selected gene pairs. RESULTS: We identified the 5q14 locus, harboring the adhesion G protein-coupled receptor V1 (ADGRV1) gene, which showed genome-wide significant association with atopy (rs4916831, meta-analysis P value = 6.8 × 10-9). Statistical filtering of GWAS results followed by text-mining filtering revealed relationships between ADGRV1 and 3 genes showing suggestive association with atopy (P ≤ 10-4). SNP-SNP interaction analysis between ADGRV1 and these 3 genes showed significant interaction between ADGRV1 rs17554723 and 2 correlated SNPs (rs2134256 and rs1354187) within the dynein axonemal heavy chain 5 (DNAH5) gene (Pmeta-int = 3.6 × 10-5 and 6.1 × 10-5, which met the multiple-testing corrected threshold of 7.3 × 10-5). Further conditional analysis indicated that rs2134256 alone accounted for the interaction signal with rs17554723. CONCLUSION: Because both DNAH5 and ADGRV1 contribute to ciliary function, this study suggests that ciliary dysfunction might represent a novel mechanism underlying atopy. Combining GWAS and epistasis analysis driven by statistical and knowledge-based evidence represents a promising approach for identifying new genes involved in complex traits.

Whole-Blood Gene Expression in Pulmonary Nontuberculous Mycobacterial Infection.

The factors predisposing toward the development of pulmonary nontuberculous mycobacterial (pNTM) disease and influencing disease progression remain unclear. Impaired immune responses have been reported in individuals with pNTM disease, but data are limited and inconsistent. In this study, we sought to use gene expression profiling to examine the host response to pNTM disease. Microarray analysis of whole-blood gene expression was performed on 25 subjects with pNTM disease and 27 uninfected control subjects with respiratory disease. Gene expression results were compared with phenotypic variables and survival data. Compared with uninfected control subjects, pNTM disease was associated with downregulation of 213 transcripts enriched for terms related to T cell signaling, including IFNG. Reduced IFNG expression was associated with more severe computed tomography changes and impaired lung function. Mortality was associated with the expression of transcripts related to the innate immune response and inflammation, whereas transcripts related to T and B cell function were associated with improved survival. These findings suggest that pNTM disease is associated with an aberrant immune response, which may reflect an underlying propensity to infection or result from NTM infection itself. There were important differences in the immune response associated with survival and mortality in pNTM disease.

Gene expression in transformed lymphocytes reveals variation in endomembrane and HLA pathways modifying cystic fibrosis pulmonary phenotypes.

Variation in cystic fibrosis (CF) phenotypes, including lung disease severity, age of onset of persistent Pseudomonas aeruginosa (P. aeruginosa) lung infection, and presence of meconium ileus (MI), has been partially explained by genome-wide association studies (GWASs). It is not expected that GWASs alone are sufficiently powered to uncover all heritable traits associated with CF phenotypic diversity. Therefore, we utilized gene expression association from lymphoblastoid cells lines from 754 p.Phe508del CF-affected homozygous individuals to identify genes and pathways. LPAR6, a G protein coupled receptor, associated with lung disease severity (false discovery rate q value = 0.0006). Additional pathway analyses, utilizing a stringent permutation-based approach, identified unique signals for all three phenotypes. Pathways associated with lung disease severity were annotated in three broad categories: (1) endomembrane function, containing p.Phe508del processing genes, providing evidence of the importance of p.Phe508del processing to explain lung phenotype variation; (2) HLA class I genes, extending previous GWAS findings in the HLA region; and (3) endoplasmic reticulum stress response genes. Expression pathways associated with lung disease were concordant for some endosome and HLA pathways, with pathways identified using GWAS associations from 1,978 CF-affected individuals. Pathways associated with age of onset of persistent P. aeruginosa infection were enriched for HLA class II genes, and those associated with MI were related to oxidative phosphorylation. Formal testing demonstrated that genes showing differential expression associated with lung disease severity were enriched for heritable genetic variation and expression quantitative traits. Gene expression provided a powerful tool to identify unrecognized heritable variation, complementing ongoing GWASs in this rare disease.

Metal worker's lung: spatial association with Mycobacterium avium.

BACKGROUND: Outbreaks of hypersensitivity pneumonitis (HP) are not uncommon in workplaces where metal working fluid (MWF) is used to facilitate metal turning. Inhalation of microbe-contaminated MWF has been assumed to be the cause, but previous investigations have failed to establish a spatial relationship between a contaminated source and an outbreak. OBJECTIVES: After an outbreak of five cases of HP in a UK factory, we carried out blinded, molecular-based microbiological investigation of MWF samples in order to identify potential links between specific microbial taxa and machines in the outbreak zone. METHODS: Custom-quantitative PCR assays, microscopy and phylogenetic analyses were performed on blinded MWF samples to quantify microbial burden and identify potential aetiological agents of HP in metal workers. MEASUREMENTS AND MAIN RESULTS: MWF from machines fed by a central sump, but not those with an isolated supply, was contaminated by mycobacteria. The factory sump and a single linked machine at the centre of the outbreak zone, known to be the workstation of the index cases, had very high levels of detectable organisms. Phylogenetic placement of mycobacterial taxonomic marker genes generated from these samples indicated that the contaminating organisms were closely related to Mycobacterium avium. CONCLUSIONS: We describe, for the first time, a close spatial relationship between the abundance of a mycobacterium-like organism, most probably M. avium, and a localised outbreak of MWF-associated HP. The further development of sequence-based analytic techniques should assist in the prevention of this important occupational disease.

Host-Microbial Interactions in Idiopathic Pulmonary Fibrosis.

RATIONALE: Changes in the respiratory microbiome are associated with disease progression in idiopathic pulmonary fibrosis (IPF). The role of the host response to the respiratory microbiome remains unknown. OBJECTIVES: To explore the host-microbial interactions in IPF. METHODS: Sixty patients diagnosed with IPF were prospectively enrolled together with 20 matched control subjects. Subjects underwent bronchoalveolar lavage (BAL), and peripheral whole blood was collected into PAXgene tubes for all subjects at baseline. For subjects with IPF, additional samples were taken at 1, 3, and 6 months and (if alive) 1 year. Gene expression profiles were generated using Affymetrix Human Gene 1.1 ST arrays. MEASUREMENTS AND MAIN RESULTS: By network analysis of gene expression data, we identified two gene modules that strongly associated with a diagnosis of IPF, BAL bacterial burden (determined by 16S quantitative polymerase chain reaction), and specific microbial operational taxonomic units, as well as with lavage and peripheral blood neutrophilia. Genes within these modules that are involved in the host defense response include NLRC4, PGLYRP1, MMP9, and DEFA4. The modules also contain two genes encoding specific antimicrobial peptides (SLPI and CAMP). Many of these particular transcripts were associated with survival and showed longitudinal overexpression in subjects experiencing disease progression, further strengthening the relationship of the transcripts with disease. CONCLUSIONS: Integrated analysis of the host transcriptome and microbial signatures demonstrated an apparent host response to the presence of an altered or more abundant microbiome. These responses remained elevated in longitudinal follow-up, suggesting that the bacterial communities of the lower airways may act as persistent stimuli for repetitive alveolar injury in IPF.

Pulmonary ORMDL3 is critical for induction of Alternaria-induced allergic airways disease.

BACKGROUND: Genome-wide association studies have identified the ORM (yeast)-like protein isoform 3 (ORMDL3) gene locus on human chromosome 17q to be a highly significant risk factor for childhood-onset asthma. OBJECTIVE: We sought to investigate in vivo the functional role of ORMDL3 in disease inception. METHODS: An Ormdl3-deficient mouse was generated and the role of ORMDL3 in the generation of allergic airways disease to the fungal aeroallergen Alternaria alternata was determined. An adeno-associated viral vector was also used to reconstitute ORMDL3 expression in airway epithelial cells of Ormdl3 knockout mice. RESULTS: Ormdl3 knockout mice were found to be protected from developing allergic airways disease and showed a marked decrease in pathophysiology, including lung function and airway eosinophilia induced by Alternaria. Alternaria is a potent inducer of cellular stress and the unfolded protein response, and ORMDL3 was found to play a critical role in driving the activating transcription factor 6-mediated arm of this response through Xbp1 and downstream activation of the endoplasmic reticulum-associated degradation pathway. In addition, ORMDL3 mediated uric acid release, another marker of cellular stress. In the knockout mice, reconstitution of Ormdl3 transcript levels specifically in the bronchial epithelium resulted in reinstatement of susceptibility to fungal allergen-induced allergic airways disease. CONCLUSIONS: This study demonstrates that ORMDL3, an asthma susceptibility gene identified by genome-wide association studies, contributes to key pathways that promote changes in airway physiology during allergic immune responses.

Bacterial microbiota of the upper respiratory tract and childhood asthma.

BACKGROUND: Patients with asthma and healthy controls differ in bacterial colonization of the respiratory tract. The upper airways have been shown to reflect colonization of the lower airways, the actual site of inflammation in asthma, which is hardly accessible in population studies. OBJECTIVE: We sought to characterize the bacterial communities at 2 sites of the upper respiratory tract obtained from children from a rural area and to relate these to asthma. METHODS: The microbiota of 327 throat and 68 nasal samples from school-age farm and nonfarm children were analyzed by 454-pyrosequencing of the bacterial 16S ribosomal RNA gene. RESULTS: Alterations in nasal microbiota but not of throat microbiota were associated with asthma. Children with asthma had lower α- and ß-diversity of the nasal microbiota as compared with healthy control children. Furthermore, asthma presence was positively associated with a specific operational taxonomic unit from the genus Moraxella in children not exposed to farming, whereas in farm children Moraxella colonization was unrelated to asthma. In nonfarm children, Moraxella colonization explained the association between bacterial diversity and asthma to a large extent. CONCLUSIONS: Asthma was mainly associated with an altered nasal microbiota characterized by lower diversity and Moraxella abundance. Children living on farms might not be susceptible to the disadvantageous effect of Moraxella. Prospective studies may clarify whether Moraxella outgrowth is a cause or a consequence of loss in diversity.

An epigenome-wide association study of total serum IgE in Hispanic children.

BACKGROUND: Total IgE is a therapeutic target in patients with allergic diseases. DNA methylation in white blood cells (WBCs) was associated with total IgE levels in an epigenome-wide association study of white subjects. Whether DNA methylation of eosinophils explains these findings is insufficiently understood. METHODS: We tested for association between genome-wide DNA methylation in WBCs and total IgE levels in 2 studies of Hispanic children: the Puerto Rico Genetics of Asthma and Lifestyle Study (PR-GOAL; n = 306) and the Genes-environments and Admixture in Latino Americans (GALA II) study (n = 573). Whole-genome methylation of DNA from WBCs was measured by using the Illumina Infinium HumanMethylation450 BeadChip. Total IgE levels were measured by using the UniCAP 100 system. In PR-GOAL WBC types (ie, neutrophils, eosinophils, basophils, lymphocytes, and monocytes) in peripheral blood were measured by using Coulter Counter techniques. In the GALA II study WBC types were imputed. Multivariable linear regression was used for the analysis of DNA methylation and total IgE levels, which was first conducted separately for each cohort, and then results from the 2 cohorts were combined in a meta-analysis. RESULTS: CpG sites in multiple genes, including novel findings and results previously reported in white subjects, were significantly associated with total IgE levels. However, adjustment for WBC types resulted in markedly fewer significant sites. Top findings from this adjusted meta-analysis were in the genes ZFPM1 (P = 1.5 × 10-12), ACOT7 (P = 2.5 × 10-11), and MND1 (P = 1.4 × 10-9). CONCLUSIONS: In an epigenome-wide association study adjusted for WBC types (including eosinophils), methylation changes in genes enriched in pathways relevant to asthma and immune responses were associated with total IgE levels among Hispanic children.

Fast and robust adjustment of cell mixtures in epigenome-wide association studies with SmartSVA.

BACKGROUND: One problem that plagues epigenome-wide association studies is the potential confounding due to cell mixtures when purified target cells are not available. Reference-free adjustment of cell mixtures has become increasingly popular due to its flexibility and simplicity. However, existing methods are still not optimal: increased false positive rates and reduced statistical power have been observed in many scenarios. METHODS: We develop SmartSVA, an optimized surrogate variable analysis (SVA) method, for fast and robust reference-free adjustment of cell mixtures. SmartSVA corrects the limitation of traditional SVA under highly confounded scenarios by imposing an explicit convergence criterion and improves the computational efficiency for large datasets. RESULTS: Compared to traditional SVA, SmartSVA achieves an order-of-magnitude speedup and better false positive control. It protects the signals when capturing the cell mixtures, resulting in significant power increase while controlling for false positives. Through extensive simulations and real data applications, we demonstrate a better performance of SmartSVA than the existing methods. CONCLUSIONS: SmartSVA is a fast and robust method for reference-free adjustment of cell mixtures for epigenome-wide association studies. As a general method, SmartSVA can be applied to other genomic studies to capture unknown sources of variability.