The power of meta-analysis in genome-wide association studies

Genome-wide association analysis using multiple Atlantic salmon populations

BACKGROUND

In a previous study, we found low persistence of linkage disequilibrium (LD) phase across breeding populations of Atlantic salmon. Accordingly, we observed no increase in accuracy from combining these populations for genomic prediction. In this study, we aimed to examine if the same were true for detection power in genome-wide association studies (GWAS), in terms of reduction in p-values, and if the precision of mapping quantitative trait loci (QTL) would improve from such analysis. Since individual records may not always be available, e.g. due to proprietorship or confidentiality, we also compared mega-analysis and meta-analysis. Mega-analysis needs access to all individual records, whereas meta-analysis utilizes parameters, such as p-values or allele substitution effects, from multiple studies or populations. Furthermore, different methods for determining the presence or absence of independent or secondary signals, such as conditional association analysis, approximate conditional and joint analysis (COJO), and the clumping approach, were assessed.

RESULTS

Mega-analysis resulted in increased detection power, in terms of reduction in p-values, and increased precision, compared to the within-population GWAS. Only one QTL was detected using conditional association analysis, both within populations and in mega-analysis, while the number of QTL detected with COJO and the clumping approach ranged from 1 to 19. The allele substitution effect and -log10p-values obtained from mega-analysis were highly correlated with the corresponding values from various meta-analysis methods. Compared to mega-analysis, a higher detection power and reduced precision were obtained with the meta-analysis methods.

CONCLUSIONS

Our results show that combining multiple datasets or populations in a mega-analysis can increase detection power and mapping precision. With meta-analysis, a higher detection power was obtained compared to mega-analysis. However, care must be taken in the interpretation of the meta-analysis results from multiple populations because their test statistics might be inflated due to population structure or cryptic relatedness.

Meta-analysis of set-based multiple phenotype association test based on GWAS summary statistics from different cohorts

Genome-wide association studies (GWAS) have emerged as popular tools for identifying genetic variants that are associated with complex diseases. Standard analysis of a GWAS involves assessing the association between each variant and a disease. However, this approach suffers from limited reproducibility and difficulties in detecting multi-variant and pleiotropic effects. Although joint analysis of multiple phenotypes for GWAS can identify and interpret pleiotropic loci which are essential to understand pleiotropy in diseases and complex traits, most of the multiple phenotype association tests are designed for a single variant, resulting in much lower power, especially when their effect sizes are small and only their cumulative effect is associated with multiple phenotypes. To overcome these limitations, set-based multiple phenotype association tests have been developed to enhance statistical power and facilitate the identification and interpretation of pleiotropic regions. In this research, we propose a new method, named Meta-TOW-S, which conducts joint association tests between multiple phenotypes and a set of variants (such as variants in a gene) utilizing GWAS summary statistics from different cohorts. Our approach applies the set-based method that Tests for the effect of an Optimal Weighted combination of variants in a gene (TOW) and accounts for sample size differences across GWAS cohorts by employing the Cauchy combination method. Meta-TOW-S combines the advantages of set-based tests and multi-phenotype association tests, exhibiting computational efficiency and enabling analysis across multiple phenotypes while accommodating overlapping samples from different GWAS cohorts. To assess the performance of Meta-TOW-S, we develop a phenotype simulator package that encompasses a comprehensive simulation scheme capable of modeling multiple phenotypes and multiple variants, including noise structures and diverse correlation patterns among phenotypes. Simulation studies validate that Meta-TOW-S maintains a desirable Type I error rate. Further simulation under different scenarios shows that Meta-TOW-S can improve power compared with other existing meta-analysis methods. When applied to four psychiatric disorders summary data, Meta-TOW-S detects a greater number of significant genes.

Emergent methods for inactivation of Cronobacter sakazakii in foods: A systematic review and meta-analysis

Cronobacter sakazakii is a potentially pathogenic bacterium that is resistant to osmotic stress and low aw, and capable of persisting in a desiccated state in powdered infant milks. It is widespread in the environment and present in various products. Despite the low incidence of cases, its high mortality rates of 40 to 80 % amongst neonates make it a microorganism of public health interest. This current study performed a comparative assessment between current reduction methods applied for C. sakazakii in various food matrices, indicating tendencies and relevant parameters for process optimization. A systematic review and meta-analysis were conducted, qualitatively identifying the main methods of inactivation and control, and quantitatively evaluating the effect of treatment factors on the reduction response. Hierarchical clustering dendrograms led to conclusions on the efficiency of each treatment. Review of recent research trend identified a focus on the potential use of alternative treatments, with most studies related to non-thermal methods and dairy products. Using random-effects meta-analysis, a summary effect-size of 4-log was estimated; however, thermal methods and treatments on dairy matrices displayed wider dispersions - of τ2 = 8.1, compared with τ2 = 4.5 for vegetal matrices and τ2 = 4.0 for biofilms. Meta-analytical models indicated that factors such as chemical concentration, energy applied, and treatment time had a more significant impact on reduction than the increase in temperature. Non-thermal treatments, synergically associated with heat, and treatments on dairy matrices were found to be the most efficient.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Empirically adjusted fixed-effects meta-analysis methods in genomic studies

In recent years, meta-analyzing summary results from multiple studies has become a common practice in genomic research, leading to a significant improvement in the power of statistical detection compared to an individual genomic study. Meta analysis methods that combine statistical estimates across studies are known to be statistically more powerful than those combining statistical significance measures. An approach combining effect size estimates based on a fixed-effects model, called METAL, has gained extreme popularity to perform the former type of meta-analysis. In this article, we discuss the limitations of METAL due to its dependence on the theoretical null distribution, leading to incorrect significance testing results. Through various simulation studies and real genomic data application, we show how modifying the z-scores in METAL, using an empirical null distribution, can significantly improve the results, especially in presence of hidden confounders. For the estimation of the null distribution, we consider two different approaches, and we highlight the scenarios when one null estimation approach outperforms the other. This article will allow researchers to gain an insight into the importance of using an empirical null distribution in the fixed-effects meta-analysis as well as in choosing the appropriate empirical null distribution estimation approach.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Federated generalized linear mixed models for collaborative genome-wide association studies

