Prediction of disease-free survival for precision medicine using cooperative learning on multi-omic data.

In precision medicine, both predicting the disease susceptibility of an individual and forecasting its disease-free survival are areas of key research. Besides the classical epidemiological predictor variables, data from multiple (omic) platforms are increasingly available. To integrate this wealth of information, we propose new methodology to combine both cooperative learning, a recent approach to leverage the predictive power of several datasets, and polygenic hazard score models. Polygenic hazard score models provide a practitioner with a more differentiated view of the predicted disease-free survival than the one given by merely a point estimate, for instance computed with a polygenic risk score. Our aim is to leverage the advantages of cooperative learning for the computation of polygenic hazard score models via Cox's proportional hazard model, thereby improving the prediction of the disease-free survival. In our experimental study, we apply our methodology to forecast the disease-free survival for Alzheimer's disease (AD) using three layers of data. One layer contains epidemiological variables such as sex, APOE (apolipoprotein E, a genetic risk factor for AD) status and 10 leading principal components. Another layer contains selected genomic loci, and the last layer contains methylation data for selected CpG sites. We demonstrate that the survival curves computed via cooperative learning yield an AUC of around $0.7$, above the state-of-the-art performance of its competitors. Importantly, the proposed methodology returns (1) a linear score that can be easily interpreted (in contrast to machine learning approaches), and (2) a weighting of the predictive power of the involved data layers, allowing for an assessment of the importance of each omic (or other) platform. Similarly to polygenic hazard score models, our methodology also allows one to compute individual survival curves for each patient.

FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies.

BACKGROUND: Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations. METHODS: We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants. RESULTS: Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22. CONCLUSION: Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

A comparison between similarity matrices for principal component analysis to assess population stratification in sequenced genetic data sets.

Genetic similarity matrices are commonly used to assess population substructure (PS) in genetic studies. Through simulation studies and by the application to whole-genome sequencing (WGS) data, we evaluate the performance of three genetic similarity matrices: the unweighted and weighted Jaccard similarity matrices and the genetic relationship matrix. We describe different scenarios that can create numerical pitfalls and lead to incorrect conclusions in some instances. We consider scenarios in which PS is assessed based on loci that are located across the genome ('globally') and based on loci from a specific genomic region ('locally'). We also compare scenarios in which PS is evaluated based on loci from different minor allele frequency bins: common (>5%), low-frequency (5-0.5%) and rare (<0.5%) single-nucleotide variations (SNVs). Overall, we observe that all approaches provide the best clustering performance when computed based on rare SNVs. The performance of the similarity matrices is very similar for common and low-frequency variants, but for rare variants, the unweighted Jaccard matrix provides preferable clustering features. Based on visual inspection and in terms of standard clustering metrics, its clusters are the densest and the best separated in the principal component analysis of variants with rare SNVs compared with the other methods and different allele frequency cutoffs. In an application, we assessed the role of rare variants on local and global PS, using WGS data from multiethnic Alzheimer's disease data sets and European or East Asian populations from the 1000 Genome Project.

Blood-based Transcriptomic and Proteomic Biomarkers of Emphysema.

Rationale: Emphysema is a chronic obstructive pulmonary disease phenotype with important prognostic implications. Identifying blood-based biomarkers of emphysema will facilitate early diagnosis and development of targeted therapies. Objectives: To discover blood omics biomarkers for chest computed tomography-quantified emphysema and develop predictive biomarker panels. Methods: Emphysema blood biomarker discovery was performed using differential gene expression, alternative splicing, and protein association analyses in a training sample of 2,370 COPDGene participants with available blood RNA sequencing, plasma proteomics, and clinical data. Internal validation was conducted in a COPDGene testing sample (n = 1,016), and external validation was done in the ECLIPSE study (n = 526). Because low body mass index (BMI) and emphysema often co-occur, we performed a mediation analysis to quantify the effect of BMI on gene and protein associations with emphysema. Elastic net models with bootstrapping were also developed in the training sample sequentially using clinical, blood cell proportions, RNA-sequencing, and proteomic biomarkers to predict quantitative emphysema. Model accuracy was assessed by the area under the receiver operating characteristic curves for subjects stratified into tertiles of emphysema severity. Measurements and Main Results: Totals of 3,829 genes, 942 isoforms, 260 exons, and 714 proteins were significantly associated with emphysema (false discovery rate, 5%) and yielded 11 biological pathways. Seventy-four percent of these genes and 62% of these proteins showed mediation by BMI. Our prediction models demonstrated reasonable predictive performance in both COPDGene and ECLIPSE. The highest-performing model used clinical, blood cell, and protein data (area under the receiver operating characteristic curve in COPDGene testing, 0.90; 95% confidence interval, 0.85-0.90). Conclusions: Blood transcriptome and proteome-wide analyses revealed key biological pathways of emphysema and enhanced the prediction of emphysema.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables.

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Fast computation of the eigensystem of genomic similarity matrices.

