Building a vertically integrated genomic learning health system: The biobank at the Colorado Center for Personalized Medicine.

Precision medicine initiatives across the globe have led to a revolution of repositories linking large-scale genomic data with electronic health records, enabling genomic analyses across the entire phenome. Many of these initiatives focus solely on research insights, leading to limited direct benefit to patients. We describe the biobank at the Colorado Center for Personalized Medicine (CCPM Biobank) that was jointly developed by the University of Colorado Anschutz Medical Campus and UCHealth to serve as a unique, dual-purpose research and clinical resource accelerating personalized medicine. This living resource currently has more than 200,000 participants with ongoing recruitment. We highlight the clinical, laboratory, regulatory, and HIPAA-compliant informatics infrastructure along with our stakeholder engagement, consent, recontact, and participant engagement strategies. We characterize aspects of genetic and geographic diversity unique to the Rocky Mountain region, the primary catchment area for CCPM Biobank participants. We leverage linked health and demographic information of the CCPM Biobank participant population to demonstrate the utility of the CCPM Biobank to replicate complex trait associations in the first 33,674 genotyped individuals across multiple disease domains. Finally, we describe our current efforts toward return of clinical genetic test results, including high-impact pathogenic variants and pharmacogenetic information, and our broader goals as the CCPM Biobank continues to grow. Bringing clinical and research interests together fosters unique clinical and translational questions that can be addressed from the large EHR-linked CCPM Biobank resource within a HIPAA- and CLIA-certified environment.

Xanthine oxidoreductase gene polymorphisms are associated with high risk of sepsis and organ failure.

BACKGROUND: Sepsis and associated organ failures confer substantial morbidity and mortality. Xanthine oxidoreductase (XOR) is implicated in the development of tissue oxidative damage in a wide variety of respiratory and cardiovascular disorders including sepsis and sepsis-associated acute respiratory distress syndrome (ARDS). We examined whether single nucleotide polymorphisms (SNPs) in the XDH gene (encoding XOR) might influence susceptibility to and outcome in patients with sepsis. METHODS: We genotyped 28 tag SNPs in XDH gene in the CELEG cohort, including 621 European American (EA) and 353 African American (AA) sepsis patients. Serum XOR activity was measured in a subset of CELEG subjects. Additionally, we assessed the functional effects of XDH variants utilizing empirical data from different integrated software tools and datasets. RESULTS: Among AA patients, six intronic variants (rs206805, rs513311, rs185925, rs561525, rs2163059, rs13387204), in a region enriched with regulatory elements, were associated with risk of sepsis (P < 0.008-0.049). Two out of six SNPs (rs561525 and rs2163059) were associated with risk of sepsis-associated ARDS in an independent validation cohort (GEN-SEP) of 590 sepsis patients of European descent. Two common SNPs (rs1884725 and rs4952085) in tight linkage disequilibrium (LD) provided strong evidence for association with increased levels of serum creatinine (Padjusted<0.0005 and 0.0006, respectively), suggesting a role in increased risk of renal dysfunction. In contrast, among EA ARDS patients, the missense variant rs17011368 (I703V) was associated with enhanced mortality at 60-days (P < 0.038). We found higher serum XOR activity in 143 sepsis patients (54.5 ± 57.1 mU/mL) compared to 31 controls (20.9 ± 12.4 mU/mL, P = 1.96 × 10- 13). XOR activity was associated with the lead variant rs185925 among AA sepsis patients with ARDS (P < 0.005 and Padjusted<0.01). Multifaceted functions of prioritized XDH variants, as suggested by various functional annotation tools, support their potential causality in sepsis. CONCLUSIONS: Our findings suggest that XOR is a novel combined genetic and biochemical marker for risk and outcome in patients with sepsis and ARDS.

Polygenic prediction of atopic dermatitis improves with atopic training and filaggrin factors.