Federated association testing is a powerful approach to conduct large-scale association studies where sites share intermediate statistics through a central server. There are, however, several standing challenges. Confounding factors like population stratification should be carefully modeled across sites. In addition, it is crucial to consider disease etiology using flexible models to prevent biases. Privacy protections for participants pose another significant challenge. Here, we propose distributed Mixed Effects Genome-wide Association study (dMEGA), a method that enables federated generalized linear mixed model-based association testing across multiple sites without explicitly sharing genotype and phenotype data. dMEGA employs a reference projection to correct for population-stratification and utilizes efficient local-gradient updates among sites, incorporating both fixed and random effects. The accuracy and efficiency of dMEGA are demonstrated through simulated and real datasets. dMEGA is publicly available at https://github.com/Li-Wentao/dMEGA.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Accuracy of haplotype estimation and whole genome imputation affects complex trait analyses in complex biobanks

Sample recruitment for research consortia, biobanks, and personal genomics companies span years, necessitating genotyping in batches, using different technologies. As marker content on genotyping arrays varies, integrating such datasets is non-trivial and its impact on haplotype estimation (phasing) and whole genome imputation, necessary steps for complex trait analysis, remains under-evaluated. Using the iPSYCH dataset, comprising 130,438 individuals, genotyped in two stages, on different arrays, we evaluated phasing and imputation performance across multiple phasing methods and data integration protocols. While phasing accuracy varied by choice of method and data integration protocol, imputation accuracy varied mostly between data integration protocols. We demonstrate an attenuation in imputation accuracy within samples of non-European origin, highlighting challenges to studying complex traits in diverse populations. Finally, imputation errors can bias association tests, reduce predictive utility of polygenic scores. Carefully optimized data integration strategies enhance accuracy and replicability of complex trait analyses in complex biobanks.

Synthesizing genome regulation data with vote-counting

The increasing availability of high-throughput datasets allows amalgamating research information across a large body of genome regulation studies. Given the recent success of meta-analyses on transcriptional regulators, epigenetic marks, and enhancer:gene associations, we expect that such surveys will continue to provide novel and reproducible insights. However, meta-analyses are severely hampered by the diversity of available data, concurring protocols, an eclectic amount of bioinformatics tools, and myriads of conceivable parameter combinations. Such factors can easily bar life scientists from synthesizing omics data and substantially curb their interpretability. Despite statistical challenges of the method, we would like to emphasize the advantages of joining data from different sources through vote-counting and showcase examples that achieve a simple but highly intuitive data integration.

Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs

Growth and fat deposition are complex traits, which can affect economical income in the pig industry. Due to the intensive artificial selection, a significant genetic improvement has been observed for growth and fat deposition in pigs. Here, we first investigated genomic-wide association studies (GWAS) and population genomics (e.g., selection signature) to explore the genetic basis of such complex traits in two Large White pig lines (n = 3,727) with the GeneSeek GGP Porcine HD array (n = 50,915 SNPs). Ten genetic variants were identified to be associated with growth and fatness traits in two Large White pig lines from different genetic backgrounds by performing both within-population GWAS and cross-population GWAS analyses. These ten significant loci represented eight candidate genes, i.e., NRG4, BATF3, IRS2, ANO1, ANO9, RNF152, KCNQ5, and EYA2. One of them, ANO1 gene was simultaneously identified for both two lines in BF100 trait. Compared to single-population GWAS, cross-population GWAS was less effective for identifying SNPs with population-specific effect, but more powerful for detecting SNPs with population-shared effects. We further detected genomic regions specifically selected in each of two populations, but did not observe a significant enrichment for the heritability of growth and backfat traits in such regions. In summary, the candidate genes will provide an insight into the understanding of the genetic architecture of growth-related traits and backfat thickness, and may have a potential use in the genomic breeding programs in pigs.

An Integrated Bioinformatics Approach to Identify Network-Derived Hub Genes in Starving Zebrafish

The present study was aimed at identifying causative hub genes within modules formed by co-expression and protein-protein interaction (PPI) networks, followed by Bayesian network (BN) construction in the liver transcriptome of starved zebrafish. To this end, the GSE11107 and GSE112272 datasets from the GEO databases were downloaded and meta-analyzed using the MetaDE package, an add-on R package. Differentially expressed genes (DEGs) were identified based upon expression intensity N(µ = 0.2, σ² = 0.4). Reconstruction of BNs was performed by the bnlearn R package on genes within modules using STRINGdb and CEMiTool. ndufs5 (shared among PPI, BN and COEX), rps26, rpl10, sdhc (shared between PPI and BN), ndufa6, ndufa10, ndufb8 (shared between PPI and COEX), skp1, atp5h, ndufb10, rpl5b, zgc:193613, zgc:123327, zgc:123178, wu:fc58f10, zgc:111986, wu:fc37b12, taldo1, wu:fb62f08, zgc:64133 and acp5a (shared between COEX and BN) were identified as causative hub genes affecting gene expression in the liver of starving zebrafish. Future work will shed light on using integrative analyses of miRNA and DNA microarrays simultaneously, and performing in silico and experimental validation of these hub-causative (CST) genes affecting starvation in zebrafish.

Reproducible acquisition, management and meta-analysis of nucleotide sequence (meta)data using q2-fondue

MOTIVATION

The volume of public nucleotide sequence data has blossomed over the past two decades and is ripe for re- and meta-analyses to enable novel discoveries. However, reproducible re-use and management of sequence datasets and associated metadata remain critical challenges. We created the open source Python package q2-fondue to enable user-friendly acquisition, re-use and management of public sequence (meta)data while adhering to open data principles.

RESULTS

q2-fondue allows fully provenance-tracked programmatic access to and management of data from the NCBI Sequence Read Archive (SRA). Unlike other packages allowing download of sequence data from the SRA, q2-fondue enables full data provenance tracking from data download to final visualization, integrates with the QIIME 2 ecosystem, prevents data loss upon space exhaustion and allows download of (meta)data given a publication library. To highlight its manifold capabilities, we present executable demonstrations using publicly available amplicon, whole genome and metagenome datasets.

AVAILABILITY AND IMPLEMENTATION

q2-fondue is available as an open-source BSD-3-licensed Python package at https://github.com/bokulich-lab/q2-fondue. Usage tutorials are available in the same repository. All Jupyter notebooks used in this article are available under https://github.com/bokulich-lab/q2-fondue-examples.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The UK Biobank: A Shining Example of Genome-Wide Association Study Science with the Power to Detect the Murky Complications of Real-World Epidemiology

Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants that are reliably associated with human traits. Although GWASs are restricted to certain variant frequencies, they have improved our understanding of the genetic architecture of complex traits and diseases. The UK Biobank (UKBB) has brought substantial analytical opportunity and performance to association studies. The dramatic expansion of many GWAS sample sizes afforded by the inclusion of UKBB data has improved the power of estimation of effect sizes but, critically, has done so in a context where phenotypic depth and precision enable outcome dissection and the application of epidemiological approaches. However, at the same time, the availability of such a large, well-curated, and deeply measured population-based collection has the capacity to increase our exposure to the many complications and inferential complexities associated with GWASs and other analyses. In this review, we discuss the impact that UKBB has had in the GWAS era, some of the opportunities that it brings, and exemplar challenges that illustrate the reality of using data from this world-leading resource.

