Population genetics and GWAS: A primer

BMI Interacts with the Genome to Regulate Gene Expression Globally, with Emphasis in the Brain and Gut

Genome-wide association studies identify common genomic variants associated with disease across a population. Individual environmental effects are often not included, despite evidence that environment mediates genomic regulation of higher order biology. Body mass index (BMI) is associated with complex disorders across clinical specialties, yet has not been modeled as a genomic environment. Here, we tested for expression quantitative trait (eQTL) loci that contextually regulate gene expression across the BMI spectrum using an interaction approach. We parsed the impact of cell type, enhancer interactions, and created novel BMI-dynamic gene expression predictor models. We found that BMI main effects associated with endocrine gene expression, while interactive variant-by-BMI effects impacted gene expression in the brain and gut. Cortical BMI-dynamic loci were experimentally dysregulated by inflammatory cytokines in an in vitro system. Using BMI-dynamic models, we identify novel genes in nitric oxide signaling pathways in the nucleus accumbens significantly associated with depression and smoking. While neither genetics nor BMI are sufficient as standalone measures to capture the complexity of downstream cellular consequences, including environment powers disease gene discovery.

Ancestry-specific high-risk gene variant profiling unmasks diabetes-associated genes

How ancestry-associated genetic variance affects disparities in the risk of polygenic diseases and influences the identification of disease-associated genes warrants a deeper understanding. We hypothesized that the discovery of genes associated with polygenic diseases may be limited by the overreliance on single-nucleotide polymorphism (SNP)-based genomic investigation, as most significant variants identified in genome-wide SNP association studies map to introns and intergenic regions of the genome. To overcome such potential limitations, we developed a gene-constrained, function-based analytical method centered on high-risk variants (hrV) that encode frameshifts, stopgains or splice site disruption. We analyzed the total number of hrV per gene in populations of different ancestry, representing a total of 185 934 subjects. Using this analysis, we developed a quantitative index of hrV (hrVI) across 20 428 genes within each population. We then applied hrVI analysis to the discovery of genes associated with type 2 diabetes mellitus (T2DM), a polygenic disease with ancestry-related disparity. HrVI profiling and gene-to-gene comparisons of ancestry-specific hrV between the case (20 781 subjects) and control (24 440 subjects) populations in the T2DM national repository identified 57 genes associated with T2DM, 40 of which were discoverable only by ancestry-specific analysis. These results illustrate how a function-based, ancestry-specific analysis of genetic variations can accelerate the identification of genes associated with polygenic diseases. Besides T2DM, such analysis may facilitate our understanding of the genetic basis for other polygenic diseases that are also greatly influenced by environmental and behavioral factors, such as obesity, hypertension and Alzheimer's disease.

The association of GNB5 with Alzheimer disease revealed by genomic analysis restricted to variants impacting gene function

Disease-associated variants identified from genome-wide association studies (GWASs) frequently map to non-coding areas of the genome such as introns and intergenic regions. An exclusive reliance on gene-agnostic methods of genomic investigation could limit the identification of relevant genes associated with polygenic diseases such as Alzheimer disease (AD). To overcome such potential restriction, we developed a gene-constrained analytical method that considers only moderate- and high-risk variants that affect gene coding sequences. We report here the application of this approach to publicly available datasets containing 181,388 individuals without and with AD and the resulting identification of 660 genes potentially linked to the higher AD prevalence among Africans/African Americans. By integration with transcriptome analysis of 23 brain regions from 2,728 AD case-control samples, we concentrated on nine genes that potentially enhance the risk of AD: AACS, GNB5, GNS, HIPK3, MED13, SHC2, SLC22A5, VPS35, and ZNF398. GNB5, the fifth member of the heterotrimeric G protein beta family encoding Gβ5, is primarily expressed in neurons and is essential for normal neuronal development in mouse brain. Homozygous or compound heterozygous loss of function of GNB5 in humans has previously been associated with a syndrome of developmental delay, cognitive impairment, and cardiac arrhythmia. In validation experiments, we confirmed that Gnb5 heterozygosity enhanced the formation of both amyloid plaques and neurofibrillary tangles in the brains of AD model mice. These results suggest that gene-constrained analysis can complement the power of GWASs in the identification of AD-associated genes and may be more broadly applicable to other polygenic diseases.