The computation of a similarity measure for genomic data is a standard tool in computational genetics. The principal components of such matrices are routinely used to correct for biases due to confounding by population stratification, for instance in linear regressions. However, the calculation of both a similarity matrix and its singular value decomposition (SVD) are computationally intensive. The contribution of this article is threefold. First, we demonstrate that the calculation of three matrices (called the covariance matrix, the weighted Jaccard matrix, and the genomic relationship matrix) can be reformulated in a unified way which allows for the application of a randomized SVD algorithm, which is faster than the traditional computation. The fast SVD algorithm we present is adapted from an existing randomized SVD algorithm and ensures that all computations are carried out in sparse matrix algebra. The algorithm only assumes that row-wise and column-wise subtraction and multiplication of a vector with a sparse matrix is available, an operation that is efficiently implemented in common sparse matrix packages. An exception is the so-called Jaccard matrix, which does not have a structure applicable for the fast SVD algorithm. Second, an approximate Jaccard matrix is introduced to which the fast SVD computation is applicable. Third, we establish guaranteed theoretical bounds on the accuracy (in [Formula: see text] norm and angle) between the principal components of the Jaccard matrix and the ones of our proposed approximation, thus putting the proposed Jaccard approximation on a solid mathematical foundation, and derive the theoretical runtime of our algorithm. We illustrate that the approximation error is low in practice and empirically verify the theoretical runtime scalings on both simulated data and data of the 1000 Genome Project.

On the effect heterogeneity of established disease susceptibility loci for Alzheimer's disease across different genetic ancestries.

INTRODUCTION: Genome-wide association studies have identified numerous disease susceptibility loci (DSLs) for Alzheimer's disease (AD). However, only a limited number of studies have investigated the dependence of the genetic effect size of established DSLs on genetic ancestry. METHODS: We utilized the whole genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) including 35,569 participants. A total of 25,459 subjects in four distinct populations (African ancestry, non-Hispanic White, admixed Hispanic, and Asian) were analyzed. RESULTS: We found that nine DSLs showed significant heterogeneity across populations. Single nucleotide polymorphism (SNP) rs2075650 in translocase of outer mitochondrial membrane 40 (TOMM40) showed the largest heterogeneity (Cochran's Q = 0.00, I2 = 90.08), followed by other SNPs in apolipoprotein C1 (APOC1) and apolipoprotein E (APOE). Two additional loci, signal-induced proliferation-associated 1 like 2 (SIPA1L2) and solute carrier 24 member 4 (SLC24A4), showed significant heterogeneity across populations. DISCUSSION: We observed substantial heterogeneity for the APOE-harboring 19q13.32 region with TOMM40/APOE/APOC1 genes. The largest risk effect was seen among African Americans, while Asians showed a surprisingly small risk effect.

Integrating Genetics, Transcriptomics, and Proteomics in Lung Tissue to Investigate Chronic Obstructive Pulmonary Disease.