MSPJ: Discovering potential biomarkers in small gene expression datasets via ensemble learning

In transcriptomics, differentially expressed genes (DEGs) provide fine-grained phenotypic resolution for comparisons between groups and insights into molecular mechanisms underlying the pathogenesis of complex diseases or phenotypes. The robust detection of DEGs from large datasets is well-established. However, owing to various limitations (e.g., the low availability of samples for some diseases or limited research funding), small sample size is frequently used in experiments. Therefore, methods to screen reliable and stable features are urgently needed for analyses with limited sample size. In this study, MSPJ, a new machine learning approach for identifying DEGs was proposed to mitigate the reduced power and improve the stability of DEG identification in small gene expression datasets. This ensemble learning-based method consists of three algorithms: an improved multiple random sampling with meta-analysis, SVM-RFE (support vector machines-recursive feature elimination), and permutation test. MSPJ was compared with ten classical methods by 94 simulated datasets and large-scale benchmarking with 165 real datasets. The results showed that, among these methods MSPJ had the best performance in most small gene expression datasets, especially those with sample size below 30. In summary, the MSPJ method enables effective feature selection for robust DEG identification in small transcriptome datasets and is expected to expand research on the molecular mechanisms underlying complex diseases or phenotypes.

Robust meta-analysis for large-scale genomic experiments based on an empirical approach

Background

Recent high-throughput technologies have opened avenues for simultaneous analyses of thousands of genes. With the availability of a multitude of public databases, one can easily access multiple genomic study results where each study comprises of significance testing results of thousands of genes. Researchers currently tend to combine this genomic information from these multiple studies in the form of a meta-analysis. As the number of genes involved is very large, the classical meta-analysis approaches need to be updated to acknowledge this large-scale aspect of the data.

Methods

In this article, we discuss how application of standard theoretical null distributional assumptions of the classical meta-analysis methods, such as Fisher’s p-value combination and Stouffer’s Z, can lead to incorrect significant testing results, and we propose a robust meta-analysis method that empirically modifies the individual test statistics and p-values before combining them.

Results

Our proposed meta-analysis method performs best in significance testing among several meta-analysis approaches, especially in presence of hidden confounders, as shown through a wide variety of simulation studies and real genomic data analysis.

Conclusion

The proposed meta-analysis method produces superior meta-analysis results compared to the standard p-value combination approaches for large-scale simultaneous testing in genomic experiments. This is particularly useful in studies with large number of genes where the standard meta-analysis approaches can result in gross false discoveries due to the presence of unobserved confounding variables.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12874-022-01530-y.

Meta-Analysis for Epigenome-Wide Association Studies

With the rapid development of methylation profiling technology, many datasets are generated to quantify genome-wide methylation patterns. Given the heavy burden of multiple testing of hundreds of thousands of DNA methylation markers, individual studies often have limited sample sizes and power. The EWAS meta-analysis is an approach that combines results from multiple studies on the same scientific question. It helps to improve statistical power by combining information from individual studies and reduce the chances of false positives. This chapter introduces commonly used meta-analysis methods and analytical tools with application to EWAS data.

Genome-wide association study detected novel susceptibility genes for social cognition impairment in people with schizophrenia

OBJECTIVES

People with schizophrenia (SCZ) present serious and generalised deficits in social cognition (SC), which affect negatively patients' functioning and treatment outcomes. The genetic background of SC has been investigated in disorders other than SCZ providing weak and sparse results. Thus, our aim was to explore possible genetic correlates of SC dysfunctions in SCZ patients with a genome-wide study (GWAS) approach.

METHODS

We performed a GWAS meta-analysis of data coming from two cohorts made of 242 and 160 SCZ patients, respectively. SC was assessed with different tools in order to cover its different domains.

RESULTS

We found GWAS significant association between the TMEM74 gene and the patients' ability in social inference as assessed by The Awareness of Social Inference Test; this association was confirmed by both SNP-based analysis (lead SNP rs3019332 p-value = 5.24 × 10^-9) and gene-based analysis (p-value = 1.09 × 10^-7). Moreover, suggestive associations of other genes with different dimensions of SC were also found.

CONCLUSIONS

Our study shows for the first time GWAS significant or suggestive associations of some gene variants with SC domains in people with SCZ. These findings should stimulate further studies to characterise the genetic underpinning of SC dysfunctions in SCZ.

An Analysis of Postharvest Fungal Pathogens Reveals Temporal-Spatial and Host-Pathogen Associations with Fungicide Resistance-Related Mutations

Fungicides are the primary tools to control a wide range of postharvest fungal pathogens. Fungicide resistance is a widespread problem that has reduced the efficacy of fungicides. Resistance to FRAC-1 (Fungicide Resistance Action Committee-1) chemistries is associated with mutations in amino acid position 198 in the β-tubulin gene. In our study, we conducted a meta-analysis of β-tubulin sequences to infer temporal, spatial, plant host, and pathogen genus patterns of fungicide resistance in postharvest fungal pathogens. In total, data were acquired from 2,647 specimens from 12 genera of fungal phytopathogens residing in 53 countries on >200 hosts collected between 1926 and 2020. The specimens containing a position 198 mutation were globally distributed in a variety of pathosystems. Analyses showed that there are associations among the mutation and the year an isolate was collected, the pathogen genus, the pathogen host, and the collection region. Interestingly, fungicide-resistant β-tubulin genotypes have been in a decline since their peak between 2005 and 2009. FRAC-1 fungicide usage data followed a similar pattern in that applications have been in a decline since their peak between 1997 and 2003. The data show that, with the reduction of selection pressure, FRAC-1 fungicide resistance in fungal populations will decline within 5 to 10 years. Based on this line of evidence, we contend that a β-tubulin position 198 mutation has uncharacterized fitness cost(s) on fungi in nature. The compiled dataset can inform end users on the regions and hosts that are most prone to contain resistant pathogens and assist decisions concerning fungicide resistance management strategies.

Variance-component-based meta-analysis of gene-environment interactions for rare variants