Pitfalls and potential of high-throughput plant phenotyping platforms

Automated high-throughput plant phenotyping (HTPP) enables non-invasive, fast and standardized evaluations of a large number of plants for size, development, and certain physiological variables. Many research groups recognize the potential of HTPP and have made significant investments in HTPP infrastructure, or are considering doing so. To make optimal use of limited resources, it is important to plan and use these facilities prudently and to interpret the results carefully. Here we present a number of points that users should consider before purchasing, building or utilizing such equipment. They relate to (1) the financial and time investment for acquisition, operation, and maintenance, (2) the constraints associated with such machines in terms of flexibility and growth conditions, (3) the pros and cons of frequent non-destructive measurements, (4) the level of information provided by proxy traits, and (5) the utilization of calibration curves. Using data from an Arabidopsis experiment, we demonstrate how diurnal changes in leaf angle can impact plant size estimates from top-view cameras, causing deviations of more than 20% over the day. Growth analysis data from another rosette species showed that there was a curvilinear relationship between total and projected leaf area. Neglecting this curvilinearity resulted in linear calibration curves that, although having a high r2 (> 0.92), also exhibited large relative errors. Another important consideration we discussed is the frequency at which calibration curves need to be generated and whether different treatments, seasons, or genotypes require distinct calibration curves. In conclusion, HTPP systems have become a valuable addition to the toolbox of plant biologists, provided that these systems are tailored to the research questions of interest, and users are aware of both the possible pitfalls and potential involved.

Balancing selection on genomic deletion polymorphisms in humans

A key question in biology is why genomic variation persists in a population for extended periods. Recent studies have identified examples of genomic deletions that have remained polymorphic in the human lineage for hundreds of millennia, ostensibly owing to balancing selection. Nevertheless, genome-wide investigation of ancient and possibly adaptive deletions remains an imperative exercise. Here, we demonstrate an excess of polymorphisms in present-day humans that predate the modern human-Neanderthal split (ancient polymorphisms), which cannot be explained solely by selectively neutral scenarios. We analyze the adaptive mechanisms that underlie this excess in deletion polymorphisms. Using a previously published measure of balancing selection, we show that this excess of ancient deletions is largely owing to balancing selection. Based on the absence of signatures of overdominance, we conclude that it is a rare mode of balancing selection among ancient deletions. Instead, more complex scenarios involving spatially and temporally variable selective pressures are likely more common mechanisms. Our results suggest that balancing selection resulted in ancient deletions harboring disproportionately more exonic variants with GWAS (genome-wide association studies) associations. We further found that ancient deletions are significantly enriched for traits related to metabolism and immunity. As a by-product of our analysis, we show that deletions are, on average, more deleterious than single nucleotide variants. We can now argue that not only is a vast majority of common variants shared among human populations, but a considerable portion of biologically relevant variants has been segregating among our ancestors for hundreds of thousands, if not millions, of years.

Deriving alpha angle from anterior-posterior dual-energy x-ray absorptiometry scans: an automated and validated approach

Introduction: Alpha angle (AA) is a widely used imaging measure of hip shape that is commonly used to define cam morphology, a bulging of the lateral aspect of the femoral head. Cam morphology has shown strong associations with hip osteoarthritis (OA) making the AA a clinically relevant measure. In both clinical practice and research studies, AA tends to be measured manually which can be inconsistent and time-consuming. Objective: We aimed to (i) develop an automated method of deriving AA from anterior-posterior dual-energy x-ray absorptiometry (DXA) scans; and (ii) validate this method against manual measures of AA. Methods: 6,807 individuals with left hip DXAs were selected from UK Biobank. Outline points were manually placed around the femoral head on 1,930 images before training a Random Forest-based algorithm to place the points on a further 4,877 images. An automatic method for calculating AA was written in Python 3 utilising these outline points. An iterative approach was taken to developing and validating the method, testing the automated measures against independent batches of manually measured images in sequential experiments. Results: Over the course of six experimental stages the concordance correlation coefficient, when comparing the automatic AA to manual measures of AA, improved from 0.28 [95% confidence interval 0.13-0.43] for the initial version to 0.88 [0.84-0.92] for the final version. The inter-rater kappa statistic comparing automatic versus manual measures of cam morphology, defined as AA ³≥60°, improved from 0.43 [80% agreement] for the initial version to 0.86 [94% agreement] for the final version. Conclusions: We have developed and validated an automated measure of AA from DXA scans, showing high agreement with manually measuring AA. The proposed method is available to the wider research community from Zenodo.