The integration of transcriptomic and proteomic data from lung tissue with chronic obstructive pulmonary disease (COPD)-associated genetic variants could provide insight into the biological mechanisms of COPD. Here, we assessed associations between lung transcriptomics and proteomics with COPD in 98 subjects from the Lung Tissue Research Consortium. Low correlations between transcriptomics and proteomics were generally observed, but higher correlations were found for COPD-associated proteins. We integrated COPD risk SNPs or SNPs near COPD-associated proteins with lung transcripts and proteins to identify regulatory cis-quantitative trait loci (QTLs). Significant expression QTLs (eQTLs) and protein QTLs (pQTLs) were found regulating multiple COPD-associated biomarkers. We investigated mediated associations from significant pQTLs through transcripts to protein levels of COPD-associated proteins. We also attempted to identify colocalized effects between COPD genome-wide association studies and eQTL and pQTL signals. Evidence was found for colocalization between COPD genome-wide association study signals and a pQTL for RHOB and an eQTL for DSP. We applied weighted gene co-expression network analysis to find consensus COPD-associated network modules. Two network modules generated by consensus weighted gene co-expression network analysis were associated with COPD with a false discovery rate lower than 0.05. One network module is related to the catenin complex, and the other module is related to plasma membrane components. In summary, multiple cis-acting determinants of transcripts and proteins associated with COPD were identified. Colocalization analysis, mediation analysis, and correlation-based network analysis of multiple omics data may identify key genes and proteins that work together to influence COPD pathogenesis.

Novel genetic variants associated with inhaled corticosteroid treatment response in older adults with asthma.

INTRODUCTION: Older adults have the greatest burden of asthma and poorest outcomes. The pharmacogenetics of inhaled corticosteroid (ICS) treatment response is not well studied in older adults. METHODS: A genome-wide association study of ICS response was performed in asthmatics of European ancestry in Genetic Epidemiology Research on Adult Health and Aging (GERA) by fitting Cox proportional hazards regression models, followed by validation in the Mass General Brigham (MGB) Biobank and Rotterdam Study. ICS response was measured using two definitions in asthmatics on ICS treatment: (1) absence of oral corticosteroid (OCS) bursts using prescription records and (2) absence of asthma-related exacerbations using diagnosis codes. A fixed-effect meta-analysis was performed for each outcome. The validated single-nucleotide polymorphisms (SNPs) were functionally annotated to standard databases. RESULTS: In 5710 subjects in GERA, 676 subjects in MGB Biobank, and 465 subjects in the Rotterdam Study, four novel SNPs on chromosome six near PTCHD4 validated across all cohorts and met genome-wide significance on meta-analysis for the OCS burst outcome. In 4541 subjects in GERA and 505 subjects in MGB Biobank, 152 SNPs with p<5 × 10-5 were validated across these two cohorts for the asthma-related exacerbation outcome. The validated SNPs included methylation and expression quantitative trait loci for CPED1, CRADD and DST for the OCS burst outcome and GM2A, SNW1, CACNA1C, DPH1, and RPS10 for the asthma-related exacerbation outcome. CONCLUSIONS: Multiple novel SNPs associated with ICS response were identified in older adult asthmatics. Several SNPs annotated to genes previously associated with asthma and other airway or allergic diseases, including PTCHD4.

Caution against examining the role of reverse causality in Mendelian Randomization.

Recently, Mendelian Randomization (MR) has gained in popularity as a concept to assess the causal relationship between phenotypes in genetic association studies. An extension of standard MR methodology, the MR Steiger approach, has recently been developed to infer the causal direction between two phenotypes in prospective studies. Through simulation studies, we examined and quantified the ability of the MR Steiger approach to determine the causal direction between two phenotypes (i.e., effect direction). Through simulation studies, our results show that the MR Steiger approach may fail to correctly identify the direction of causality. This is true, especially in the presence of pleiotropy. We also applied the MR Steiger method to the COPDGene study, a case-control study of chronic obstructive pulmonary disease (COPD) in current and former smokers, to examine the role of smoking on lung function. We have created an R package on Github called reverseDirection which runs simulations for user-specified scenarios to examine when the MR Steiger approach can correctly determine the causal direction between two phenotypes in any user specified scenario. In summary, our results emphasize the importance of caution when the MR Steiger approach is used in to infer the direction of causality.