Complex diseases are often caused by interplay between genetic and environmental factors. Existing gene-environment interaction (G × E) tests for rare variants largely focus on detecting gene-based G × E effects in a single study; thus, their statistical power is limited by the sample size of the study. Meta-analysis methods that synthesize summary statistics of G × E effects from multiple studies for rare variants are still limited. Based on variance component models, we propose four meta-analysis methods of testing G × E effects for rare variants: HOM-INT-FIX, HET-INT-FIX, HOM-INT-RAN, and HET-INT-RAN. Our methods consider homogeneous or heterogeneous G × E effects across studies and treat the main genetic effect as either fixed or random. Through simulations, we show that the empirical distributions of the four meta-statistics under the null hypothesis align with their expected theoretical distributions. When the interaction effect is homogeneous across studies, HOM-INT-FIX and HOM-INT-RAN have as much statistical power as a pooled analysis conducted on a single interaction test with individual-level data from all studies. When the interaction effect is heterogeneous across studies, HET-INT-FIX and HET-INT-RAN provide higher power than pooled analysis. Our methods are further validated via testing 12 candidate gene-age interactions in blood pressure traits using whole-exome sequencing data from UK Biobank.

Advancing drug discovery using the power of the human genome

Human genetics plays an increasingly important role in drug development and population health. Here we review the history of human genetics in the context of accelerating the discovery of therapies, present examples of how human genetics evidence supports successful drug targets, and discuss how polygenic risk scores could be beneficial in various clinical settings. We highlight the value of direct-to-consumer platforms in the era of fast-paced big data biotechnology, and how diverse genetic and health data can benefit society. © 2021 23andMe, Inc. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Meta-analysis of genome-wide association studies reveal common loci controlling agronomic and quality traits in a wide range of normal and heat stressed environments

Several stable QTL were detected using metaGWAS analysis for different agronomic and quality traits under 26 normal and heat stressed environments. Heat stress, exacerbated by global warming, has a negative influence on wheat production worldwide and climate resilient cultivars can help mitigate these impacts. Selection decisions should therefore depend on multi-environment experiments representing a range of temperatures at critical stages of development. Here, we applied a meta-genome wide association analysis (metaGWAS) approach to detect stable QTL with significant effects across multiple environments. The metaGWAS was applied to 11 traits scored in 26 trials that were sown at optimal or late times of sowing (TOS1 and TOS2, respectively) at five locations. A total of 2571 unique wheat genotypes (13,959 genotypes across all environments) were included and the analysis conducted on TOS1, TOS2 and both times of sowing combined (TOS1&2). The germplasm was genotyped using a 90 k Infinium chip and imputed to exome sequence level, resulting in 341,195 single nucleotide polymorphisms (SNPs). The average accuracy across all imputed SNPs was high (92.4%). The three metaGWAS analyses revealed 107 QTL for the 11 traits, of which 16 were detected in all three analyses and 23 were detected in TOS1&2 only. The remaining QTL were detected in either TOS1 or TOS2 with or without TOS1&2, reflecting the complex interactions between the environments and the detected QTL. Eight QTL were associated with grain yield and seven with multiple traits. The identified QTL provide an important resource for gene enrichment and fine mapping to further understand the mechanisms of gene × environment interaction under both heat stressed and unstressed conditions.

Integrative High Dimensional Multiple Testing with Heterogeneity under Data Sharing Constraints

Identifying informative predictors in a high dimensional regression model is a critical step for association analysis and predictive modeling. Signal detection in the high dimensional setting often fails due to the limited sample size. One approach to improving power is through meta-analyzing multiple studies which address the same scientific question. However, integrative analysis of high dimensional data from multiple studies is challenging in the presence of between-study heterogeneity. The challenge is even more pronounced with additional data sharing constraints under which only summary data can be shared across different sites. In this paper, we propose a novel data shielding integrative large-scale testing (DSILT) approach to signal detection allowing between-study heterogeneity and not requiring the sharing of individual level data. Assuming the underlying high dimensional regression models of the data differ across studies yet share similar support, the proposed method incorporates proper integrative estimation and debiasing procedures to construct test statistics for the overall effects of specific covariates. We also develop a multiple testing procedure to identify significant effects while controlling the false discovery rate (FDR) and false discovery proportion (FDP). Theoretical comparisons of the new testing procedure with the ideal individual-level meta-analysis (ILMA) approach and other distributed inference methods are investigated. Simulation studies demonstrate that the proposed testing procedure performs well in both controlling false discovery and attaining power. The new method is applied to a real example detecting interaction effects of the genetic variants for statins and obesity on the risk for type II diabetes.

Genome-wide association study and Mendelian randomization analysis provide insights for improving rice yield potential

Rice yield per plant has a complex genetic architecture, which is mainly determined by its three component traits: the number of grains per panicle (GPP), kilo-grain weight (KGW), and tillers per plant (TP). Exploring ideotype breeding based on selection for genetically less complex component traits is an alternative route for further improving rice production. To understand the genetic basis of the relationship between rice yield and component traits, we investigated the four traits of two rice hybrid populations (575 + 1495 F₁) in different environments and conducted meta-analyses of genome-wide association study (meta-GWAS). In total, 3589 significant loci for three components traits were detected, while only 3 loci for yield were detected. It indicated that rice yield is mainly controlled by minor-effect loci and hardly to be identified. Selecting quantitative trait locus/gene affected component traits to further enhance yield is recommended. Mendelian randomization design is adopted to investigate the genetic effects of loci on yield through component traits and estimate the genetic relationship between rice yield and its component traits by these loci. The loci for GPP or TP mainly had a positive genetic effect on yield, but the loci for KGW with different direction effects (positive effect or negative effect). Additionally, TP (Beta = 1.865) has a greater effect on yield than KGW (Beta = 1.016) and GPP (Beta = 0.086). Five significant loci for component traits that had an indirect effect on yield were identified. Pyramiding superior alleles of the five loci revealed improved yield. A combination of direct and indirect effects may better contribute to the yield potential of rice. Our findings provided a rationale for using component traits as indirect indices to enhanced rice yield, which will be helpful for further understanding the genetic basis of yield and provide valuable information for improving rice yield potential.

Mendelian Randomization in Cardiovascular Research: Establishing Causality When There Are Unmeasured Confounders