Evolutionary Ecology of Plant-Arthropod Interactions in Light of the "Omics" Sciences: A Broad Guide

Aboveground plant-arthropod interactions are typically complex, involving herbivores, predators, pollinators, and various other guilds that can strongly affect plant fitness, directly or indirectly, and individually, synergistically, or antagonistically. However, little is known about how ongoing natural selection by these interacting guilds shapes the evolution of plants, i.e., how they affect the differential survival and reproduction of genotypes due to differences in phenotypes in an environment. Recent technological advances, including next-generation sequencing, metabolomics, and gene-editing technologies along with traditional experimental approaches (e.g., quantitative genetics experiments), have enabled far more comprehensive exploration of the genes and traits involved in complex ecological interactions. Connecting different levels of biological organization (genes to communities) will enhance the understanding of evolutionary interactions in complex communities, but this requires a multidisciplinary approach. Here, we review traditional and modern methods and concepts, then highlight future avenues for studying the evolution of plant-arthropod interactions (e.g., plant-herbivore-pollinator interactions). Besides promoting a fundamental understanding of plant-associated arthropod communities' genetic background and evolution, such knowledge can also help address many current global environmental challenges.

Preparation and Curation of Omics Data for Genome-Wide Association Studies

With the development of large-scale molecular phenotyping platforms, genome-wide association studies have greatly developed, being no longer limited to the analysis of classical agronomic traits, such as yield or flowering time, but also embracing the dissection of the genetic basis of molecular traits. Data generated by omics platforms, however, pose some technical and statistical challenges to the classical methodology and assumptions of an association study. Although genotyping data are subject to strict filtering procedures, and several advanced statistical approaches are now available to adjust for population structure, less attention has been instead devoted to the preparation of omics data prior to GWAS. In the present chapter, we briefly present the methods to acquire profiling data from transcripts, proteins, and small molecules, and discuss the tools and possibilities to clean, normalize, and remove the unwanted variation from large datasets of molecular phenotypic traits prior to their use in GWAS.

Meta-analysis: Congruence of genomic and phenotypic differentiation across diverse natural study systems

Linking genotype to phenotype is a primary goal for understanding the genomic underpinnings of evolution. However, little work has explored whether patterns of linked genomic and phenotypic differentiation are congruent across natural study systems and traits. Here, we investigate such patterns with a meta-analysis of studies examining population-level differentiation at subsets of loci and traits putatively responding to divergent selection. We show that across the 31 studies (88 natural population-level comparisons) we examined, there was a moderate (R 2 = 0.39) relationship between genomic differentiation (F ST ) and phenotypic differentiation (P ST ) for loci and traits putatively under selection. This quantitative relationship between P ST and F ST for loci under selection in diverse taxa provides broad context and cross-system predictions for genomic and phenotypic adaptation by natural selection in natural populations. This context may eventually allow for more precise ideas of what constitutes "strong" differentiation, predictions about the effect size of loci, comparisons of taxa evolving in nonparallel ways, and more. On the other hand, links between P ST and F ST within studies were very weak, suggesting that much work remains in linking genomic differentiation to phenotypic differentiation at specific phenotypes. We suggest that linking genotypes to specific phenotypes can be improved by correlating genomic and phenotypic differentiation across a spectrum of diverging populations within a taxon and including wide coverage of both genomes and phenomes.

Deriving alpha angle from anterior-posterior dual-energy x-ray absorptiometry scans: an automated and validated approach