BACKGROUND: While numerous genetic loci associated with atopic dermatitis (AD) have been discovered, to date, work leveraging the combined burden of AD risk variants across the genome to predict disease risk has been limited. OBJECTIVES: This study aims to determine whether polygenic risk scores (PRSs) relying on genetic determinants for AD provide useful predictions for disease occurrence and severity. It also explicitly tests the value of including genome-wide association studies of related allergic phenotypes and known FLG loss-of-function (LOF) variants. METHODS: AD PRSs were constructed for 1619 European American individuals from the Atopic Dermatitis Research Network using an AD training dataset and an atopic training dataset including AD, childhood onset asthma, and general allergy. Additionally, whole genome sequencing data were used to explore genetic scoring specific to FLG LOF mutations. RESULTS: Genetic scores derived from the AD-only genome-wide association studies were predictive of AD cases (PRSAD: odds ratio [OR], 1.70; 95% CI, 1.49-1.93). Accuracy was first improved when PRSs were built off the larger atopy genome-wide association studies (PRSAD+: OR, 2.16; 95% CI, 1.89-2.47) and further improved when including FLG LOF mutations (PRSAD++: OR, 3.23; 95% CI, 2.57-4.07). Importantly, while all 3 PRSs correlated with AD severity, the best prediction was from PRSAD++, which distinguished individuals with severe AD from control subjects with OR of 3.86 (95% CI, 2.77-5.36). CONCLUSIONS: This study demonstrates how PRSs for AD that include genetic determinants across atopic phenotypes and FLG LOF variants may be a promising tool for identifying individuals at high risk for developing disease and specifically severe disease.

Genome-wide association study of asthma, total IgE, and lung function in a cohort of Peruvian children.

BACKGROUND: Genetic ancestry plays a role in asthma health disparities. OBJECTIVE: Our aim was to evaluate the impact of ancestry on and identify genetic variants associated with asthma, total serum IgE level, and lung function. METHODS: A total of 436 Peruvian children (aged 9-19 years) with asthma and 291 without asthma were genotyped by using the Illumina Multi-Ethnic Global Array. Genome-wide proportions of indigenous ancestry populations from continental America (NAT) and European ancestry from the Iberian populations in Spain (IBS) were estimated by using ADMIXTURE. We assessed the relationship between ancestry and the phenotypes and performed a genome-wide association study. RESULTS: The mean ancestry proportions were 84.7% NAT (case patients, 84.2%; controls, 85.4%) and 15.3% IBS (15.8%; 14.6%). With adjustment for asthma, NAT was associated with higher total serum IgE levels (P < .001) and IBS was associated with lower total serum IgE levels (P < .001). NAT was associated with higher FEV1 percent predicted values (P < .001), whereas IBS was associated with lower FEV1 values in the controls but not in the case patients. The HLA-DR/DQ region on chromosome 6 (Chr6) was strongly associated with total serum IgE (rs3135348; P = 3.438 × 10-10) and was independent of an association with the haplotype HLA-DQA1∼HLA-DQB1:04.01∼04.02 (P = 1.55 × 10-05). For lung function, we identified a locus (rs4410198; P = 5.536 × 10-11) mapping to Chr19, near a cluster of zinc finger interacting genes that colocalizes to the long noncoding RNA CTD-2537I9.5. This novel locus was replicated in an independent sample of pediatric case patients with asthma with similar admixture from Brazil (P = .005). CONCLUSION: This study confirms the role of HLA in atopy, and identifies a novel locus mapping to a long noncoding RNA for lung function that may be specific to children with NAT.

Multiethnic genome-wide and HLA association study of total serum IgE level.