Mendelian randomization is an epidemiological approach to making causal inferences using observational data. It makes use of the natural randomization that occurs in the generation of an individual's genetic makeup in a way that is analogous to the study design of a randomized controlled trial and uses instrumental variable analysis where the genetic variant(s) are the instrument (analogous to random allocation to treatment group in an randomized controlled trial). As with any instrumental variable, there are 3 assumptions that must be made about the genetic instrument: (1) it is associated (not necessarily causally) with the exposure (relevance condition); (2) it is associated with the outcome only through the exposure (exclusion restriction condition); and (3) it does not share a common cause with the outcome (ie, no confounders of the genetic instrument and outcome, independence condition). Using the example of type II diabetes and coronary artery disease, we demonstrate how the method may be used to investigate causality and discuss potential benefits and pitfalls. We conclude that although Mendelian randomization studies can usually not establish causality on their own, they may usefully contribute to the evidence base and increase our certainty about the effectiveness (or otherwise) of interventions to reduce cardiovascular disease.

The Effects of Genetic Background for Diurnal Preference on Sleep Development in Early Childhood

PURPOSE

No previous research has examined the impact of the genetic background of diurnal preference on children´s sleep. Here, we examined the effects of genetic risk score for the liability of diurnal preference on sleep development in early childhood in two population-based cohorts from Finland.

PARTICIPANTS AND METHODS

The primary sample (CHILD-SLEEP, CS) comprised 1420 infants (695 girls), and the replication sample (FinnBrain, FB; 962 girls) 2063 infants. Parent-reported sleep duration, sleep-onset latency and bedtime were assessed at three, eight, 18 and 24 months in CS, and at six, 12 and 24 months in FB. Actigraphy-based sleep latency and efficiency were measured in CS in 365 infants at eight months (168 girls), and in 197 infants at 24 months (82 girls). Mean standard scores for each sleep domain were calculated in both samples. Polygenic risk scores (PRS) were used to quantitate the genetic risk for eveningness (PRSBestFit) and morningness (PRS10kBest).

RESULTS

PRSBestFit associated with longer sleep-onset latency and later bedtime, and PRS10kBest related to shorter sleep-onset latency in CS. The link between genetic risk for diurnal preference and sleep-onset latency was replicated in FB, and meta-analysis resulted in associations (P<0.0005) with both PRS-values (PRSBestFit: Z=3.55; and PRS10kBest: Z=-3.68). Finally, PRSBestFit was related to actigraphy-based lower sleep efficiency and longer sleep latency at eight months.

CONCLUSION

Genetic liability to diurnal preference for eveningness relates to longer sleep-onset during the first two years of life, and to objectively measured lowered sleep efficiency. These findings enhance our understanding on the biological factors affecting sleep development, and contribute to clarify the physiological sleep architecture in early childhood.

Genetic variation in the body mass index of adult survivors of childhood acute lymphoblastic leukemia: A report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort

BACKGROUND

Treatment characteristics such as cranial radiation therapy (CRT) do not fully explain adiposity risk in childhood acute lymphoblastic leukemia (ALL) survivors. This study was aimed at characterizing genetic variation related to adult body mass index (BMI) among survivors of childhood ALL.

METHODS

Genetic associations of BMI among 1458 adult survivors of childhood ALL (median time from diagnosis, 20 years) were analyzed by multiple approaches. A 2-stage genome-wide association study in the Childhood Cancer Survivor Study (CCSS) and the St. Jude Lifetime Cohort Study (SJLIFE) was performed. BMI was a highly polygenic trait in the general population. Within the known loci, the BMI percent variance explained was estimated, and additive interactions (chi-square test) with CRT in the CCSS were evaluated. The role of DNA methylation in CRT interaction was further evaluated in a subsample of ALL survivors.

RESULTS

In a meta-analysis of the CCSS and SJLIFE, 2 novel loci associated with adult BMI among survivors of childhood ALL (LINC00856 rs575792008 and EMR1 rs62123082; P_Meta < 5E-8) were identified. It was estimated that the more than 700 known loci explained 6.2% of the variation in adult BMI in childhood ALL survivors. Within the known loci, significant main effects for 23 loci and statistical interactions with CRT at 9 loci (P < 7.0E-5) were further identified. At 2 CRT-interacting loci, DNA methylation patterns may have differed by age.

CONCLUSIONS

Adult survivors of childhood ALL have genetic heritability for BMI similar to that observed in the general population. This study provides evidence that treatment with CRT can modify the effect of genetic variants on adult BMI in childhood ALL survivors.

Critical Analysis of Genome-Wide Association Studies: Triple Negative Breast Cancer Quae Exempli Causa

Genome-wide association studies (GWAS) are useful in assessing and analyzing either differences or variations in DNA sequences across the human genome to detect genetic risk factors of diseases prevalent within a target population under study. The ultimate goal of GWAS is to predict either disease risk or disease progression by identifying genetic risk factors. These risk factors will define the biological basis of disease susceptibility for the purposes of developing innovative, preventative, and therapeutic strategies. As single nucleotide polymorphisms (SNPs) are often used in GWAS, their relevance for triple negative breast cancer (TNBC) will be assessed in this review. Furthermore, as there are different levels and patterns of linkage disequilibrium (LD) present within different human subpopulations, a plausible strategy to evaluate known SNPs associated with incidence of breast cancer in ethnically different patient cohorts will be presented and discussed. Additionally, a description of GWAS for TNBC will be presented, involving various identified SNPs correlated with miRNA sites to determine their efficacies on either prognosis or progression of TNBC in patients. Although GWAS have identified multiple common breast cancer susceptibility variants that individually would result in minor risks, it is their combined effects that would likely result in major risks. Thus, one approach to quantify synergistic effects of such common variants is to utilize polygenic risk scores. Therefore, studies utilizing predictive risk scores (PRSs) based on known breast cancer susceptibility SNPs will be evaluated. Such PRSs are potentially useful in improving stratification for screening, particularly when combining family history, other risk factors, and risk prediction models. In conclusion, although interpretation of the results from GWAS remains a challenge, the use of SNPs associated with TNBC may elucidate and better contextualize these studies.

Genome-wide association study of emotional empathy in children

The genetic contribution to different aspects of empathy is now established, although the exact loci are unknown. We undertook a genome-wide association study of emotional empathy (EE) as measured by emotion recognition skills in 4,780 8-year old children from the ALSPAC cohort who were genotyped and imputed to Phase 1 version 3 of the 1000 Genomes Project. We failed to find any genome-wide significant signal in either our unstratified analysis or analysis stratified according to sex. A gene-based association analysis similarly failed to find any significant loci. In contrast, our transcriptome-wide association study (TWAS) with a whole blood reference panel identified two significant loci in the unstratified analysis, residualised for the effects of age, sex and IQ. One signal was for CD93 on chromosome 20; this gene is not strongly expressed in the brain, however. The other signal was for AL118508, a non-protein coding pseudogene, which completely lies within CD93’s genomic coordinates, thereby explaining its signal. Neither are obvious candidates for involvement in the brain processes that underlie emotion recognition and its developmental pathways.