Introduction: Alpha angle (AA) is a widely used measure of hip shape that is commonly used to define cam morphology, a bulging of the lateral aspect of the femoral head. Cam morphology has shown strong associations with hip osteoarthritis (OA) making the AA a clinically relevant measure. In both clinical practice and research studies, AA tends to be measured manually which can be inconsistent and time-consuming. Objective: We aimed to (i) develop an automated method of deriving AA from anterior-posterior dual-energy x-ray absorptiometry (DXA) scans; and (ii) validate this method against manual measures of AA. Methods: 6,807 individuals with left hip DXAs were selected from UK Biobank. Outline points were manually placed around the femoral head on 1,930 images before training a Random Forest-based algorithm to place the points on a further 4,877 images. An automatic method for calculating AA was written in Python 3 utilising these outline points. An iterative approach was taken to developing and validating the method, testing the automated measures against independent batches of manually measured images in sequential experiments. Results: Over the course of six experimental stages the concordance correlation coefficient, when comparing the automatic AA to manual measures of AA, improved from 0.28 [95% confidence interval 0.13-0.43] for the initial version to 0.88 [0.84-0.92] for the final version. The inter-rater kappa statistic comparing automatic versus manual measures of cam morphology, defined as AA ³≥60°, improved from 0.43 [80% agreement] for the initial version to 0.86 [94% agreement] for the final version. Conclusions: We have developed and validated an automated measure of AA from DXA scans, showing high agreement with manually measuring AA. The proposed method is available to the wider research community from Zenodo.

Individual Differences in Cognitive Performance Are Better Predicted by Global Rather Than Localized BOLD Activity Patterns Across the Cortex

Despite its central role in revealing the neurobiological mechanisms of behavior, neuroimaging research faces the challenge of producing reliable biomarkers for cognitive processes and clinical outcomes. Statistically significant brain regions, identified by mass univariate statistical models commonly used in neuroimaging studies, explain minimal phenotypic variation, limiting the translational utility of neuroimaging phenotypes. This is potentially due to the observation that behavioral traits are influenced by variations in neuroimaging phenotypes that are globally distributed across the cortex and are therefore not captured by thresholded, statistical parametric maps commonly reported in neuroimaging studies. Here, we developed a novel multivariate prediction method, the Bayesian polyvertex score, that turns a unthresholded statistical parametric map into a summary score that aggregates the many but small effects across the cortex for behavioral prediction. By explicitly assuming a globally distributed effect size pattern and operating on the mass univariate summary statistics, it was able to achieve higher out-of-sample variance explained than mass univariate and popular multivariate methods while still preserving the interpretability of a generative model. Our findings suggest that similar to the polygenicity observed in the field of genetics, the neural basis of complex behaviors may rest in the global patterning of effect size variation of neuroimaging phenotypes, rather than in localized, candidate brain regions and networks.

Hierarchical Modelling of Haplotype Effects on a Phylogeny

We introduce a hierarchical model to estimate haplotype effects based on phylogenetic relationships between haplotypes and their association with observed phenotypes. In a population there are many, but not all possible, distinct haplotypes and few observations per haplotype. Further, haplotype frequencies tend to vary substantially. Such data structure challenge estimation of haplotype effects. However, haplotypes often differ only due to few mutations, and leveraging similarities can improve the estimation of effects. We build on extensive literature and develop an autoregressive model of order one that models haplotype effects by leveraging phylogenetic relationships described with a directed acyclic graph. The phylogenetic relationships can be either in a form of a tree or a network, and we refer to the model as the haplotype network model. The model can be included as a component in a phenotype model to estimate associations between haplotypes and phenotypes. Our key contribution is that we obtain a sparse model, and by using hierarchical autoregression, the flow of information between similar haplotypes is estimated from the data. A simulation study shows that the hierarchical model can improve estimates of haplotype effects compared to an independent haplotype model, especially with few observations for a specific haplotype. We also compared it to a mutation model and observed comparable performance, though the haplotype model has the potential to capture background specific effects. We demonstrate the model with a study of mitochondrial haplotype effects on milk yield in cattle. We provide R code to fit the model with the INLA package.

Genome-Wide Association for Itraconazole Sensitivity in Non-resistant Clinical Isolates of Aspergillus fumigatus