locStra: Fast analysis of regional/global stratification in whole-genome sequencing studies.

locStra is an R -package for the analysis of regional and global population stratification in whole-genome sequencing (WGS) studies, where regional stratification refers to the substructure defined by the loci in a particular region on the genome. Population substructure can be assessed based on the genetic covariance matrix, the genomic relationship matrix, and the unweighted/weighted genetic Jaccard similarity matrix. Using a sliding window approach, the regional similarity matrices are compared with the global ones, based on user-defined window sizes and metrics, for example, the correlation between regional and global eigenvectors. An algorithm for the specification of the window size is provided. As the implementation fully exploits sparse matrix algebra and is written in C++, the analysis is highly efficient. Even on single cores, for realistic study sizes (several thousand subjects, several million rare variants per subject), the runtime for the genome-wide computation of all regional similarity matrices does typically not exceed one hour, enabling an unprecedented investigation of regional stratification across the entire genome. The package is applied to three WGS studies, illustrating the varying patterns of regional substructure across the genome and its beneficial effects on association testing.

Genome-wide association analysis of COVID-19 mortality risk in SARS-CoV-2 genomes identifies mutation in the SARS-CoV-2 spike protein that colocalizes with P.1 of the Brazilian strain.

SARS-CoV-2 mortality has been extensively studied in relation to host susceptibility. How sequence variations in the SARS-CoV-2 genome affect pathogenicity is poorly understood. Starting in October 2020, using the methodology of genome-wide association studies (GWAS), we looked at the association between whole-genome sequencing (WGS) data of the virus and COVID-19 mortality as a potential method of early identification of highly pathogenic strains to target for containment. Although continuously updating our analysis, in December 2020, we analyzed 7548 single-stranded SARS-CoV-2 genomes of COVID-19 patients in the GISAID database and associated variants with mortality using a logistic regression. In total, evaluating 29,891 sequenced loci of the viral genome for association with patient/host mortality, two loci, at 12,053 and 25,088 bp, achieved genome-wide significance (p values of 4.09e-09 and 4.41e-23, respectively), though only 25,088 bp remained significant in follow-up analyses. Our association findings were exclusively driven by the samples that were submitted from Brazil (p value of 4.90e-13 for 25,088 bp). The mutation frequency of 25,088 bp in the Brazilian samples on GISAID has rapidly increased from about 0.4 in October/December 2020 to 0.77 in March 2021. Although GWAS methodology is suitable for samples in which mutation frequencies varies between geographical regions, it cannot account for mutation frequencies that change rapidly overtime, rendering a GWAS follow-up analysis of the GISAID samples that have been submitted after December 2020 as invalid. The locus at 25,088 bp is located in the P.1 strain, which later (April 2021) became one of the distinguishing loci (precisely, substitution V1176F) of the Brazilian strain as defined by the Centers for Disease Control. Specifically, the mutations at 25,088 bp occur in the S2 subunit of the SARS-CoV-2 spike protein, which plays a key role in viral entry of target host cells. Since the mutations alter amino acid coding sequences, they potentially imposing structural changes that could enhance viral infectivity and symptom severity. Our analysis suggests that GWAS methodology can provide suitable analysis tools for the real-time detection of new more transmissible and pathogenic viral strains in databases such as GISAID, though new approaches are needed to accommodate rapidly changing mutation frequencies over time, in the presence of simultaneously changing case/control ratios. Improvements of the associated metadata/patient information in terms of quality and availability will also be important to fully utilize the potential of GWAS methodology in this field.