Flaws (and quality) in research today: can artificial intelligence intervene?

The existing flaws in both conducting and reporting of research have been outlined and criticized in the past. Weak research design, poor methodology, lack of fresh ideas and poor reporting are the main points to blame. Issues have been continually raised on the types of results published, review process, sponsorship, notion, ethics, and incentives in publishing, the role of regulatory agencies and stakeholders, the role of funding, and the cooperation between funders and academic institutions and the training of both clinicians and methodologists or statisticians. As a result, there is loss of the utmost goal: the production of robust research to form recommendations to support pragmatic decision in a real-world context. We propose the construction of a model based on artificial intelligence that could assist stakeholders, clinicians, and patients to guide conducting the best quality of research. We briefly describe the levels of the workflow, including the input and output data collection, the feature extraction/selection, the architecture, and parameterization of the model, along with its training, operation, and refinement.

Systematic Identification of Housekeeping Genes Possibly Used as References in Caenorhabditis elegans by Large-Scale Data Integration

For accurate gene expression quantification, normalization of gene expression data against reliable reference genes is required. It is known that the expression levels of commonly used reference genes vary considerably under different experimental conditions, and therefore, their use for data normalization is limited. In this study, an unbiased identification of reference genes in Caenorhabditis elegans was performed based on 145 microarray datasets (2296 gene array samples) covering different developmental stages, different tissues, drug treatments, lifestyle, and various stresses. As a result, thirteen housekeeping genes (rps-23, rps-26, rps-27, rps-16, rps-2, rps-4, rps-17, rpl-24.1, rpl-27, rpl-33, rpl-36, rpl-35, and rpl-15) with enhanced stability were comprehensively identified by using six popular normalization algorithms and RankAggreg method. Functional enrichment analysis revealed that these genes were significantly overrepresented in GO terms or KEGG pathways related to ribosomes. Validation analysis using recently published datasets revealed that the expressions of newly identified candidate reference genes were more stable than the commonly used reference genes. Based on the results, we recommended using rpl-33 and rps-26 as the optimal reference genes for microarray and rps-2 and rps-4 for RNA-sequencing data validation. More importantly, the most stable rps-23 should be a promising reference gene for both data types. This study, for the first time, successfully displays a large-scale microarray data driven genome-wide identification of stable reference genes for normalizing gene expression data and provides a potential guideline on the selection of universal internal reference genes in C. elegans, for quantitative gene expression analysis.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

A novel estimator of between-study variance in random-effects models

BACKGROUND

With the rapid development of high-throughput sequencing technologies, many datasets on the same biological subject are generated. A meta-analysis is an approach that combines results from different studies on the same topic. The random-effects model in a meta-analysis enables the modeling of differences between studies by incorporating the between-study variance.

RESULTS

This paper proposes a moments estimator of the between-study variance that represents the across-study variation. A new random-effects method (DSLD2), which involves two-step estimation starting with the DSL estimate and the [Formula: see text] in the second step, is presented. The DSLD2 method is compared with 6 other meta-analysis methods based on effect sizes across 8 aspects under three hypothesis settings. The results show that DSLD2 is a suitable method for identifying differentially expressed genes under the first hypothesis. The DSLD2 method is also applied to Alzheimer's microarray datasets. The differentially expressed genes detected by the DSLD2 method are significantly enriched in neurological diseases.

CONCLUSIONS

The results from both simulationes and an application show that DSLD2 is a suitable method for detecting differentially expressed genes under the first hypothesis.

Genetically Determined Height and Risk of Non-hodgkin Lymphoma

Although the evidence is not consistent, epidemiologic studies have suggested that taller adult height may be associated with an increased risk of some non-Hodgkin lymphoma (NHL) subtypes. Height is largely determined by genetic factors, but how these genetic factors may contribute to NHL risk is unknown. We investigated the relationship between genetic determinants of height and NHL risk using data from eight genome-wide association studies (GWAS) comprising 10,629 NHL cases, including 3,857 diffuse large B-cell lymphoma (DLBCL), 2,847 follicular lymphoma (FL), 3,100 chronic lymphocytic leukemia (CLL), and 825 marginal zone lymphoma (MZL) cases, and 9,505 controls of European ancestry. We evaluated genetically predicted height by constructing polygenic risk scores using 833 height-associated SNPs. We used logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) for association between genetically determined height and the risk of four NHL subtypes in each GWAS and then used fixed-effect meta-analysis to combine subtype results across studies. We found suggestive evidence between taller genetically determined height and increased CLL risk (OR = 1.08, 95% CI = 1.00-1.17, p = 0.049), which was slightly stronger among women (OR = 1.15, 95% CI: 1.01-1.31, p = 0.036). No significant associations were observed with DLBCL, FL, or MZL. Our findings suggest that there may be some shared genetic factors between CLL and height, but other endogenous or environmental factors may underlie reported epidemiologic height associations with other subtypes.

Methodological challenges in constructing DNA methylation risk scores

Polygenic approaches often access more variance of complex traits than is possible by single variant approaches. For genotype data, genetic risk scores (GRS) are widely used for risk prediction as well as in association and interaction studies. Recently, interest has been growing in transferring GRS approaches to DNA methylation data (methylation risk scores, MRS), which can be used 1) as biomarkers for environmental exposures, 2) in association analyses in which single CpG sites do not achieve significance, 3) as dimension reduction approach in interaction and mediation analyses, and 4) to predict individual risks of disease or treatment success. Most GRS approaches can directly be transferred to methylation data. However, since methylation data is more sensitive to confounding, e.g. by age and tissue, it is more complex to find appropriate external weights. In this review, we will outline the adaption of current GRS approaches to methylation data and highlight occurring challenges.

There is no association between rs6296 and alcoholism: a meta-analysis