Aspergillus fumigatus is a potentially lethal opportunistic pathogen that infects over ~200,000 people and causes ~100,000 deaths per year globally. Treating A. fumigatus infections is particularly challenging because of the recent emergence of azole-resistance. The majority of studies focusing on the molecular mechanisms underlying azole resistance have examined azole-resistant isolates. However, isolates that are susceptible to azoles also display variation in their sensitivity, presenting a unique opportunity to identify genes contributing to azole sensitivity. Here, we used genome-wide association (GWA) analysis to identify loci involved in azole sensitivity by analyzing the association between 68,853 SNPs and itraconazole (ITCZ) minimum inhibitory concentration (MIC) in 76 clinical isolates of A. fumigatus from Japan. Population structure analysis suggests the presence of four distinct populations, with ITCZ MICs distributed relatively evenly across populations. We independently conducted GWA when treating ITCZ MIC as a quantitative trait and a binary trait, and identified two SNPs with strong associations in both analyses. These SNPs fell within the coding regions of Afu2g02220 and Afu2g02140. We functionally validated Afu2g02220 by knocking it out using a CRISPR/Cas9 approach, because orthologs of this gene are involved in sterol modification and ITCZ targets the ergosterol biosynthesis pathway. Knockout strains displayed no difference in growth compared to the parent strain in minimal media, yet a minor but consistent inhibition of growth in the presence of 0.15 μg/ml ITCZ. Our results suggest that GWA paired with efficient gene deletion is a powerful and unbiased strategy for identifying the genetic basis of complex traits in A. fumigatus.

Population Genomics of Filamentous Plant Pathogens-A Brief Overview of Research Questions, Approaches, and Pitfalls

With ever-decreasing sequencing costs, research on the population biology of plant pathogens is transitioning from population genetics-using dozens of genetic markers or polymorphism data of several genes-to population genomics-using several hundred to tens of thousands of markers or whole-genome sequence data. The field of population genomics is characterized by rapid theoretical and methodological advances and by numerous steps and pitfalls in its technical and analytical workflow. In this article, we aim to provide a brief overview of topics relevant to the study of population genomics of filamentous plant pathogens and direct readers to more extensive reviews for in-depth understanding. We briefly discuss different types of population genomics-inspired research questions and give insights into the sampling strategies that can be used to answer such questions. We then consider different sequencing strategies, the various options available for data processing, and some of the currently available tools for population genomic data analysis. We conclude by highlighting some of the hurdles along the population genomic workflow, providing cautionary warnings relative to assumptions and technical challenges, and presenting our own future perspectives of the field of population genomics for filamentous plant pathogens.

Why Tumor Genetic Heterogeneity May Require Rethinking Cancer Genesis and Treatment

Tumor genetic heterogeneity, in which individual tumors contain both multiple variant cancer-associated and normal genes, has been widely reported, although its significance has yet to be fully understood. We propose a genetic heterogeneity-based selection-centric hypothesis in which genetic heterogeneity, caused by the temporary reduction of DNA repair efficiency, occurs very early in human development, resulting in a small minority of cells in normal tissues acquiring cancer-associated genes that remain dormant. Cancer develops when precancer cells are selected for by altered tissue microenvironments; similar scenarios occur with development of metastases and therapeutic resistance in established cancer. This suggests that a normal cell selection treatment approach based on preferentially selecting normal cells within tumors may be effective in treating cancer.

Accuracy of genomic selection and long-term genetic gain for resistance to Verticillium wilt in strawberry

Verticillium wilt, a soil-borne disease caused by the fungal pathogen Verticillium dahliae, threatens strawberry (Fragaria × ananassa) production worldwide. The development of resistant cultivars has been a persistent challenge, in part because the genetics of resistance is complex. The heritability of resistance and genetic gains in breeding for resistance to this pathogen have not been well documented. To elucidate the genetics, assess long-term genetic gains, and estimate the accuracy of genomic selection for resistance to Verticillium wilt, we analyzed a genetically diverse population of elite and exotic germplasm accessions (n = 984), including 245 cultivars developed since 1854. We observed a full range of phenotypes, from highly susceptible to highly resistant: < 3% were classified as highly resistant, whereas > 50% were classified as moderately to highly susceptible. Broad-sense heritability estimates ranged from 0.70-0.76, whereas narrow-sense genomic heritability estimates ranged from 0.33-0.45. We found that genetic gains in breeding for resistance to Verticillium wilt have been negative over the last 165 years (mean resistance has decreased over time). We identified several highly resistant accessions that might harbor favorable alleles that are either rare or non-existent in modern populations. We did not observe the segregation of large-effect loci. The accuracy of genomic predictions ranged from 0.38-0.53 among years and whole-genome regression methods. We show that genomic selection has promise for increasing genetic gains and accelerating the development of resistant cultivars in strawberry by shortening selection cycles and enabling selection in early developmental stages without phenotyping.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Optimizing the Power to Identify the Genetic Basis of Complex Traits with Evolve and Resequence Studies

