Epigenetic responses to rhinovirus exposure in airway epithelial cells are correlated with key transcriptional pathways in chronic rhinosinusitis.

BACKGROUND: Viruses may drive immune mechanisms responsible for chronic rhinosinusitis with nasal polyposis (CRSwNP), but little is known about the underlying molecular mechanisms. OBJECTIVES: To identify epigenetic and transcriptional responses to a common upper respiratory pathogen, rhinovirus (RV), that are specific to patients with CRSwNP using a primary sinonasal epithelial cell culture model. METHODS: Airway epithelial cells were collected at surgery from patients with CRSwNP (cases) and from controls without sinus disease, cultured, and then exposed to RV or vehicle for 48 h. Differential gene expression and DNA methylation (DNAm) between cases and controls in response to RV were determined using linear mixed models. Weighted gene co-expression analysis (WGCNA) was used to identify (a) co-regulated gene expression and DNAm signatures, and (b) genes, pathways, and regulatory mechanisms specific to CRSwNP. RESULTS: We identified 5585 differential transcriptional and 261 DNAm responses (FDR <0.10) to RV between CRSwNP cases and controls. These differential responses formed three co-expression/co-methylation modules that were related to CRSwNP and three that were related to RV (Bonferroni corrected p < .01). Most (95%) of the differentially methylated CpGs (DMCs) were in modules related to CRSwNP, whereas the differentially expressed genes (DEGs) were more equally distributed between the CRSwNP- and RV-related modules. Genes in the CRSwNP-related modules were enriched in known CRS and/or viral response immune pathways. CONCLUSION: RV activates specific epigenetic programs and correlated transcriptional networks in the sinonasal epithelium of individuals with CRSwNP. These novel observations suggest epigenetic signatures specific to patients with CRSwNP modulate response to viral pathogens at the mucosal environmental interface. Determining how viral response pathways are involved in epithelial inflammation in CRSwNP could lead to therapeutic targets for this burdensome airway disorder.

Fine-mapping studies distinguish genetic risks for childhood- and adult-onset asthma in the HLA region.

BACKGROUND: Genome-wide association studies of asthma have revealed robust associations with variation across the human leukocyte antigen (HLA) complex with independent associations in the HLA class I and class II regions for both childhood-onset asthma (COA) and adult-onset asthma (AOA). However, the specific variants and genes contributing to risk are unknown. METHODS: We used Bayesian approaches to perform genetic fine-mapping for COA and AOA (n=9432 and 21,556, respectively; n=318,167 shared controls) in White British individuals from the UK Biobank and to perform expression quantitative trait locus (eQTL) fine-mapping in immune (lymphoblastoid cell lines, n=398; peripheral blood mononuclear cells, n=132) and airway (nasal epithelial cells, n=188) cells from ethnically diverse individuals. We also examined putatively causal protein coding variation from protein crystal structures and conducted replication studies in independent multi-ethnic cohorts from the UK Biobank (COA n=1686; AOA n=3666; controls n=56,063). RESULTS: Genetic fine-mapping revealed both shared and distinct causal variation between COA and AOA in the class I region but only distinct causal variation in the class II region. Both gene expression levels and amino acid variation contributed to risk. Our results from eQTL fine-mapping and amino acid visualization suggested that the HLA-DQA1*03:01 allele and variation associated with expression of the nonclassical HLA-DQA2 and HLA-DQB2 genes accounted entirely for the most significant association with AOA in GWAS. Our studies also suggested a potentially prominent role for HLA-C protein coding variation in the class I region in COA. We replicated putatively causal variant associations in a multi-ethnic cohort. CONCLUSIONS: We highlight roles for both gene expression and protein coding variation in asthma risk and identified putatively causal variation and genes in the HLA region. A convergence of genomic, transcriptional, and protein coding evidence implicates the HLA-DQA2 and HLA-DQB2 genes and HLA-DQA1*03:01 allele in AOA.