Previous studies have reported controversial results about the association between rs6296 and alcoholism. Thus, a meta-analysis was performed to further explore this association. A comprehensive search was conducted to identify relevant case-control or cohort studies (up to December 1, 2017). A fixed- or random-effect model was selected as a pooling method depending on the heterogeneity among studies. The heterogeneity was measured by Q test and I² statistic. The Harbord and Peters test was used to estimate publication bias. Fifteen English articles with 16 outcomes and 5,429 participants were included in this meta-analysis. A fixed-effect model was chosen, and the pooled result showed that rs6296 was not related to alcoholism (z = 1.93, p = .053). The Harbord and Peters test showed that there was no publication bias. This meta-analysis indicated that rs6296 may be not be significantly associated with alcoholism, which needs to be further confirmed by future research.

Meta-MultiSKAT: Multiple phenotype meta-analysis for region-based association test

The power of genetic association analyses can be increased by jointly meta-analyzing multiple correlated phenotypes. Here, we develop a meta-analysis framework, Meta-MultiSKAT, that uses summary statistics to test for association between multiple continuous phenotypes and variants in a region of interest. Our approach models the heterogeneity of effects between studies through a kernel matrix and performs a variance component test for association. Using a genotype kernel, our approach can test for rare-variants and the combined effects of both common and rare-variants. To achieve robust power, within Meta-MultiSKAT, we developed fast and accurate omnibus tests combining different models of genetic effects, functional genomic annotations, multiple correlated phenotypes, and heterogeneity across studies. In addition, Meta-MultiSKAT accommodates situations where studies do not share exactly the same set of phenotypes or have differing correlation patterns among the phenotypes. Simulation studies confirm that Meta-MultiSKAT can maintain the type-I error rate at the exome-wide level of 2.5 × 10^-6 . Further simulations under different models of association show that Meta-MultiSKAT can improve the power of detection from 23% to 38% on average over single phenotype-based meta-analysis approaches. We demonstrate the utility and improved power of Meta-MultiSKAT in the meta-analyses of four white blood cell subtype traits from the Michigan Genomics Initiative (MGI) and SardiNIA studies.

Scientific Utopia III: Crowdsourcing Science

Most scientific research is conducted by small teams of investigators who together formulate hypotheses, collect data, conduct analyses, and report novel findings. These teams operate independently as vertically integrated silos. Here we argue that scientific research that is horizontally distributed can provide substantial complementary value, aiming to maximize available resources, promote inclusiveness and transparency, and increase rigor and reliability. This alternative approach enables researchers to tackle ambitious projects that would not be possible under the standard model. Crowdsourced scientific initiatives vary in the degree of communication between project members from largely independent work curated by a coordination team to crowd collaboration on shared activities. The potential benefits and challenges of large-scale collaboration span the entire research process: ideation, study design, data collection, data analysis, reporting, and peer review. Complementing traditional small science with crowdsourced approaches can accelerate the progress of science and improve the quality of scientific research.

Meta-analysis of genome-wide association studies provides insights into genetic control of tomato flavor

Tomato flavor has changed over the course of long-term domestication and intensive breeding. To understand the genetic control of flavor, we report the meta-analysis of genome-wide association studies (GWAS) using 775 tomato accessions and 2,316,117 SNPs from three GWAS panels. We discover 305 significant associations for the contents of sugars, acids, amino acids, and flavor-related volatiles. We demonstrate that fruit citrate and malate contents have been impacted by selection during domestication and improvement, while sugar content has undergone less stringent selection. We suggest that it may be possible to significantly increase volatiles that positively contribute to consumer preferences while reducing unpleasant volatiles, by selection of the relevant allele combinations. Our results provide genetic insights into the influence of human selection on tomato flavor and demonstrate the benefits obtained from meta-analysis.

Flavor is one of the most important traits for improving tomato sensory quality and consumer acceptability. Here, the authors report meta-analysis of genome-wide association studies of flavor related traits and show genetic insights into the influence of human selection during domestication and improvement.

Identifying Multi-Omics Causers and Causal Pathways for Complex Traits

The central dogma of molecular biology delineates a unidirectional causal flow, i.e., DNA → RNA → protein → trait. Genome-wide association studies, next-generation sequencing association studies, and their meta-analyses have successfully identified ~12,000 susceptibility genetic variants that are associated with a broad array of human physiological traits. However, such conventional association studies ignore the mediate causers (i.e., RNA, protein) and the unidirectional causal pathway. Such studies may not be ideally powerful; and the genetic variants identified may not necessarily be genuine causal variants. In this article, we model the central dogma by a mediate causal model and analytically prove that the more remote an omics level is from a physiological trait, the smaller the magnitude of their correlation is. Under both random and extreme sampling schemes, we numerically demonstrate that the proteome-trait correlation test is more powerful than the transcriptome-trait correlation test, which in turn is more powerful than the genotype-trait association test. In conclusion, integrating RNA and protein expressions with DNA data and causal inference are necessary to gain a full understanding of how genetic causal variants contribute to phenotype variations.

GWAS of Behavioral Traits

Over the past decade, genome-wide association studies (GWAS) have evolved into a powerful tool to investigate genetic risk factors for human diseases via a hypothesis-free scan of the genome. The success of GWAS for psychiatric disorders and behavioral traits have been somewhat mixed, partly owing to the complexity and heterogeneity of these traits. Significant progress has been made in the last few years in the development and implementation of complex statistical methods and algorithms incorporating GWAS. Such advanced statistical methods applied to GWAS hits in combination with incorporation of different layers of genomics data have catapulted the search for novel genes for behavioral traits and improved our understanding of the complex polygenic architecture of these traits.This chapter will give a brief overview on GWAS and statistical methods currently used in GWAS. The chapter will focus on reviewing the current literature and highlight some of the most important GWAS on psychiatric and other behavioral traits and will conclude with a discussion on future directions.

Variability of Two Metabolomic Platforms in CKD

BACKGROUND AND OBJECTIVES

Nontargeted metabolomics can measure thousands of low-molecular-weight biochemicals, but important gaps limit its utility for biomarker discovery in CKD. These include the need to characterize technical and intraperson analyte variation, to pool data across platforms, and to outline analyte relationships with eGFR.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

Plasma samples from 49 individuals with CKD (eGFR<60 ml/min per 1.73 m² and/or ≥1 g proteinuria) were examined from two study visits; 20 samples were repeated as blind replicates. To enable comparison across two nontargeted platforms, samples were profiled at Metabolon and the Broad Institute.

RESULTS