Evolve and resequence (E&R) studies are frequently used to dissect the genetic basis of quantitative traits. By subjecting a population to truncating selection for several generations and estimating the allele frequency differences between selected and nonselected populations using next-generation sequencing (NGS), the loci contributing to the selected trait may be identified. The role of different parameters, such as, the population size or the number of replicate populations has been examined in previous works. However, the influence of the selection regime, that is the strength of truncating selection during the experiment, remains little explored. Using whole genome, individual based forward simulations of E&R studies, we found that the power to identify the causative alleles may be maximized by gradually increasing the strength of truncating selection during the experiment. Notably, such an optimal selection regime comes at no or little additional cost in terms of sequencing effort and experimental time. Interestingly, we also found that a selection regime which optimizes the power to identify the causative loci is not necessarily identical to a regime that maximizes the phenotypic response. Finally, our simulations suggest that an E&R study with an optimized selection regime may have a higher power to identify the genetic basis of quantitative traits than a genome-wide association study, highlighting that E&R is a powerful approach for finding the loci underlying complex traits.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Host-pathogen coevolution increases genetic variation in susceptibility to infection

It is common to find considerable genetic variation in susceptibility to infection in natural populations. We have investigated whether natural selection increases this variation by testing whether host populations show more genetic variation in susceptibility to pathogens that they naturally encounter than novel pathogens. In a large cross-infection experiment involving four species of Drosophila and four host-specific viruses, we always found greater genetic variation in susceptibility to viruses that had coevolved with their host. We went on to examine the genetic architecture of resistance in one host species, finding that there are more major-effect genetic variants in coevolved host-pathogen interactions. We conclude that selection by pathogens has increased genetic variation in host susceptibility, and much of this effect is caused by the occurrence of major-effect resistance polymorphisms within populations.

Manhattan Harvester and Cropper: a system for GWAS peak detection

Background

Selection of interesting regions from genome wide association studies (GWAS) is typically performed by eyeballing of Manhattan Plots. This is no longer possible with thousands of different phenotypes. There is a need for tools that can automatically detect genomic regions that correspond to what the experienced researcher perceives as peaks worthwhile of further study.

Results

We developed Manhattan Harvester, a tool designed for “peak extraction” from GWAS summary files and computation of parameters characterizing various aspects of individual peaks. We present the algorithms used and a model for creating a general quality score that evaluates peaks similarly to that of a human researcher. Our tool Cropper utilizes a graphical interface for inspecting, cropping and subsetting Manhattan Plot regions. Cropper is used to validate and visualize the regions detected by Manhattan Harvester.

Conclusions

We conclude that our tools fill the current void in automatically screening large number of GWAS output files in batch mode. The interesting regions are detected and quantified by various parameters by Manhattan Harvester. Cropper offers graphical tools for in-depth inspection of the regions. The tools are open source and freely available.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2600-4) contains supplementary material, which is available to authorized users.

Gene-Environment-Time Interactions in Neurodegenerative Diseases: Hypotheses and Research Approaches

Background

Amyotrophic lateral sclerosis (ALS), Alzheimer's, and Parkinson's diseases are age-related neurodegenerative diseases. ALS is not a single entity but a syndrome with many different causes. In all 3 diseases, gene mutations account for only 10-15% of cases. Many environmental and lifestyle factors have been implicated as risk factors for ALS, though none have been proven to cause the disease. It is generally believed that ALS results from interactions between environmental risk factors and genetic predisposing factors. The advent of next-generation sequencing and recent advances in research into environmental risk factors offer the opportunity to investigate these interactions.

Summary

We propose a hypothesis to explain the syndrome of ALS based on the interaction of many individual environmental risk factors with many individual genetic predisposing factors. We hypothesize that there are many such combinations of individual, specific, genetic, and environmental factors, and that each combination can lead to the development of the syndrome of ALS. We also propose a hypothesis that explains the overlap between the age-related neurodegenerations and their genetic underpinnings. Age and duration of exposure are crucial factors in these age-related neurodegenerative diseases, and we consider how these may relate to gene-environment interactions.

Key Messages

To date, genetic studies and environmental studies have investigated the causes of ALS separately. We argue that this univariate approach will not lead to discoveries of important gene-environment interactions. We propose new research approaches to investigating gene-environment interactions based on these hypotheses.