Multi-omics colocalization with genome-wide association studies reveals a context-specific genetic mechanism at a childhood onset asthma risk locus.

BACKGROUND: Genome-wide association studies (GWASs) have identified thousands of variants associated with asthma and other complex diseases. However, the functional effects of most of these variants are unknown. Moreover, GWASs do not provide context-specific information on cell types or environmental factors that affect specific disease risks and outcomes. To address these limitations, we used an upper airway epithelial cell (AEC) culture model to assess transcriptional and epigenetic responses to rhinovirus (RV), an asthma-promoting pathogen, and provide context-specific functional annotations to variants discovered in GWASs of asthma. METHODS: Genome-wide genetic, gene expression, and DNA methylation data in vehicle- and RV-treated upper AECs were collected from 104 individuals who had a diagnosis of airway disease (n=66) or were healthy participants (n=38). We mapped cis expression and methylation quantitative trait loci (cis-eQTLs and cis-meQTLs, respectively) in each treatment condition (RV and vehicle) in AECs from these individuals. A Bayesian test for colocalization between AEC molecular QTLs and adult onset asthma and childhood onset asthma GWAS SNPs, and a multi-ethnic GWAS of asthma, was used to assign the function to variants associated with asthma. We used Mendelian randomization to demonstrate DNA methylation effects on gene expression at asthma colocalized loci. RESULTS: Asthma and allergic disease-associated GWAS SNPs were specifically enriched among molecular QTLs in AECs, but not in GWASs from non-immune diseases, and in AEC eQTLs, but not among eQTLs from other tissues. Colocalization analyses of AEC QTLs with asthma GWAS variants revealed potential molecular mechanisms of asthma, including QTLs at the TSLP locus that were common to both the RV and vehicle treatments and to both childhood onset and adult onset asthma, as well as QTLs at the 17q12-21 asthma locus that were specific to RV exposure and childhood onset asthma, consistent with clinical and epidemiological studies of these loci. CONCLUSIONS: This study provides evidence of functional effects for asthma risk variants in AECs and insight into RV-mediated transcriptional and epigenetic response mechanisms that modulate genetic effects in the airway and risk for asthma.

Chromosome 17q12-21 Variants Are Associated with Multiple Wheezing Phenotypes in Childhood.