The Metabolon platform reported 837 known metabolites and 483 unnamed compounds (selected from 44,953 unknown ion features). The Broad Institute platform reported 594 known metabolites and 26,106 unknown ion features. Median coefficients of variation (CVs) across blind replicates were 14.6% (Metabolon) and 6.3% (Broad Institute) for known metabolites, and 18.9% for (Metabolon) unnamed compounds and 24.5% for (Broad Institute) unknown ion features. Median CVs for day-to-day variability were 29.0% (Metabolon) and 24.9% (Broad Institute) for known metabolites, and 41.8% for (Metabolon) unnamed compounds and 40.9% for (Broad Institute) unknown ion features. A total of 381 known metabolites were shared across platforms (median correlation 0.89). Many metabolites were negatively correlated with eGFR at P<0.05, including 35.7% (Metabolon) and 18.9% (Broad Institute) of known metabolites.

CONCLUSIONS

Nontargeted metabolomics quantifies >1000 analytes with low technical CVs, and agreement for overlapping metabolites across two leading platforms is excellent. Many metabolites demonstrate substantial intraperson variation and correlation with eGFR.

GWAS for BMI: a treasure trove of fundamental insights into the genetic basis of obesity

Muller et al. [1] have provided a strong critique of the Genome-Wide Association Studies (GWAS) of body-mass index (BMI), arguing that the GWAS approach for the study of BMI is flawed, and has provided us with few biological insights. They suggest that what is needed instead is a new start, involving GWAS for more complex energy balance related traits. In this invited counter-point, we highlight the substantial advances that have occurred in the obesity field, directly stimulated by the GWAS of BMI. We agree that GWAS for BMI is not perfect, but consider that the best route forward for additional discoveries will likely be to expand the search for common and rare variants linked to BMI and other easily obtained measures of obesity, rather than attempting to perform new, much smaller GWAS for energy balance traits that are complex and expensive to measure. For GWAS in general, we emphasise that the power from increasing the sample size of a crude but easily measured phenotype outweighs the benefits of better phenotyping.

Biomedical informatics and machine learning for clinical genomics

While tens of thousands of pathogenic variants are used to inform the many clinical applications of genomics, there remains limited information on quantitative disease risk for the majority of variants used in clinical practice. At the same time, rising demand for genetic counselling has prompted a growing need for computational approaches that can help interpret genetic variation. Such tasks include predicting variant pathogenicity and identifying variants that are too common to be penetrant. To address these challenges, researchers are increasingly turning to integrative informatics approaches. These approaches often leverage vast sources of data, including electronic health records and population-level allele frequency databases (e.g. gnomAD), as well as machine learning techniques such as support vector machines and deep learning. In this review, we highlight recent informatics and machine learning approaches that are improving our understanding of pathogenic variation and discuss obstacles that may limit their emerging role in clinical genomics.

Genome-wide association study identifies the GLDC/IL33 locus associated with survival of osteosarcoma patients

Survival rates for osteosarcoma, the most common primary bone cancer, have changed little over the past three decades and are particularly low for patients with metastatic disease. We conducted a multi-institutional genome-wide association study (GWAS) to identify germline genetic variants associated with overall survival in 632 patients with osteosarcoma, including 523 patients of European ancestry and 109 from Brazil. We conducted a time-to-event analysis and estimated hazard ratios (HR) and 95% confidence intervals (CI) using Cox proportional hazards models, with and without adjustment for metastatic disease. The results were combined across the European and Brazilian case sets using a random-effects meta-analysis. The strongest association after meta-analysis was for rs3765555 at 9p24.1, which was inversely associated with overall survival (HR = 1.76; 95% CI 1.41-2.18, p = 4.84 × 10-7 ). After imputation across this region, the combined analysis identified two SNPs that reached genome-wide significance. The strongest single association was with rs55933544 (HR = 1.9; 95% CI 1.5-2.4; p = 1.3 × 10-8 ), which localizes to the GLDC gene, adjacent to the IL33 gene and was consistent across both the European and Brazilian case sets. Using publicly available data, the risk allele was associated with lower expression of IL33 and low expression of IL33 was associated with poor survival in an independent set of patients with osteosarcoma. In conclusion, we have identified the GLDC/IL33 locus on chromosome 9p24.1 as associated with overall survival in patients with osteosarcoma. Further studies are needed to confirm this association and shed light on the biological underpinnings of this susceptibility locus.

Beyond genomics: understanding exposotypes through metabolomics

BACKGROUND

Over the past 20 years, advances in genomic technology have enabled unparalleled access to the information contained within the human genome. However, the multiple genetic variants associated with various diseases typically account for only a small fraction of the disease risk. This may be due to the multifactorial nature of disease mechanisms, the strong impact of the environment, and the complexity of gene-environment interactions. Metabolomics is the quantification of small molecules produced by metabolic processes within a biological sample. Metabolomics datasets contain a wealth of information that reflect the disease state and are consequent to both genetic variation and environment. Thus, metabolomics is being widely adopted for epidemiologic research to identify disease risk traits. In this review, we discuss the evolution and challenges of metabolomics in epidemiologic research, particularly for assessing environmental exposures and providing insights into gene-environment interactions, and mechanism of biological impact.

MAIN TEXT

Metabolomics can be used to measure the complex global modulating effect that an exposure event has on an individual phenotype. Combining information derived from all levels of protein synthesis and subsequent enzymatic action on metabolite production can reveal the individual exposotype. We discuss some of the methodological and statistical challenges in dealing with this type of high-dimensional data, such as the impact of study design, analytical biases, and biological variance. We show examples of disease risk inference from metabolic traits using metabolome-wide association studies. We also evaluate how these studies may drive precision medicine approaches, and pharmacogenomics, which have up to now been inefficient. Finally, we discuss how to promote transparency and open science to improve reproducibility and credibility in metabolomics.

CONCLUSIONS

Comparison of exposotypes at the human population level may help understanding how environmental exposures affect biology at the systems level to determine cause, effect, and susceptibilities. Juxtaposition and integration of genomics and metabolomics information may offer additional insights. Clinical utility of this information for single individuals and populations has yet to be routinely demonstrated, but hopefully, recent advances to improve the robustness of large-scale metabolomics will facilitate clinical translation.

Runs of homozygosity: windows into population history and trait architecture

Long runs of homozygosity (ROH) arise when identical haplotypes are inherited from each parent and thus a long tract of genotypes is homozygous. Cousin marriage or inbreeding gives rise to such autozygosity; however, genome-wide data reveal that ROH are universally common in human genomes even among outbred individuals. The number and length of ROH reflect individual demographic history, while the homozygosity burden can be used to investigate the genetic architecture of complex disease. We discuss how to identify ROH in genome-wide microarray and sequence data, their distribution in human populations and their application to the understanding of inbreeding depression and disease risk.