Increased chest CT derived bone and muscle measures capture markers of improved morbidity and mortality in COPD.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a disease of accelerated aging and is associated with comorbid conditions including osteoporosis and sarcopenia. These extrapulmonary conditions are highly prevalent yet frequently underdiagnosed and overlooked by pulmonologists in COPD treatment and management. There is evidence supporting a role for bone-muscle crosstalk which may compound osteoporosis and sarcopenia risk in COPD. Chest CT is commonly utilized in COPD management, and we evaluated its utility to identify low bone mineral density (BMD) and reduced pectoralis muscle area (PMA) as surrogates for osteoporosis and sarcopenia. We then tested whether BMD and PMA were associated with morbidity and mortality in COPD. METHODS: BMD and PMA were analyzed from chest CT scans of 8468 COPDGene participants with COPD and controls (smoking and non-smoking). Multivariable regression models tested the relationship of BMD and PMA with measures of function (6-min walk distance (6MWD), handgrip strength) and disease severity (percent emphysema and lung function). Multivariable Cox proportional hazards models were used to evaluate the relationship between sex-specific quartiles of BMD and/or PMA derived from non-smoking controls with all-cause mortality. RESULTS: COPD subjects had significantly lower BMD and PMA compared with controls. Higher BMD and PMA were associated with increased physical function and less disease severity. Participants with the highest BMD and PMA quartiles had a significantly reduced mortality risk (36% and 46%) compared to the lowest quartiles. CONCLUSIONS: These findings highlight the potential for CT-derived BMD and PMA to characterize osteoporosis and sarcopenia using equipment available in the pulmonary setting.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease.

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008-2010) and SpiroMeta Consortium (multiple countries, 1947-2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic ß4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

A polygenic risk score for asthma in a large racially diverse population.

BACKGROUND: Polygenic risk scores (PRSs) will have important utility for asthma and other chronic diseases as a tool for predicting disease incidence and subphenotypes. OBJECTIVE: We utilized findings from a large multiancestry GWAS of asthma to compute a PRS for asthma with relevance for racially diverse populations. METHODS: We derived two PRSs for asthma using a standard approach (based on genome-wide significant variants) and a lasso sum regression approach (allowing all genetic variants to potentially contribute). We used data from the racially diverse Kaiser Permanente GERA cohort (68 638 non-Hispanic Whites, 5874 Hispanics, 6870 Asians and 2760 Blacks). Race was self-reported by questionnaire. RESULTS: For the standard PRS, non-Hispanic Whites showed the highest odds ratio for a standard deviation increase in PRS for asthma (OR = 1.16 (95% CI 1.14-1.18)). The standard PRS was also associated with asthma in Hispanic (OR = 1.12 (95% CI 1.05-1.19)) and Asian (OR = 1.10 (95% CI 1.04-1.17)) subjects, with a trend towards increased risk in Blacks (OR = 1.05 (95% CI 0.97-1.15)). We detected an interaction by sex, with men showing a higher risk of asthma with an increase in PRS as compared to women. The lasso sum regression-derived PRS showed stronger associations with asthma in non-Hispanic White subjects (OR = 1.20 (95% CI 1.18-1.23)), Hispanics (OR = 1.17 (95% 1.10-1.26)), Asians (OR = 1.18 (95% CI 1.10-1.27)) and Blacks (OR = 1.10 (95% CI 0.99-1.22)). CONCLUSION: Polygenic risk scores across multiple racial/ethnic groups were associated with increased asthma risk, suggesting that PRSs have potential as a tool for predicting disease development.

Pharmacogenetics of Bronchodilator Response: Future Directions.

PURPOSE OF REVIEW: Several genome-wide association studies (GWASs) of bronchodilator response (BDR) to albuterol have been published over the past decade. This review describes current knowledge gaps, including pharmacogenetic studies of albuterol response in minority populations, effect modification of pharmacogenetic associations by age, and relevance of BDR phenotype characterization to pharmacogenetic findings. New approaches, such as leveraging additional "omics" data to focus pharmacogenetic interrogation, as well as developing polygenic risk scores in asthma treatment responses, are also discussed. RECENT FINDINGS: Recent pharmacogenetic studies of albuterol response in minority populations have identified genetic polymorphisms in loci (DNAH5, NFKB1, PLCB1, ADAMTS3, COX18, and PRKG1), that are associated with BDR. Additional studies are needed to replicate these findings. Modification of the pharmacogenetic associations for SPATS2L and ASB3 polymorphisms by age has also been published. Evidence from metabolomic and epigenomic studies of BDR may point to new pharmacogenetic targets. Lastly, a polygenic risk score for response to albuterol has been developed but requires validation in additional cohorts. In order to expand our knowledge of pharmacogenetics of BDR, additional studies in minority populations are needed. Consideration of effect modification by age and leverage of other "omics" data beyond genomics may also help uncover novel pharmacogenetic loci for use in precision medicine for asthma treatment.

Proper conditional analysis in the presence of missing data: Application to large scale meta-analysis of tobacco use phenotypes.

Meta-analysis of genetic association studies increases sample size and the power for mapping complex traits. Existing methods are mostly developed for datasets without missing values, i.e. the summary association statistics are measured for all variants in contributing studies. In practice, genotype imputation is not always effective. This may be the case when targeted genotyping/sequencing assays are used or when the un-typed genetic variant is rare. Therefore, contributed summary statistics often contain missing values. Existing methods for imputing missing summary association statistics and using imputed values in meta-analysis, approximate conditional analysis, or simple strategies such as complete case analysis all have theoretical limitations. Applying these approaches can bias genetic effect estimates and lead to seriously inflated type-I or type-II errors in conditional analysis, which is a critical tool for identifying independently associated variants. To address this challenge and complement imputation methods, we developed a method to combine summary statistics across participating studies and consistently estimate joint effects, even when the contributed summary statistics contain large amounts of missing values. Based on this estimator, we proposed a score statistic called PCBS (partial correlation based score statistic) for conditional analysis of single-variant and gene-level associations. Through extensive analysis of simulated and real data, we showed that the new method produces well-calibrated type-I errors and is substantially more powerful than existing approaches. We applied the proposed approach to one of the largest meta-analyses to date for the cigarettes-per-day phenotype. Using the new method, we identified multiple novel independently associated variants at known loci for tobacco use, which were otherwise missed by alternative methods. Together, the phenotypic variance explained by these variants was 1.1%, improving that of previously reported associations by 71%. These findings illustrate the extent of locus allelic heterogeneity and can help pinpoint causal variants.

Effect of population stratification on SNP-by-environment interaction.

Proportions of false-positive rates in genome-wide association analysis are affected by population stratification, and if it is not correctly adjusted, the statistical analysis can produce the large false-negative finding. Therefore various approaches have been proposed to adjust such problems in genome-wide association studies. However, in spite of its importance, a few studies have been conducted in genome-wide single nucleotide polymorphism (SNP)-by-environment interaction studies. In this report, we illustrate in which scenarios can lead to the false-positive rates in association mapping and approach to maintaining the overall type-1 error rate.

Assessing pleiotropy and mediation in genetic loci associated with chronic obstructive pulmonary disease.

Genetic association studies have increasingly recognized variant effects on multiple phenotypes. Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with environmental and genetic causes. Multiple genetic variants have been associated with COPD, many of which show significant associations to additional phenotypes. However, it is unknown if these associations represent biological pleiotropy or if they exist through correlation of related phenotypes ("mediated pleiotropy"). Using 6,670 subjects from the COPDGene study, we describe the association of known COPD susceptibility loci with other COPD-related phenotypes and distinguish if these act directly on the phenotypes (i.e., biological pleiotropy) or if the association is due to correlation (i.e., mediated pleiotropy). We identified additional associated phenotypes for 13 of 25 known COPD loci. Tests for pleiotropy between genotype and associated outcomes were significant for all loci. In cases of significant pleiotropy, we performed mediation analysis to test if SNPs had a direct association to phenotype. Most loci showed a mediated effect through the hypothesized causal pathway. However, many loci also had direct associations, suggesting causal explanations (i.e., emphysema leading to reduced lung function) are incomplete. Our results highlight the high degree of pleiotropy in complex disease-associated loci and provide novel insights into the mechanisms underlying COPD.

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy.

BACKGROUND: Global gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context. METHODS: We used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age. RESULTS: The SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10- 9), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10- 8). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells. CONCLUSION: Focusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids.