From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

Clinical predictors of pulmonary tuberculosis among South African adults with HIV.

Background: Tuberculosis (TB) clinical prediction rules rely on presence of symptoms, however many undiagnosed cases in the community are asymptomatic. This study aimed to explore the utility of clinical factors in predicting TB among people with HIV not seeking care. Methods: Baseline data were analysed from an observational cohort of ambulant adults with HIV in South Africa. Participants were tested for Mycobacterium tuberculosis (Mtb) sensitisation (interferon-γ release assay, IGRA) and microbiologically-confirmed prevalent pulmonary TB disease at baseline, and actively surveilled for incident TB through 15 months. Multivariable LASSO regression with post-selection inference was used to test associations with Mtb sensitisation and TB disease. Findings: Between March 22, 2017, and May 15, 2018, 861 participants were enrolled; Among 851 participants included in the analysis, 94·5% were asymptomatic and 45·9% sensitised to Mtb. TB prevalence was 2·0% at baseline and incidence 2·3/100 person-years through 15 months follow-up. Study site was associated with baseline Mtb sensitisation (p < 0·001), prevalent (p < 0·001), and incident TB disease (p = 0·037). Independent of site, higher CD4 counts (per 50 cells/mm3, aOR 1·48, 95%CI 1·12-1·77, p = 0·006) were associated with increased IGRA positivity, and participants without TB disease (aOR 0·80, 95%CI 0·69-0·94, p = 0·006) had reduced IGRA positivity; no variables were independently associated with prevalent TB. Mixed ancestry (aHR 1·49, 95%CI 1·30->1000, p = 0·005) and antiretroviral initiation (aHR 1·48, 95%CI 1·01-929·93, p = 0·023) were independently associated with incident TB. Models incorporating clinical features alone performed poorly in diagnosing prevalent (AUC 0·65, 95%CI 0·44-0·85) or predicting progression to incident (0·67, 0·46-0·88) TB. Interpretation: CD4 count and antiretroviral initiation, proxies for immune status and HIV stage, were associated with Mtb sensitisation and TB disease. Inadequate performance of clinical prediction models may reflect predominantly subclinical disease diagnosed in this setting and unmeasured local site factors affecting transmission and progression. Funding: The CORTIS-HR study was funded by the Bill & Melinda Gates Foundation (OPP1151915) and by the Strategic Health Innovation Partnerships Unit of the South African Medical Research Council with funds received from the South African Department of Science and Technology. The regulatory sponsor was the University of Cape Town.

Heritability jointly explained by host genotype and microbiome: will improve traits prediction?

As we observe the $70$th anniversary of the publication by Robertson that formalized the notion of 'heritability', geneticists remain puzzled by the problem of missing/hidden heritability, where heritability estimates from genome-wide association studies (GWASs) fall short of that from twin-based studies. Many possible explanations have been offered for this discrepancy, including existence of genetic variants poorly captured by existing arrays, dominance, epistasis and unaccounted-for environmental factors; albeit these remain controversial. We believe a substantial part of this problem could be solved or better understood by incorporating the host's microbiota information in the GWAS model for heritability estimation and may also increase human traits prediction for clinical utility. This is because, despite empirical observations such as (i) the intimate role of the microbiome in many complex human phenotypes, (ii) the overlap between genetic variants associated with both microbiome attributes and complex diseases and (iii) the existence of heritable bacterial taxa, current GWAS models for heritability estimate do not take into account the contributory role of the microbiome. Furthermore, heritability estimate from twin-based studies does not discern microbiome component of the observed total phenotypic variance. Here, we summarize the concept of heritability in GWAS and microbiome-wide association studies, focusing on its estimation, from a statistical genetics perspective. We then discuss a possible statistical method to incorporate the microbiome in the estimation of heritability in host GWAS.

Simulation of African and non-African low and high coverage whole genome sequence data to assess variant calling approaches.

Current variant calling (VC) approaches have been designed to leverage populations of long-range haplotypes and were benchmarked using populations of European descent, whereas most genetic diversity is found in non-European such as Africa populations. Working with these genetically diverse populations, VC tools may produce false positive and false negative results, which may produce misleading conclusions in prioritization of mutations, clinical relevancy and actionability of genes. The most prominent question is which tool or pipeline has a high rate of sensitivity and precision when analysing African data with either low or high sequence coverage, given the high genetic diversity and heterogeneity of this data. Here, a total of 100 synthetic Whole Genome Sequencing (WGS) samples, mimicking the genetics profile of African and European subjects for different specific coverage levels (high/low), have been generated to assess the performance of nine different VC tools on these contrasting datasets. The performances of these tools were assessed in false positive and false negative call rates by comparing the simulated golden variants to the variants identified by each VC tool. Combining our results on sensitivity and positive predictive value (PPV), VarDict [PPV = 0.999 and Matthews correlation coefficient (MCC) = 0.832] and BCFtools (PPV = 0.999 and MCC = 0.813) perform best when using African population data on high and low coverage data. Overall, current VC tools produce high false positive and false negative rates when analysing African compared with European data. This highlights the need for development of VC approaches with high sensitivity and precision tailored for populations characterized by high genetic variations and low linkage disequilibrium.

Dissecting genome-wide studies for microbiome-related metabolic diseases.

Despite the meteoric rise in genome-wide association studies for metabolic diseases (MetD) over the last few years, our understanding of the pathogenesis of these diseases is still far from complete. Recent developments have established that MetD arises from complex interactions between host genetics, the gut microbiome and the environment. However, our knowledge of the genetic and microbiome components involved and the underlying molecular mechanisms remains limited. Here, we review and summarize recent studies investigating the genetic and microbiome basis of MetD. Then, given the critical importance of study-individual's ancestry in these studies, we leverage 4932 whole-genome sequence samples from 18 worldwide ethnic groups to examine genetic diversity in currently reported variants associated with MetD. The analyses show marked differences in gene-specific proportion of pathogenic single-nucleotide polymorphisms (SNPs) and gene-specific SNPs MAFs across ethnic groups, highlighting the importance of population- and ethnic-specific investigations in pinpointing the causative factors for MetD. We conclude with a discussion of research areas where further investigation on interactions between host genetics, microbiome and the environment is needed.

Dating admixture events is unsolved problem in multi-way admixed populations.

Advances in human sequencing technologies, coupled with statistical and computational tools, have fostered the development of methods for dating admixture events. These methods have merits and drawbacks in estimating admixture events in multi-way admixed populations. Here, we first provide a comprehensive review and comparison of current methods pertinent to dating admixture events. Second, we assess various admixture dating tools. We do so by performing various simulations. Third, we apply the top two assessed methods to real data of a uniquely admixed population from South Africa. Results reveal that current dating admixture models are not sufficiently equipped to estimate ancient admixtures events and to identify multi-faceted admixture events in complex multi-way admixed populations. We conclude with a discussion of research areas where further work on dating admixture-based methods is needed.

A broad survey of DNA sequence data simulation tools.

In silico DNA sequence generation is a powerful technology to evaluate and validate bioinformatics tools, and accordingly more than 35 DNA sequence simulation tools have been developed. With such a diverse array of tools to choose from, an important question is: Which tool should be used for a desired outcome? This question is largely unanswered as documentation for many of these DNA simulation tools is sparse. To address this, we performed a review of DNA sequence simulation tools developed to date and evaluated 20 state-of-art DNA sequence simulation tools on their ability to produce accurate reads based on their implemented sequence error model. We provide a succinct description of each tool and suggest which tool is most appropriate for the given different scenarios. Given the multitude of similar yet non-identical tools, researchers can use this review as a guide to inform their choice of DNA sequence simulation tool. This paves the way towards assessing existing tools in a unified framework, as well as enabling different simulation scenario analysis within the same framework.

Tantalizing dilemma in risk prediction from disease scoring statistics.

Over the past decade, human host genome-wide association studies (GWASs) have contributed greatly to our understanding of the impact of host genetics on phenotypes. Recently, the microbiome has been recognized as a complex trait in host genetic variation, leading to microbiome GWAS (mGWASs). For these, many different statistical methods and software tools have been developed for association mapping. Applications of these methods and tools have revealed several important findings; however, the establishment of causal factors and the direction of causality in the interactive role between human genetic polymorphisms, the microbiome and the host phenotypes are still a huge challenge. Here, we review disease scoring approaches in host and mGWAS and their underlying statistical methods and tools. We highlight the challenges in pinpointing the genetic-associated causal factors in host and mGWAS and discuss the role of multi-omic approach in disease scoring statistics that may provide a better understanding of human phenotypic variation by enabling further system biological experiment to establish causality.

Host and Microbiome Genome-Wide Association Studies: Current State and Challenges.

The involvement of the microbiome in health and disease is well established. Microbiome genome-wide association studies (mGWAS) are used to elucidate the interaction of host genetic variation with the microbiome. The emergence of this relatively new field has been facilitated by the advent of next generation sequencing technologies that enable the investigation of the complex interaction between host genetics and microbial communities. In this paper, we review recent studies investigating host-microbiome interactions using mGWAS. Additionally, we highlight the marked disparity in the sampling population of mGWAS carried out to date and draw attention to the critical need for inclusion of diverse populations.