BACKGROUND: Total serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE. OBJECTIVE: We aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901). METHODS: We performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data. RESULTS: We identified 6 loci at genome-wide significance (P < 5 × 10-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (Pdiscovery = 2 × 10-4; Preplication = 5 × 10-4; Pdiscovery+replication = 4 × 10-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (PHispanic/Latino discovery+replication = 8 × 10-8). CONCLUSION: We performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.

Whole genome sequencing identifies novel genetic mutations in patients with eczema herpeticum.

BACKGROUND: Eczema herpeticum (EH) is a rare complication of atopic dermatitis (AD) caused by disseminated herpes simplex virus (HSV) infection. The role of rare and/or deleterious genetic variants in disease etiology is largely unknown. This study aimed to identify genes that harbor damaging genetic variants associated with HSV infection in AD with a history of recurrent eczema herpeticum (ADEH+). METHODS: Whole genome sequencing (WGS) was performed on 49 recurrent ADEH+ (≥3 EH episodes), 491 AD without a history of eczema herpeticum (ADEH-) and 237 non-atopic control (NA) subjects. Variants were annotated, and a gene-based approach (SKAT-O) was used to identify genes harboring damaging genetic variants associated with ADEH+. Genes identified through WGS were studied for effects on HSV responses and keratinocyte differentiation. RESULTS: Eight genes were identified in the comparison of recurrent ADEH+to ADEH-and NA subjects: SIDT2, CLEC7A, GSTZ1, TPSG1, SP110, RBBP8NL, TRIM15, and FRMD3. Silencing SIDT2 and RBBP8NL in normal human primary keratinocytes (NHPKs) led to significantly increased HSV-1 replication. SIDT2-silenced NHPKs had decreased gene expression of IFNk and IL1b in response to HSV-1 infection. RBBP8NL-silenced NHPKs had decreased gene expression of IFNk, but increased IL1b. Additionally, silencing SIDT2 and RBBP8NL also inhibited gene expression of keratinocyte differentiation markers keratin 10 (KRT10) and loricrin (LOR). CONCLUSION: SIDT2 and RBBP8NL participate in keratinocyte's response to HSV-1 infection. SIDT2 and RBBP8NL also regulate expression of keratinocyte differentiation genes of KRT10 and LOR.

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.

The pharmacogenomics of inhaled corticosteroids and lung function decline in COPD.

Inhaled corticosteroids (ICS) are widely prescribed for patients with chronic obstructive pulmonary disease (COPD), yet have variable outcomes and adverse reactions, which may be genetically determined. The primary aim of the study was to identify the genetic determinants for forced expiratory volume in 1 s (FEV1) changes related to ICS therapy.In the Lung Health Study (LHS)-2, 1116 COPD patients were randomised to the ICS triamcinolone acetonide (n=559) or placebo (n=557) with spirometry performed every 6 months for 3 years. We performed a pharmacogenomic genome-wide association study for the genotype-by-ICS treatment effect on 3 years of FEV1 changes (estimated as slope) in 802 genotyped LHS-2 participants. Replication was performed in 199 COPD patients randomised to the ICS, fluticasone or placebo.A total of five loci showed genotype-by-ICS interaction at p<5×10-6; of these, single nucleotide polymorphism (SNP) rs111720447 on chromosome 7 was replicated (discovery p=4.8×10-6, replication p=5.9×10-5) with the same direction of interaction effect. ENCODE (Encyclopedia of DNA Elements) data revealed that in glucocorticoid-treated (dexamethasone) A549 alveolar cell line, glucocorticoid receptor binding sites were located near SNP rs111720447. In stratified analyses of LHS-2, genotype at SNP rs111720447 was significantly associated with rate of FEV1 decline in patients taking ICS (C allele ß 56.36 mL·year-1, 95% CI 29.96-82.76 mL·year-1) and in patients who were assigned to placebo, although the relationship was weaker and in the opposite direction to that in the ICS group (C allele ß -27.57 mL·year-1, 95% CI -53.27- -1.87 mL·year-1).The study uncovered genetic factors associated with FEV1 changes related to ICS in COPD patients, which may provide new insight on the potential biology of steroid responsiveness in COPD.

The genetics of smoking in individuals with chronic obstructive pulmonary disease.

BACKGROUND: Smoking is the principal modifiable environmental risk factor for chronic obstructive pulmonary disease (COPD) which affects 300 million people and is the 3rd leading cause of death worldwide. Most of the genetic studies of smoking have relied on self-reported smoking status which is vulnerable to reporting and recall bias. Using data from the Lung Health Study (LHS), we sought to identify genetic variants associated with quantitative smoking and cessation in individuals with mild to moderate COPD. METHODS: The LHS is a longitudinal multicenter study of mild-to-moderate COPD subjects who were all smokers at recruitment. We performed genome-wide association studies (GWASs) for salivary cotinine (n = 4024), exhaled carbon monoxide (eCO) (n = 2854), cigarettes per day (CPD) (n = 2706) and smoking cessation at year 5 follow-up (n = 717 quitters and 2175 smokers). The GWAS analyses were adjusted for age, gender, and genetic principal components. RESULTS: For cotinine levels, SNPs near UGT2B10 gene achieved genome-wide significance (i.e. P < 5 × 10- 8) with top SNP rs10023464, P = 1.27 × 10- 11. For eCO levels, one significant SNP was identified which mapped to the CHRNA3 gene (rs12914385, P = 2.38 × 10- 8). A borderline region mapping to KCNMA1 gene was associated with smoking cessation (rs207675, P = 5.95 × 10- 8). Of the identified loci, only the CHRNA3/5 locus showed significant associations with lung function but only in heavy smokers. No regions met genome-wide significance for CPD. CONCLUSION: The study demonstrates that using objective measures of smoking such as eCO and/or salivary cotinine can more precisely capture the genetic contribution to multiple aspects of smoking behaviour. The KCNMA1 gene association with smoking cessation may represent a potential therapeutic target and warrants further studies. TRIAL REGISTRATION: The Lung Health Study ClinicalTrials.gov Identifier: NCT00000568 . Date of registration: October 28, 1999.

The role of ST2 and ST2 genetic variants in schistosomiasis.

BACKGROUND: Chronic schistosomiasis and its severe complication, periportal fibrosis, are characterized by a predominant Th2 response. To date, specific single nucleotide polymorphisms in ST2 have been some of the most consistently associated genetic variants for asthma. OBJECTIVE: We investigated the role of ST2 (a receptor for the Th2 cytokine IL-33) in chronic and late-stage schistosomiasis caused by Schistosoma japonicum and the potential effect of ST2 genetic variants on stage of disease and ST2 expression. METHODS: We recruited 947 adult participants (339 with end-stage schistosomiasis and liver cirrhosis, 307 with chronic infections without liver fibrosis, and 301 health controls) from a S japonicum-endemic area (Hubei, China). Six ST2 single nucleotide polymorphisms were genotyped. Serum soluble ST2 (sST2) was measured by ELISA, and ST2 expression in normal liver tissues, Hepatitis B virus-induced fibrotic liver tissues, and S japonicum-induced fibrotic liver tissues was measured by immunohistochemistry. RESULTS: We found sST2 levels were significantly higher in the end-stage group (36.04 [95% CI, 33.85-38.37]) compared with chronic cases and controls (22.7 [95% CI, 22.0-23.4], P < 1E-10). In addition, S japonicum-induced fibrotic liver tissues showed increased ST2 staining compared with normal liver tissues (P = .0001). Markers rs12712135, rs1420101, and rs6543119 were strongly associated with sST2 levels (P = 2E-10, 5E-05, and 6E-05, respectively), and these results were replicated in an independent cohort from Brazil living in a S mansoni endemic region. CONCLUSIONS: We demonstrate for the first time that end-stage schistosomiasis is associated with elevated sST2 levels and show that ST2 genetic variants are associated with sST2 levels in patients with schistosomiasis.

Genetic variants affecting cross-sectional lung function in adults show little or no effect on longitudinal lung function decline.

BACKGROUND: Genome-wide association studies have identified numerous genetic regions that influence cross-sectional lung function. Longitudinal decline in lung function also includes a heritable component but the genetic determinants have yet to be defined. OBJECTIVES: We aimed to determine whether regions associated with cross-sectional lung function were also associated with longitudinal decline and to seek novel variants which influence decline. METHODS: We analysed genome-wide data from 4167 individuals from the Busselton Health Study cohort, who had undergone spirometry (12 695 observations across eight time points). A mixed model was fitted and weighted risk scores were calculated for the joint effect of 26 known regions on baseline and longitudinal changes in FEV1 and FEV1/FVC. Potential additional regions of interest were identified and followed up in two independent cohorts. RESULTS: The 26 regions previously associated with cross-sectional lung function jointly showed a strong effect on baseline lung function (p=4.44×10-16 for FEV1/FVC) but no effect on longitudinal decline (p=0.160 for FEV1/FVC). This was replicated in an independent cohort. 39 additional regions of interest (48 variants) were identified; these associations were not replicated in two further cohorts. CONCLUSIONS: Previously identified genetic variants jointly have a strong effect on cross-sectional lung function in adults but little or no effect on the rate of decline of lung function. It is possible that they influence COPD risk through lung development. Although no genetic variants have yet been associated with lung function decline at stringent genome-wide significance, longitudinal change in lung function is heritable suggesting that there is scope for future discoveries.

Do COPD subtypes really exist? COPD heterogeneity and clustering in 10 independent cohorts.

BACKGROUND: COPD is a heterogeneous disease, but there is little consensus on specific definitions for COPD subtypes. Unsupervised clustering offers the promise of 'unbiased' data-driven assessment of COPD heterogeneity. Multiple groups have identified COPD subtypes using cluster analysis, but there has been no systematic assessment of the reproducibility of these subtypes. OBJECTIVE: We performed clustering analyses across 10 cohorts in North America and Europe in order to assess the reproducibility of (1) correlation patterns of key COPD-related clinical characteristics and (2) clustering results. METHODS: We studied 17 146 individuals with COPD using identical methods and common COPD-related characteristics across cohorts (FEV1, FEV1/FVC, FVC, body mass index, Modified Medical Research Council score, asthma and cardiovascular comorbid disease). Correlation patterns between these clinical characteristics were assessed by principal components analysis (PCA). Cluster analysis was performed using k-medoids and hierarchical clustering, and concordance of clustering solutions was quantified with normalised mutual information (NMI), a metric that ranges from 0 to 1 with higher values indicating greater concordance. RESULTS: The reproducibility of COPD clustering subtypes across studies was modest (median NMI range 0.17-0.43). For methods that excluded individuals that did not clearly belong to any cluster, agreement was better but still suboptimal (median NMI range 0.32-0.60). Continuous representations of COPD clinical characteristics derived from PCA were much more consistent across studies. CONCLUSIONS: Identical clustering analyses across multiple COPD cohorts showed modest reproducibility. COPD heterogeneity is better characterised by continuous disease traits coexisting in varying degrees within the same individual, rather than by mutually exclusive COPD subtypes.

Surfactant protein D is a causal risk factor for COPD: results of Mendelian randomisation.

Surfactant protein D (SP-D) is produced primarily in the lung and is involved in regulating pulmonary surfactants, lipid homeostasis and innate immunity. Circulating SP-D levels in blood are associated with chronic obstructive pulmonary disease (COPD), although causality remains elusive.In 4061 subjects with COPD, we identified genetic variants associated with serum SP-D levels. We then determined whether these variants affected lung tissue gene expression in 1037 individuals. A Mendelian randomisation framework was then applied, whereby serum SP-D-associated variants were tested for association with COPD risk in 11 157 cases and 36 699 controls and with 11 years decline of lung function in the 4061 individuals.Three regions on chromosomes 6 (human leukocyte antigen region), 10 (SFTPD gene) and 16 (ATP2C2 gene) were associated with serum SP-D levels at genome-wide significance. In Mendelian randomisation analyses, variants associated with increased serum SP-D levels decreased the risk of COPD (estimate -0.19, p=6.46×10-03) and slowed the lung function decline (estimate=0.0038, p=7.68×10-3).Leveraging genetic variation effect on protein, lung gene expression and disease phenotypes provided novel insights into SP-D biology and established a causal link between increased SP-D levels and protection against COPD risk and progression. SP-D represents a very promising biomarker and therapeutic target for COPD.

Effect of polymorphisms on TGFB1 on allergic asthma and helminth infection in an African admixed population.

BACKGROUND: Allergic asthma is a complex disorder that results from a combination of genetic and environmental factors. Studies suggest that helminth infections can activate a regulatory network characterized by the production of regulatory cytokines, such as interleukin 10 and transforming growth factor ß1 (TGF-ß1) and subsequently protect against immune-mediated diseases, such as asthma. On the other hand, TGF-ß1 is increased in the lungs of individuals with asthma and may modulate airway inflammation. The role of TGF- ß 1 single-nucleotide polymorphisms (SNPs) in allergic disease remains inconclusive. OBJECTIVE: To evaluate the effects of genetic variations in the TGF-ß1 on allergy and helminths infections in children. METHODS: We tested for association among 4 TGF-ß1 SNPs and allergic asthma, specific IgE, skin prick test result, and IL-10 production in 1,335 Brazilians. In addition, we analyzed the association with markers of helminth infection (parasite burden, anti-Ascaris IgE, and worm specific IgG4). The polymorphisms were genotyped using Taq Man probes. RESULTS: We found an association between rs1800470 (C allele) and atopic wheezing (odds ratio [OR], 0.60; 95% confidence interval [CI], 0.37-0.95) and markers of allergy (OR, 0.41; 95% CI, 0.22-0.79). In contrast, a positive association was observed between the haplotype ACCA and Trichuris trichiura infection (OR, 1.85; P = .003) and Ascaris lumbricoides infection (OR, 2.01; P < .001). This haplotype was also associated with increased IL-10 production (ß = 50.7; P < .001). CONCLUSION: Individuals with TGF-ß1 polymorphisms have an increased susceptibility to helminth infections and a lower risk of developing allergy. These studies suggest that immune modulation of allergic disease results not only from environmental factors but also from genetic susceptibility and IL-10 production.

Serum endostatin is a genetically determined predictor of survival in pulmonary arterial hypertension.

RATIONALE: Pulmonary arterial hypertension (PAH) is a medically incurable disease resulting in death from right ventricular (RV) failure. Both pulmonary vascular and RV remodeling are linked to dynamic changes in the microvasculature. Therefore, we hypothesized that circulating angiostatic factors could be linked to outcomes and represent novel biomarkers of disease severity in PAH. OBJECTIVES: We sought to determine the relationship of a potent angiostatic factor, endostatin (ES), with disease severity and mortality in PAH. Furthermore, we assessed genetic predictors of ES expression and/or function and their association with outcomes in PAH. METHODS: We measured levels of serum ES in two independent cohorts of patients with PAH. Contemporaneous clinical data included New York Heart Association functional class, 6-minute-walk distance, invasive hemodynamics, and laboratory chemistries. MEASUREMENTS AND MAIN RESULTS: Serum ES correlated with poor functional status, decreased exercise tolerance, and invasive hemodynamics variables. Furthermore, serum ES was a strong predictor of mortality. A loss-of-function, missense variant in the gene encoding ES, Col18a1, was linked to lower circulating protein and was independently associated with reduced mortality. CONCLUSIONS: Our data link increased expression of ES to disease severity in PAH and demonstrate a significant relationship with adverse outcomes. Circulating ES levels can be genetically influenced, implicating ES as a genetically determined modifier of disease severity impacting on survival. These observations support serum ES as a potential biomarker in PAH with the capacity to predict poor outcomes. More importantly, this study implicates Col18a1/ES as a potential new therapeutic target in PAH.

Targeted deep sequencing identifies rare loss-of-function variants in IFNGR1 for risk of atopic dermatitis complicated by eczema herpeticum.

BACKGROUND: A subset of atopic dermatitis is associated with increased susceptibility to eczema herpeticum (ADEH+). We previously reported that common single nucleotide polymorphisms (SNPs) in the IFN-γ (IFNG) and IFN-γ receptor 1 (IFNGR1) genes were associated with the ADEH+ phenotype. OBJECTIVE: We sought to interrogate the role of rare variants in interferon pathway genes for the risk of ADEH+. METHODS: We performed targeted sequencing of interferon pathway genes (IFNG, IFNGR1, IFNAR1, and IL12RB1) in 228 European American patients with AD selected according to their eczema herpeticum status, and severity was measured by using the Eczema Area and Severity Index. Replication genotyping was performed in independent samples of 219 European American and 333 African American subjects. Functional investigation of loss-of-function variants was conducted by using site-directed mutagenesis. RESULTS: We identified 494 single nucleotide variants encompassing 105 kb of sequence, including 145 common, 349 (70.6%) rare (minor allele frequency <5%), and 86 (17.4%) novel variants, of which 2.8% were coding synonymous, 93.3% were noncoding (64.6% intronic), and 3.8% were missense. We identified 6 rare IFNGR1 missense variants, including 3 damaging variants (Val14Met [V14M], Val61Ile, and Tyr397Cys [Y397C]) conferring a higher risk for ADEH+ (P = .031). Variants V14M and Y397C were confirmed to be deleterious, leading to partial IFNGR1 deficiency. Seven common IFNGR1 SNPs, along with common protective haplotypes (2-7 SNPs), conferred a reduced risk of ADEH+ (P = .015-.002 and P = .0015-.0004, respectively), and both SNP and haplotype associations were replicated in an independent African American sample (P = .004-.0001 and P = .001-.0001, respectively). CONCLUSION: Our results provide evidence that both genetic variants in the gene encoding IFNGR1 are implicated in susceptibility to the ADEH+ phenotype.

Genome-Wide Association Study Identification of Novel Loci Associated with Airway Responsiveness in Chronic Obstructive Pulmonary Disease.

Increased airway responsiveness is linked to lung function decline and mortality in subjects with chronic obstructive pulmonary disease (COPD); however, the genetic contribution to airway responsiveness remains largely unknown. A genome-wide association study (GWAS) was performed using the Illumina (San Diego, CA) Human660W-Quad BeadChip on European Americans with COPD from the Lung Health Study. Linear regression models with correlated meta-analyses, including data from baseline (n = 2,814) and Year 5 (n = 2,657), were used to test for common genetic variants associated with airway responsiveness. Genotypic imputation was performed using reference 1000 Genomes Project data. Expression quantitative trait loci (eQTL) analyses in lung tissues were assessed for the top 10 markers identified, and immunohistochemistry assays assessed protein staining for SGCD and MYH15. Four genes were identified within the top 10 associations with airway responsiveness. Markers on chromosome 9p21.2 flanked by LINGO2 met a predetermined threshold of genome-wide significance (P < 9.57 × 10(-8)). Markers on chromosomes 3q13.1 (flanked by MYH15), 5q33 (SGCD), and 6q21 (PDSS2) yielded suggestive evidence of association (9.57 × 10(-8) < P ≤ 4.6 × 10(-6)). Gene expression studies in lung tissue showed single nucleotide polymorphisms on chromosomes 5 and 3 to act as eQTL for SGCD (P = 2.57 × 10(-9)) and MYH15 (P = 1.62 × 10(-6)), respectively. Immunohistochemistry confirmed localization of SGCD protein to airway smooth muscle and vessels and MYH15 to airway epithelium, vascular endothelium, and inflammatory cells. We identified novel loci associated with airway responsiveness in a GWAS among smokers with COPD. Risk alleles on chromosomes 5 and 3 acted as eQTLs for SGCD and MYH15 messenger RNA, and these proteins were expressed in lung cells relevant to the development of airway responsiveness.

APOε4 is associated with enhanced in vivo innate immune responses in human subjects.

BACKGROUND: The genetic determinants of the human innate immune response are poorly understood. Apolipoprotein (Apo) E, a lipid-trafficking protein that affects inflammation, has well-described wild-type (ε3) and disease-associated (ε2 and ε4) alleles, but its connection to human innate immunity is undefined. OBJECTIVE: We sought to define the relationship of APOε4 to the human innate immune response. METHODS: We evaluated APOε4 in several functional models of the human innate immune response, including intravenous LPS challenge in human subjects, and assessed APOε4 association to organ injury in patients with severe sepsis, a disease driven by dysregulated innate immunity. RESULTS: Whole blood from healthy APOε3/APOε4 volunteers induced higher cytokine levels on ex vivo stimulation with Toll-like receptor (TLR) 2, TLR4, or TLR5 ligands than blood from APOε3/APOε3 patients, whereas TLR7/8 responses were similar. This was associated with increased lipid rafts in APOε3/APOε4 monocytes. By contrast, APOε3/APOε3 and APOε3/APOε4 serum neutralized LPS equivalently and supported similar LPS responses in Apoe-deficient macrophages, arguing against a differential role for secretory APOE4 protein. After intravenous LPS, APOε3/APOε4 patients had higher hyperthermia and plasma TNF-α levels and earlier plasma IL-6 than APOε3/APOε3 patients. APOE4-targeted replacement mice displayed enhanced hypothermia, plasma cytokines, and hepatic injury and altered splenic lymphocyte apoptosis after systemic LPS compared with APOE3 counterparts. In a cohort of 828 patients with severe sepsis, APOε4 was associated with increased coagulation system failure among European American patients. CONCLUSIONS: APOε4 is a determinant of the human innate immune response to multiple TLR ligands and associates with altered patterns of organ injury in human sepsis.

Genome-wide association study of lung function phenotypes in a founder population.

BACKGROUND: Lung function is a long-term predictor of mortality and morbidity. OBJECTIVE: We sought to identify single nucleotide polymorphisms (SNPs) associated with lung function. METHODS: We performed a genome-wide association study (GWAS) of FEV1, forced vital capacity (FVC), and FEV1/FVC in 1144 Hutterites aged 6 to 89 years, who are members of a founder population of European descent. We performed least absolute shrinkage and selection operation regression to select the minimum set of SNPs that best predict FEV1/FVC in the Hutterites and used the GRAIL algorithm to mine the Gene Ontology database for evidence of functional connections between genes near the predictive SNPs. RESULTS: Our GWAS identified significant associations between FEV1/FVC and SNPs at the THSD4-UACA-TLE3 locus on chromosome 15q23 (P = 5.7 × 10(-8) to 3.4 × 10(-9)). Nine SNPs at or near 4 additional loci had P < 10(-5) with FEV1/FVC. Only 2 SNPs were found with P < 10(-5) for FEV1 or FVC. We found nominal levels of significance with SNPs at 9 of the 27 previously reported loci associated with lung function measures. Among a predictive set of 80 SNPs, 6 loci were identified that had a significant degree of functional connectivity (GRAIL P < .05), including 3 clusters of ß-defensin genes, 2 chemokine genes (CCL18 and CXCL12), and TNFRSF13B. CONCLUSION: This study identifies genome-wide significant associations and replicates results of previous GWASs. Multimarker modeling implicated for the first time common variation in genes involved in antimicrobial immunity in airway mucosa that influences lung function.

African ancestry is a risk factor for asthma and high total IgE levels in African admixed populations.

Characterization of genetic admixture of populations in the Americas and the Caribbean is of interest for anthropological, epidemiological, and historical reasons. Asthma has a higher prevalence and is more severe in populations with a high African component. Association of African ancestry with asthma has been demonstrated. We estimated admixture proportions of samples from six trihybrid populations of African descent and determined the relationship between African ancestry and asthma and total serum IgE levels (tIgE). We genotyped 237 ancestry informative markers in asthmatics and nonasthmatic controls from Barbados (190/277), Jamaica (177/529), Brazil (40/220), Colombia (508/625), African Americans from New York (207/171), and African Americans from Baltimore/Washington, D.C. (625/757). We estimated individual ancestries and evaluated genetic stratification using Structure and principal component analysis. Association of African ancestry and asthma and tIgE was evaluated by regression analysis. Mean ± SD African ancestry ranged from 0.76 ± 0.10 among Barbadians to 0.33 ± 0.13 in Colombians. The European component varied from 0.14 ± 0.05 among Jamaicans and Barbadians to 0.26 ± 0.08 among Colombians. African ancestry was associated with risk for asthma in Colombians (odds ratio (OR) = 4.5, P = 0.001) Brazilians (OR = 136.5, P = 0.003), and African Americans of New York (OR: 4.7; P = 0.040). African ancestry was also associated with higher tIgE levels among Colombians (ß = 1.3, P = 0.04), Barbadians (ß = 3.8, P = 0.03), and Brazilians (ß = 1.6, P = 0.03). Our findings indicate that African ancestry can account for, at least in part, the association between asthma and its associated trait, tIgE levels.