Rationale: Birth cohort studies have identified several temporal patterns of wheezing, only some of which are associated with asthma. Whether 17q12-21 genetic variants, which are closely associated with asthma, are also associated with childhood wheezing phenotypes remains poorly explored.Objectives: To determine whether wheezing phenotypes, defined by latent class analysis (LCA), are associated with nine 17q12-21 SNPs and if so, whether these relationships differ by race/ancestry.Methods: Data from seven U.S. birth cohorts (n = 3,786) from the CREW (Children's Respiratory Research and Environment Workgroup) were harmonized to represent whether subjects wheezed in each year of life from birth until age 11 years. LCA was then performed to identify wheeze phenotypes. Genetic associations between SNPs and wheeze phenotypes were assessed separately in European American (EA) (n = 1,308) and, for the first time, in African American (AA) (n = 620) children.Measurements and Main Results: The LCA best supported four latent classes of wheeze: infrequent, transient, late-onset, and persistent. Odds of belonging to any of the three wheezing classes (vs. infrequent) increased with the risk alleles for multiple SNPs in EA children. Only one SNP, rs2305480, showed increased odds of belonging to any wheezing class in both AA and EA children.Conclusions: These results indicate that 17q12-21 is a "wheezing locus," and this association may reflect an early life susceptibility to respiratory viruses common to all wheezing children. Which children will have their symptoms remit or reoccur during childhood may be independent of the influence of rs2305480.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma.

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

Epigenetic landscape links upper airway microbiota in infancy with allergic rhinitis at 6 years of age.

BACKGROUND: The upper airways present a barrier to inhaled allergens and microbes, which alter immune responses and subsequent risk for diseases, such as allergic rhinitis (AR). OBJECTIVE: We tested the hypothesis that early-life microbial exposures leave a lasting signature in DNA methylation that ultimately influences the development of AR in children. METHODS: We studied upper airway microbiota at 1 week, 1 month, and 3 months of life, and measured DNA methylation and gene expression profiles in upper airway mucosal cells and assessed AR at age 6 years in children in the Copenhagen Prospective Studies on Asthma in Childhood birth cohort. RESULTS: We identified 956 AR-associated differentially methylated CpGs in upper airway mucosal cells at age 6 years, 792 of which formed 3 modules of correlated differentially methylated CpGs. The eigenvector of 1 module was correlated with the expression of genes enriched for lysosome and bacterial invasion of epithelial cell pathways. Early-life microbial diversity was lower at 1 week (richness P = .0079) in children with AR at age 6 years, and reduced diversity at 1 week was also correlated with the same module's eigenvector (ρ = -0.25; P = 3.3 × 10-5). We show that the effect of microbiota richness at 1 week on risk for AR at age 6 years was mediated in part by the epigenetic signature of this module. CONCLUSIONS: Our results suggest that upper airway microbial composition in infancy contributes to the development of AR during childhood, and this trajectory is mediated, at least in part, through altered DNA methylation patterns in upper airway mucosal cells.

Expression quantitative trait locus fine mapping of the 17q12-21 asthma locus in African American children: a genetic association and gene expression study.

BACKGROUND: African ancestry is associated with a higher prevalence and greater severity of asthma than European ancestries, yet genetic studies of the most common locus associated with childhood-onset asthma, 17q12-21, in African Americans have been inconclusive. The aim of this study was to leverage both the phenotyping of the Children's Respiratory and Environmental Workgroup (CREW) birth cohort consortium, and the reduced linkage disequilibrium in African Americans, to fine map the 17q12-21 locus. METHODS: We first did a genetic association study and meta-analysis using 17q12-21 tag single-nucleotide polymorphisms (SNPs) for childhood-onset asthma in 1613 European American and 870 African American children from the CREW consortium. Nine tag SNPs were selected based on linkage disequilibrium patterns at 17q12-21 and their association with asthma, considering the effect allele under an additive model (0, 1, or 2 effect alleles). Results were meta-analysed with publicly available summary data from the EVE consortium (on 4303 European American and 3034 African American individuals) for seven of the nine SNPs of interest. Subsequently, we tested for expression quantitative trait loci (eQTLs) among the SNPs associated with childhood-onset asthma and the expression of 17q12-21 genes in resting peripheral blood mononuclear cells (PBMCs) from 85 African American CREW children and in upper airway epithelial cells from 246 African American CREW children; and in lower airway epithelial cells from 44 European American and 72 African American adults from a case-control study of asthma genetic risk in Chicago (IL, USA). FINDINGS: 17q12-21 SNPs were broadly associated with asthma in European Americans. Only two SNPs (rs2305480 in gasdermin-B [GSDMB] and rs8076131 in ORMDL sphingolipid biosynthesis regulator 3 [ORMDL3]) were associated with asthma in African Americans, at a Bonferroni-corrected threshold of p<0·0055 (for rs2305480_G, odds ratio [OR] 1·36 [95% CI 1·12-1·65], p=0·0014; and for rs8076131_A, OR 1·37 [1·13-1·67], p=0·0010). In upper airway epithelial cells from African American children, genotype at rs2305480 was the most significant eQTL for GSDMB (eQTL effect size [ß] 1·35 [95% CI 1·25-1·46], p<0·0001), and to a lesser extent showed an eQTL effect for post-GPI attachment to proteins phospholipase 3 (ß 1·15 [1·08-1·22], p<0·0001). No SNPs were eQTLs for ORMDL3. By contrast, in PBMCs, the five core SNPs were associated only with expression of GSDMB and ORMDL3. Genotype at rs12936231 (in zona pellucida binding protein 2) showed the strongest associations across both genes (for GSDMB, eQTLß 1·24 [1·15-1·32], p<0·0001; and for ORMDL3 (ß 1·19 [1·12-1·24], p<0·0001). The eQTL effects of rs2305480 on GSDMB expression were replicated in lower airway cells from African American adults (ß 1·29 [1·15-1·44], p<0·0001). INTERPRETATION: Our study suggests that SNPs regulating GSDMB expression in airway epithelial cells have a major role in childhood-onset asthma, whereas SNPs regulating the expression levels of 17q12-21 genes in resting blood cells are not central to asthma risk. Our genetic and gene expression data in African Americans and European Americans indicated GSDMB to be the leading candidate gene at this important asthma locus. FUNDING: National Institutes of Health, Office of the Director.

Association of HLA-DRB1∗09:01 with tIgE levels among African-ancestry individuals with asthma.

BACKGROUND: Asthma is a complex chronic inflammatory disease of the airways. Association studies between HLA and asthma were first reported in the 1970s, and yet, the precise role of HLA alleles in asthma is not fully understood. Numerous genome-wide association studies were recently conducted on asthma, but were always limited to simple genetic markers (single nucleotide polymorphisms) and not complex HLA gene polymorphisms (alleles/haplotypes), therefore not capturing the biological relevance of this complex locus for asthma pathogenesis. OBJECTIVE: To run the first HLA-centric association study with asthma and specific asthma-related phenotypes in a large cohort of African-ancestry individuals. METHODS: We collected high-density genomics data for the Consortium on Asthma among African-ancestry Populations in the Americas (N = 4993) participants. Using computer-intensive machine-learning attribute bagging methods to infer HLA alleles, and Easy-HLA to infer HLA 5-gene haplotypes, we conducted a high-throughput HLA-centric association study of asthma susceptibility and total serum IgE (tIgE) levels in subjects with and without asthma. RESULTS: Among the 1607 individuals with asthma, 972 had available tIgE levels, with a mean tIgE level of 198.7 IU/mL. We could not identify any association with asthma susceptibility. However, we showed that HLA-DRB1∗09:01 was associated with increased tIgE levels (P = 8.5 × 10-4; weighted effect size, 0.51 [0.15-0.87]). CONCLUSIONS: We identified for the first time an HLA allele associated with tIgE levels in African-ancestry individuals with asthma. Our report emphasizes that by leveraging powerful computational machine-learning methods, specific/extreme phenotypes, and population diversity, we can explore HLA gene polymorphisms in depth and reveal the full extent of complex disease associations.

Author Correction: Association study in African-admixed populations across the Americas recapitulates asthma risk loci in non-African populations.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Shared and distinct genetic risk factors for childhood-onset and adult-onset asthma: genome-wide and transcriptome-wide studies.

BACKGROUND: Childhood-onset and adult-onset asthma differ with respect to severity and comorbidities. Whether they also differ with respect to genetic risk factors has not been previously investigated in large samples. The goals of this study were to identify shared and distinct genetic risk loci for childhood-onset and adult-onset asthma, and to identify the genes that might mediate the effects of associated variation. METHODS: We did genome-wide and transcriptome-wide studies, using data from the UK Biobank, in individuals with asthma, including adults with childhood-onset asthma (onset before 12 years of age), adults with adult-onset asthma (onset between 26 and 65 years of age), and adults without asthma (controls; aged older than 38 years). We did genome-wide association studies (GWAS) for childhood-onset asthma and adult-onset asthma each compared with shared controls, and for age of asthma onset in all asthma cases, with a genome-wide significance threshold of p<5 × 10-8. Enrichment studies determined the tissues in which genes at GWAS loci were most highly expressed, and PrediXcan, a transcriptome-wide gene-based test, was used to identify candidate risk genes. FINDINGS: Of 376 358 British white individuals from the UK Biobank, we included 37 846 with self-reports of doctor-diagnosed asthma: 9433 adults with childhood-onset asthma; 21 564 adults with adult-onset asthma; and an additional 6849 young adults with asthma with onset between 12 and 25 years of age. For the first and second GWAS analyses, 318 237 individuals older than 38 years without asthma were used as controls. We detected 61 independent asthma loci: 23 were childhood-onset specific, one was adult-onset specific, and 37 were shared. 19 loci were associated with age of asthma onset. The most significant asthma-associated locus was at 17q12 (odds ratio 1·406, 95% CI 1·365-1·448; p=1·45 × 10-111) in the childhood-onset GWAS. Genes at the childhood onset-specific loci were most highly expressed in skin, blood, and small intestine; genes at the adult onset-specific loci were most highly expressed in lung, blood, small intestine, and spleen. PrediXcan identified 113 unique candidate genes at 22 of the 61 GWAS loci. Single-nucleotide polymorphism-based heritability estimates were more than three times larger for childhood-onset asthma (0·327) than for adult-onset disease (0·098). The onset of disease in childhood was associated with additional genes with relatively large effect sizes, with the largest odds ratio observed at the FLG locus at 1q21.3 (1·970, 95% CI 1·823-2·129). INTERPRETATION: Genetic risk factors for adult-onset asthma are largely a subset of the genetic risk for childhood-onset asthma but with overall smaller effects, suggesting a greater role for non-genetic risk factors in adult-onset asthma. Combined with gene expression and tissue enrichment patterns, we suggest that the establishment of disease in children is driven more by dysregulated allergy and epithelial barrier function genes, whereas the cause of adult-onset asthma is more lung-centred and environmentally determined, but with immune-mediated mechanisms driving disease progression in both children and adults. FUNDING: US National Institutes of Health.

Association study in African-admixed populations across the Americas recapitulates asthma risk loci in non-African populations.

Asthma is a complex disease with striking disparities across racial and ethnic groups. Despite its relatively high burden, representation of individuals of African ancestry in asthma genome-wide association studies (GWAS) has been inadequate, and true associations in these underrepresented minority groups have been inconclusive. We report the results of a genome-wide meta-analysis from the Consortium on Asthma among African Ancestry Populations (CAAPA; 7009 asthma cases, 7645 controls). We find strong evidence for association at four previously reported asthma loci whose discovery was driven largely by non-African populations, including the chromosome 17q12-q21 locus and the chr12q13 region, a novel (and not previously replicated) asthma locus recently identified by the Trans-National Asthma Genetic Consortium (TAGC). An additional seven loci reported by TAGC show marginal evidence for association in CAAPA. We also identify two novel loci (8p23 and 8q24) that may be specific to asthma risk in African ancestry populations.

Parent-of-origin effects on quantitative phenotypes in a large Hutterite pedigree.

The impact of the parental origin of associated alleles in GWAS has been largely ignored. Yet sequence variants could affect traits differently depending on whether they are inherited from the mother or the father, as in imprinted regions, where identical inherited DNA sequences can have different effects based on the parental origin. To explore parent-of-origin effects (POEs), we studied 21 quantitative phenotypes in a large Hutterite pedigree to identify variants with single parent (maternal-only or paternal-only) effects, and then variants with opposite parental effects. Here we show that POEs, which can be opposite in direction, are relatively common in humans, have potentially important clinical effects, and will be missed in traditional GWAS. We identified POEs with 11 phenotypes, most of which are risk factors for cardiovascular disease. Many of the loci identified are characteristic of imprinted regions and are associated with the expression of nearby genes.

Author Correction: Assembly of a pan-genome from deep sequencing of 910 humans of African descent.

In the version of this article initially published, the statement "there are no pan-genomes for any other animal or plant species" was incorrect. The statement has been corrected to "there are no reported pan-genomes for any other animal species, to our knowledge." We thank David Edwards for bringing this error to our attention. The error has been corrected in the HTML and PDF versions of the article.

Integrating predicted transcriptome from multiple tissues improves association detection.

Integration of genome-wide association studies (GWAS) and expression quantitative trait loci (eQTL) studies is needed to improve our understanding of the biological mechanisms underlying GWAS hits, and our ability to identify therapeutic targets. Gene-level association methods such as PrediXcan can prioritize candidate targets. However, limited eQTL sample sizes and absence of relevant developmental and disease context restrict our ability to detect associations. Here we propose an efficient statistical method (MultiXcan) that leverages the substantial sharing of eQTLs across tissues and contexts to improve our ability to identify potential target genes. MultiXcan integrates evidence across multiple panels using multivariate regression, which naturally takes into account the correlation structure. We apply our method to simulated and real traits from the UK Biobank and show that, in realistic settings, we can detect a larger set of significantly associated genes than using each panel separately. To improve applicability, we developed a summary result-based extension called S-MultiXcan, which we show yields highly concordant results with the individual level version when LD is well matched. Our multivariate model-based approach allowed us to use the individual level results as a gold standard to calibrate and develop a robust implementation of the summary-based extension. Results from our analysis as well as software and necessary resources to apply our method are publicly available.

An admixture mapping meta-analysis implicates genetic variation at 18q21 with asthma susceptibility in Latinos.

BACKGROUND: Asthma is a common but complex disease with racial/ethnic differences in prevalence, morbidity, and response to therapies. OBJECTIVE: We sought to perform an analysis of genetic ancestry to identify new loci that contribute to asthma susceptibility. METHODS: We leveraged the mixed ancestry of 3902 Latinos and performed an admixture mapping meta-analysis for asthma susceptibility. We replicated associations in an independent study of 3774 Latinos, performed targeted sequencing for fine mapping, and tested for disease correlations with gene expression in the whole blood of more than 500 subjects from 3 racial/ethnic groups. RESULTS: We identified a genome-wide significant admixture mapping peak at 18q21 in Latinos (P = 6.8 × 10-6), where Native American ancestry was associated with increased risk of asthma (odds ratio [OR], 1.20; 95% CI, 1.07-1.34; P = .002) and European ancestry was associated with protection (OR, 0.86; 95% CI, 0.77-0.96; P = .008). Our findings were replicated in an independent childhood asthma study in Latinos (P = 5.3 × 10-3, combined P = 2.6 × 10-7). Fine mapping of 18q21 in 1978 Latinos identified a significant association with multiple variants 5' of SMAD family member 2 (SMAD2) in Mexicans, whereas a single rare variant in the same window was the top association in Puerto Ricans. Low versus high SMAD2 blood expression was correlated with case status (13.4% lower expression; OR, 3.93; 95% CI, 2.12-7.28; P < .001). In addition, lower expression of SMAD2 was associated with more frequent exacerbations among Puerto Ricans with asthma. CONCLUSION: Ancestry at 18q21 was significantly associated with asthma in Latinos and implicated multiple ancestry-informative noncoding variants upstream of SMAD2 with asthma susceptibility. Furthermore, decreased SMAD2 expression in blood was strongly associated with increased asthma risk and increased exacerbations.

Assembly of a pan-genome from deep sequencing of 910 humans of African descent.

We used a deeply sequenced dataset of 910 individuals, all of African descent, to construct a set of DNA sequences that is present in these individuals but missing from the reference human genome. We aligned 1.19 trillion reads from the 910 individuals to the reference genome (GRCh38), collected all reads that failed to align, and assembled these reads into contiguous sequences (contigs). We then compared all contigs to one another to identify a set of unique sequences representing regions of the African pan-genome missing from the reference genome. Our analysis revealed 296,485,284 bp in 125,715 distinct contigs present in the populations of African descent, demonstrating that the African pan-genome contains ~10% more DNA than the current human reference genome. Although the functional significance of nearly all of this sequence is unknown, 387 of the novel contigs fall within 315 distinct protein-coding genes, and the rest appear to be intergenic.

Gene Coexpression Networks in Whole Blood Implicate Multiple Interrelated Molecular Pathways in Obesity in People with Asthma.

OBJECTIVE: Asthmatic children who develop obesity through adolescence have poorer disease outcomes compared with those who do not. This study aimed to characterize the biology of childhood asthma complicated by adult obesity. METHODS: Gene expression networks are powerful statistical tools for characterizing human disease that leverage the putative coregulatory relationships of genes to infer relevant biological pathways. Weighted gene coexpression network analysis of gene expression data was performed in whole blood from 514 adult asthmatic subjects. Then, module preservation and association replication analyses were performed in 418 subjects from two independent asthma cohorts (one pediatric and one adult). RESULTS: A multivariate model was identified in which three gene coexpression network modules were associated with incident obesity in the discovery cohort (each P < 0.05). Two module memberships were enriched for genes in pathways related to platelets, integrins, extracellular matrix, smooth muscle, NF-κB signaling, and Hedgehog signaling. The network structures of each of the obesity modules were significantly preserved in both replication cohorts (permutation P = 9.999E-05). The corresponding module gene sets were significantly enriched for differential expression in subjects with obesity in both replication cohorts (each P < 0.05). CONCLUSIONS: The gene coexpression network profiles thus implicate multiple interrelated pathways in the biology of an important endotype of asthma with obesity.

Parent of origin gene expression in a founder population identifies two new candidate imprinted genes at known imprinted regions.

Genomic imprinting is the phenomena that leads to silencing of one copy of a gene inherited from a specific parent. Mutations in imprinted regions have been involved in diseases showing parent of origin effects. Identifying genes with evidence of parent of origin expression patterns in family studies allows the detection of more subtle imprinting. Here, we use allele specific expression in lymphoblastoid cell lines from 306 Hutterites related in a single pedigree to provide formal evidence for parent of origin effects. We take advantage of phased genotype data to assign parent of origin to RNA-seq reads in individuals with gene expression data. Our approach identified known imprinted genes, two putative novel imprinted genes, PXDC1 and PWAR6, and 14 genes with asymmetrical parent of origin gene expression. We used gene expression in peripheral blood leukocytes (PBL) to validate our findings, and then confirmed imprinting control regions (ICRs) using DNA methylation levels in the PBLs.

Using an atlas of gene regulation across 44 human tissues to inform complex disease- and trait-associated variation.

We apply integrative approaches to expression quantitative loci (eQTLs) from 44 tissues from the Genotype-Tissue Expression project and genome-wide association study data. About 60% of known trait-associated loci are in linkage disequilibrium with a cis-eQTL, over half of which were not found in previous large-scale whole blood studies. Applying polygenic analyses to metabolic, cardiovascular, anthropometric, autoimmune, and neurodegenerative traits, we find that eQTLs are significantly enriched for trait associations in relevant pathogenic tissues and explain a substantial proportion of the heritability (40-80%). For most traits, tissue-shared eQTLs underlie a greater proportion of trait associations, although tissue-specific eQTLs have a greater contribution to some traits, such as blood pressure. By integrating information from biological pathways with eQTL target genes and applying a gene-based approach, we validate previously implicated causal genes and pathways, and propose new variant and gene associations for several complex traits, which we replicate in the UK BioBank and BioVU.

Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics.

Scalable, integrative methods to understand mechanisms that link genetic variants with phenotypes are needed. Here we derive a mathematical expression to compute PrediXcan (a gene mapping approach) results using summary data (S-PrediXcan) and show its accuracy and general robustness to misspecified reference sets. We apply this framework to 44 GTEx tissues and 100+ phenotypes from GWAS and meta-analysis studies, creating a growing public catalog of associations that seeks to capture the effects of gene expression variation on human phenotypes. Replication in an independent cohort is shown. Most of the associations are tissue specific, suggesting context specificity of the trait etiology. Colocalized significant associations in unexpected tissues underscore the need for an agnostic scanning of multiple contexts to improve our ability to detect causal regulatory mechanisms. Monogenic disease genes are enriched among significant associations for related traits, suggesting that smaller alterations of these genes may cause a spectrum of milder phenotypes.