Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank

Genome-wide variance quantitative trait locus analysis suggests small interaction effects in blood pressure traits

Genome-wide variance quantitative trait loci (vQTL) analysis complements genome-wide association study (GWAS) and has the potential to identify novel variants associated with the trait, explain additional trait variance and lead to the identification of factors that modulate the genetic effects. I conducted genome-wide analysis of the UK Biobank data and identified 27 vQTLs associated with systolic blood pressure (SBP), diastolic blood pressure (DBP) and pulse pressure (PP). The top single-nucleotide polymorphisms (SNPs) are enriched for expression QTLs (eQTLs) or splicing QTLs (sQTLs) annotated by GTEx, suggesting their regulatory roles in mediating the associations with blood pressure (BP). Of the 27 vQTLs, 14 are known BP-associated QTLs discovered by GWASs. The heteroscedasticity effects of the 13 novel vQTLs are larger than their genetic main effects, which were not detected by existing GWASs. The total R-squared of the 27 top SNPs due to variance heteroscedasticity is 0.28%, compared with 0.50% owing to their main effects. The overall effect size of the variance heteroscedasticity is small in GWAS SNPs compared with their main effects. For the 411, 384 and 285 GWAS SNPs associated with SBP, DBP and PP, respectively, their heteroscedasticity effects were 0.52%, 0.43%, and 0.16%, and their main effects were 5.13%, 5.61%, and 3.75%, respectively. The number and effects of the vQTLs are small, which suggests that the effects of gene-environment and gene-gene interactions are small. The main effects of the SNPs remain the major source of genetic variance for BP, which would probably be true for other complex traits as well.

Age and diet shape the genetic architecture of body weight in diversity outbred mice

Understanding how genetic variation shapes a complex trait relies on accurately quantifying both the additive genetic and genotype–environment interaction effects in an age-dependent manner. We used a linear mixed model to quantify diet-dependent genetic contributions to body weight measured through adulthood in diversity outbred female mice under five diets. We observed that heritability of body weight declined with age under all diets, except the 40% calorie restriction diet. We identified 14 loci with age-dependent associations and 19 loci with age- and diet-dependent associations, with many diet-dependent loci previously linked to neurological function and behavior in mice or humans. We found their allelic effects to be dynamic with respect to genomic background, age, and diet, identifying several loci where distinct alleles affect body weight at different ages. These results enable us to more fully understand and predict the effectiveness of dietary intervention on overall health throughout age in distinct genetic backgrounds.

Body weight is one trait influenced by genes, age and environmental factors. Both internal and external environmental pressures are known to affect genetic variation over time. However, it is largely unknown how all factors – including age – interact to shape metabolism and bodyweight.

Wright et al. set out to quantify the interactions between genes and diet in ageing mice and found that the effect of genetics on mouse body weight changes with age. In the experiments, Wright et al. weighed 960 female mice with diverse genetic backgrounds, starting at two months of age into adulthood. The animals were randomized to different diets at six months of age. Some mice had unlimited food access, others received 20% or 40% less calories than a typical mouse diet, and some fasted one or two days per week.

Variations in their genetic background explained about 80% of differences in mice’s weight, but the influence of genetics relative to non-genetic factors decreased as they aged. Mice on the 40% calorie restriction diet were an exception to this rule and genetics accounted for 80% of their weight throughout adulthood, likely due to reduced influence from diet and reduced interactions between diet and genes. Several genes involved in metabolism, neurological function, or behavior, were associated with mouse weight.

The experiments highlight the importance of considering interactions between genetics, environment, and age in determining complex traits like body weight. The results and the approaches used by Wright et al. may help other scientists learn more about how the genetic predisposition to disease changes with environmental stimuli and age.

Variance-quantitative trait loci enable systematic discovery of gene-environment interactions for cardiometabolic serum biomarkers

Gene-environment interactions represent the modification of genetic effects by environmental exposures and are critical for understanding disease and informing personalized medicine. These often induce differential phenotypic variance across genotypes; these variance-quantitative trait loci can be prioritized in a two-stage interaction detection strategy to greatly reduce the computational and statistical burden and enable testing of a broader range of exposures. We perform genome-wide variance-quantitative trait locus analysis for 20 serum cardiometabolic biomarkers by multi-ancestry meta-analysis of 350,016 unrelated participants in the UK Biobank, identifying 182 independent locus-biomarker pairs (p < 4.5×10-9). Most are concentrated in a small subset (4%) of loci with genome-wide significant main effects, and 44% replicate (p < 0.05) in the Women's Genome Health Study (N = 23,294). Next, we test each locus-biomarker pair for interaction across 2380 exposures, identifying 847 significant interactions (p < 2.4×10-7), of which 132 are independent (p < 0.05) after accounting for correlation between exposures. Specific examples demonstrate interaction of triglyceride-associated variants with distinct body mass- versus body fat-related exposures as well as genotype-specific associations between alcohol consumption and liver stress at the ADH1B gene. Our catalog of variance-quantitative trait loci and gene-environment interactions is publicly available in an online portal.

Leveraging the local genetic structure for trans-ancestry association mapping

Over the past two decades, genome-wide association studies (GWASs) have successfully advanced our understanding of the genetic basis of complex traits. Despite the fruitful discovery of GWASs, most GWAS samples are collected from European populations, and these GWASs are often criticized for their lack of ancestry diversity. Trans-ancestry association mapping (TRAM) offers an exciting opportunity to fill the gap of disparities in genetic studies between non-Europeans and Europeans. Here, we propose a statistical method, LOG-TRAM, to leverage the local genetic architecture for TRAM. By using biobank-scale datasets, we showed that LOG-TRAM can greatly improve the statistical power of identifying risk variants in under-represented populations while producing well-calibrated p values. We applied LOG-TRAM to the GWAS summary statistics of various complex traits/diseases from BioBank Japan, UK Biobank, and African populations. We obtained substantial gains in power and achieved effective correction of confounding biases in TRAM. Finally, we showed that LOG-TRAM can be successfully applied to identify ancestry-specific loci and the LOG-TRAM output can be further used for construction of more accurate polygenic risk scores in under-represented populations.

Polygenic Scores for Plasticity: A New Tool for Studying Gene-Environment Interplay

Fertility, health, education, and other outcomes of interest to demographers are the product of an individual's genetic makeup and their social environment. Yet, gene × environment (G×E) research deploys a limited toolkit on the genetic side to study the gene-environment interplay, relying on polygenic scores (PGSs) that reflect the influence of genetics on levels of an outcome. In this article, we develop a genetic summary measure better suited for G×E research: variance polygenic scores (vPGSs), which are PGSs that reflect genetic contributions to plasticity in outcomes. First, we use the UK Biobank (N ∼ 408,000 in the analytic sample) and the Health and Retirement Study (N ∼ 5,700 in the analytic sample) to compare four approaches to constructing PGSs for plasticity. The results show that widely used methods for discovering which genetic variants affect outcome variability fail to serve as distinctive new tools for G×E. Second, using the PGSs that do capture distinctive genetic contributions to plasticity, we analyze heterogeneous effects of a UK education reform on health and educational attainment. The results show the properties of a useful new tool for population scientists studying the interplay of nature and nurture and for population-based studies that are releasing PGSs to applied researchers.

A compressed variance component mixed model for detecting QTNs and QTN-by-environment and QTN-by-QTN interactions in genome-wide association studies

Although genome-wide association studies are widely used to mine genes for quantitative traits, the effects to be estimated are confounded, and the methodologies for detecting interactions are imperfect. To address these issues, the mixed model proposed here first estimates the genotypic effects for AA, Aa, and aa, and the genotypic polygenic background replaces additive and dominance polygenic backgrounds. Then, the estimated genotypic effects are partitioned into additive and dominance effects using a one-way analysis of variance model. This strategy was further expanded to cover QTN-by-environment interactions (QEIs) and QTN-by-QTN interactions (QQIs) using the same mixed-model framework. Thus, a three-variance-component mixed model was integrated with our multi-locus random-SNP-effect mixed linear model (mrMLM) method to establish a new methodological framework, 3VmrMLM, that detects all types of loci and estimates their effects. In Monte Carlo studies, 3VmrMLM correctly detected all types of loci and almost unbiasedly estimated their effects, with high powers and accuracies and a low false positive rate. In re-analyses of 10 traits in 1439 rice hybrids, detection of 269 known genes, 45 known gene-by-environment interactions, and 20 known gene-by-gene interactions strongly validated 3VmrMLM. Further analyses of known genes showed more small (67.49%), minor-allele-frequency (35.52%), and pleiotropic (30.54%) genes, with higher repeatability across datasets (54.36%) and more dominance loci. In addition, a heteroscedasticity mixed model in multiple environments and dimension reduction methods in quite a number of environments were developed to detect QEIs, and variable selection under a polygenic background was proposed for QQI detection. This study provides a new approach for revealing the genetic architecture of quantitative traits.

A robust mean and variance test with application to high-dimensional phenotypes

Most studies of continuous health-related outcomes examine differences in mean levels (location) of the outcome by exposure. However, identifying effects on the variability (scale) of an outcome, and combining tests of mean and variability (location-and-scale), could provide additional insights into biological mechanisms. A joint test could improve power for studies of high-dimensional phenotypes, such as epigenome-wide association studies of DNA methylation at CpG sites. One possible cause of heterogeneity of variance is a variable interacting with exposure in its effect on outcome, so a joint test of mean and variability could help in the identification of effect modifiers. Here, we review a scale test, based on the Brown-Forsythe test, for analysing variability of a continuous outcome with respect to both categorical and continuous exposures, and develop a novel joint location-and-scale score (JLSsc) test. These tests were compared to alternatives in simulations and used to test associations of mean and variability of DNA methylation with gender and gestational age using data from the Accessible Resource for Integrated Epigenomics Studies (ARIES). In simulations, the Brown-Forsythe and JLSsc tests retained correct type I error rates when the outcome was not normally distributed in contrast to the other approaches tested which all had inflated type I error rates. These tests also identified > 7500 CpG sites for which either mean or variability in cord blood methylation differed according to gender or gestational age. The Brown-Forsythe test and JLSsc are robust tests that can be used to detect associations not solely driven by a mean effect.

IPRS: Leveraging Gene-Environment Interaction to Reconstruct Polygenic Risk Score

Background: Polygenic risk score (PRS) is widely regarded as a predictor of genetic susceptibility to disease, applied to individuals to predict the risk of disease occurrence. When the gene-environment (G×E) interaction is considered, the traditional PRS prediction model directly uses PRS to interact with the environment without considering the interactions between each variant and environment, which may lead to prediction performance and risk stratification of complex diseases are not promising.

Methods: We developed a method called interaction PRS (iPRS), reconstructing PRS by leveraging G×E interactions. Two extensive simulations evaluated prediction performance, risk stratification, and calibration performance of the iPRS prediction model, and compared it with the traditional PRS prediction model. Real data analysis was performed using existing data from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial study to predict genetic susceptibility, pack-years of smoking history, and G×E interactions in patients with lung cancer.

Results: Two extensive simulations indicated iPRS prediction model could improve the prediction performance of disease risk, the accuracy of risk stratification, and clinical calibration performance compared with the traditional PRS prediction model, especially when antagonism accounted for the majority of the interaction. PLCO real data analysis also suggested that the iPRS prediction model was superior to the PRS prediction model in predictive effect (p = 0.0205).

Conclusion: IPRS prediction model could have a good application prospect in predicting disease risk, optimizing the screening of high-risk populations, and improving the clinical benefits of preventive interventions among populations.

GWAS of longitudinal trajectories at biobank scale

Biobanks linked to massive, longitudinal electronic health record (EHR) data make numerous new genetic research questions feasible. One among these is the study of biomarker trajectories. For example, high blood pressure measurements over visits strongly predict stroke onset, and consistently high fasting glucose and Hb1Ac levels define diabetes. Recent research reveals that not only the mean level of biomarker trajectories but also their fluctuations, or within-subject (WS) variability, are risk factors for many diseases. Glycemic variation, for instance, is recently considered an important clinical metric in diabetes management. It is crucial to identify the genetic factors that shift the mean or alter the WS variability of a biomarker trajectory. Compared to traditional cross-sectional studies, trajectory analysis utilizes more data points and captures a complete picture of the impact of time-varying factors, including medication history and lifestyle. Currently, there are no efficient tools for genome-wide association studies (GWASs) of biomarker trajectories at the biobank scale, even for just mean effects. We propose TrajGWAS, a linear mixed effect model-based method for testing genetic effects that shift the mean or alter the WS variability of a biomarker trajectory. It is scalable to biobank data with 100,000 to 1,000,000 individuals and many longitudinal measurements and robust to distributional assumptions. Simulation studies corroborate that TrajGWAS controls the type I error rate and is powerful. Analysis of eleven biomarkers measured longitudinally and extracted from UK Biobank primary care data for more than 150,000 participants with 1,800,000 observations reveals loci that significantly alter the mean or WS variability.

Nonreciprocal Transport in a Bilayer of MnBi2Te4 and Pt

MnBi₂Te₄ (MBT) is the first intrinsic magnetic topological insulator with the interaction of spin-momentum locked surface electrons and intrinsic magnetism, and it exhibits novel magnetic and topological phenomena. Recent studies suggested that the interaction of electrons and magnetism can be affected by the Mn-doped Bi₂Te₃ phase at the surface due to inevitable structural defects. Here, we report an observation of nonreciprocal transport, that is, current-direction-dependent resistance, in a bilayer composed of antiferromagnetic MBT and nonmagnetic Pt. The emergence of the nonreciprocal response below the Néel temperature confirms a correlation between nonreciprocity and intrinsic magnetism in the surface state of MBT. The angular dependence of the nonreciprocal transport indicates that nonreciprocal response originates from the asymmetry scattering of electrons at the surface of MBT mediated by magnon. Our work provides an insight into nonreciprocity arising from the correlation between magnetism and Dirac surface electrons in intrinsic magnetic topological insulators.

Phenome-wide association study of the major histocompatibility complex region in the Korean population identifies novel association signals

Human leukocyte antigen (HLA) gene variants in the major histocompatibility complex (MHC) region are associated with numerous complex human diseases and quantitative traits. Previous phenome-wide association studies (PheWAS) for this region demonstrated that HLA association patterns to the phenome have both population-specific and population-shared components. We performed MHC PheWAS in the Korean population by analyzing associations between phenotypes and genetic variants in the MHC region using the Korea Biobank Array project data samples from the Korean Genome and Epidemiology Study (KoGES) cohorts. Using this single-population dataset, we curated and analyzed 82 phenotypes for 125 673 Korean individuals after imputing HLA using CookHLA, a recently developed imputation framework. More than one-third of these phenotypes showed significant associations, confirming 56 known associations and discovering 13 novel association signals that were not reported previously. In addition, we analyzed heritability explained by the variants in the MHC region and genetic correlations among phenotypes based on the MHC variants.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.

Differential effects of steroid hormones on levels of broad-sense heritability in a wild bird: possible mechanism of environment × genetic variance interaction?

Genetic variation is one of the key concepts in evolutionary biology and an important prerequisite of evolutionary change. However, we know very little about processes that modulate its levels in wild populations. In particular, we still are to understand why genetic variances often depend on environmental conditions. One of possible environment-sensitive modulators of observed levels of genetic variance are maternal effects. In this study we attempt to experimentally test the hypothesis that maternally transmitted agents (e.g. hormones) may influence the expression of genetic variance in quantitative traits in the offspring. We manipulated the levels of steroid hormones (testosterone and corticosterone) in eggs laid by blue tits in a wild population. Our experimental setup allowed for full crossing of genetic and rearing effects with the experimental manipulation. We observed that birds treated with corticosterone exhibited a significant decrease in broad-sense genetic variance of tarsus length, and an increase in this component in body mass on the 2nd day post-hatching. Our study indicates, that maternally transmitted substances such as hormones may have measurable impact on the levels of genetic variance and hence, on the evolutionary potential of quantitative traits.

Gene by Environment Interactions reveal new regulatory aspects of signaling network plasticity

Phenotypes can change during exposure to different environments through the regulation of signaling pathways that operate in integrated networks. How signaling networks produce different phenotypes in different settings is not fully understood. Here, Gene by Environment Interactions (GEIs) were used to explore the regulatory network that controls filamentous/invasive growth in the yeast Saccharomyces cerevisiae. GEI analysis revealed that the regulation of invasive growth is decentralized and varies extensively across environments. Different regulatory pathways were critical or dispensable depending on the environment, microenvironment, or time point tested, and the pathway that made the strongest contribution changed depending on the environment. Some regulators even showed conditional role reversals. Ranking pathways' roles across environments revealed an under-appreciated pathway (OPI1) as the single strongest regulator among the major pathways tested (RAS, RIM101, and MAPK). One mechanism that may explain the high degree of regulatory plasticity observed was conditional pathway interactions, such as conditional redundancy and conditional cross-pathway regulation. Another mechanism was that different pathways conditionally and differentially regulated gene expression, such as target genes that control separate cell adhesion mechanisms (FLO11 and SFG1). An exception to decentralized regulation of invasive growth was that morphogenetic changes (cell elongation and budding pattern) were primarily regulated by one pathway (MAPK). GEI analysis also uncovered a round-cell invasion phenotype. Our work suggests that GEI analysis is a simple and powerful approach to define the regulatory basis of complex phenotypes and may be applicable to many systems.

Large uncertainty in individual polygenic risk score estimation impacts PRS-based risk stratification

Although the cohort-level accuracy of polygenic risk scores (PRSs)-estimates of genetic value at the individual level-has been widely assessed, uncertainty in PRSs remains underexplored. In the present study, we show that Bayesian PRS methods can estimate the variance of an individual's PRS and can yield well-calibrated credible intervals via posterior sampling. For 13 real traits in the UK Biobank (n = 291,273 unrelated 'white British'), we observe large variances in individual PRS estimates which impact interpretation of PRS-based stratification; averaging across traits, only 0.8% (s.d. = 1.6%) of individuals with PRS point estimates in the top decile have corresponding 95% credible intervals fully contained in the top decile. We provide an analytical estimator for the expectation of individual PRS variance as a function of SNP heritability, number of causal SNPs and sample size. Our results showcase the importance of incorporating uncertainty in individual PRS estimates into subsequent analyses.

Appraisal of Gene-Environment Interactions in GWAS for Evidence-Based Precision Nutrition Implementation

PURPOSE OF REVIEW

This review aims to analyse the currently reported gene-environment (G × E) interactions in genome-wide association studies (GWAS), involving environmental factors such as lifestyle and dietary habits related to metabolic syndrome phenotypes. For this purpose, the present manuscript reviews the available GWAS registered on the GWAS Catalog reporting the interaction between environmental factors and metabolic syndrome traits.

RECENT FINDINGS

Advances in omics-related analytical and computational approaches in recent years have led to a better understanding of the biological processes underlying these G × E interactions. A total of 42 GWAS were analysed, reporting over 300 loci interacting with environmental factors. Alcohol consumption, sleep time, smoking habit and physical activity were the most studied environmental factors with significant G × E interactions. The implementation of more comprehensive GWAS will provide a better understanding of the metabolic processes that determine individual responses to environmental exposures and their association with the development of chronic diseases such as obesity and the metabolic syndrome. This will facilitate the development of precision approaches for better prevention, management and treatment of these diseases.

Shared genetic liability between major depressive disorder and osteoarthritis

Aims

Deciphering the genetic relationships between major depressive disorder (MDD) and osteoarthritis (OA) may facilitate an understanding of their biological mechanisms, as well as inform more effective treatment regimens. We aim to investigate the mechanisms underlying relationships between MDD and OA in the context of common genetic variations.

Methods

Linkage disequilibrium score regression was used to test the genetic correlation between MDD and OA. Polygenic analysis was performed to estimate shared genetic variations between the two diseases. Two-sample bidirectional Mendelian randomization analysis was used to investigate causal relationships between MDD and OA. Genomic loci shared between MDD and OA were identified using cross-trait meta-analysis. Fine-mapping of transcriptome-wide associations was used to prioritize putatively causal genes for the two diseases.

Results

MDD has a significant genetic correlation with OA (r_g = 0.29) and the two diseases share a considerable proportion of causal variants. Mendelian randomization analysis indicates that genetic liability to MDD has a causal effect on OA (b_xy = 0.24) and genetic liability to OA conferred a causal effect on MDD (b_xy = 0.20). Cross-trait meta-analyses identified 29 shared genomic loci between MDD and OA. Together with fine-mapping of transcriptome-wide association signals, our results suggest that Estrogen Receptor 1 (ESR1), SRY-Box Transcription Factor 5 (SOX5), and Glutathione Peroxidase 1 (GPX1) may have therapeutic implications for both MDD and OA.

Conclusion

The study reveals substantial shared genetic liability between MDD and OA, which may confer risk for one another. Our findings provide a novel insight into phenotypic relationships between MDD and OA.

Cite this article: Bone Joint Res 2022;11(1):12–22.

FunGraph: A statistical protocol to reconstruct omnigenic multilayer interactome networks for complex traits

We describe a statistical protocol of how to reconstruct and dissect functional omnigenic multilayer interactome networks that mediate complex dynamic traits in a genome-wide association study (GWAS). This protocol, named FunGraph, can analyze how each locus affects phenotypic variation through its own direct effect and a complete set of indirect effects due to regulation by other loci co-existing in large-scale networks. FunGraph is applicable to any GWAS aimed to characterize the genetic architecture of dynamic phenotypic traits.

For complete details on the use and execution of this protocol, please refer to Wang et al. (2021).

•

TurboID enabled biotin-based proximity labeling protocol for C. elegans

•

Experimental design guidelines for proximity labeling in C. elegans

•

A step-by-step TurboID protocol from transgene design to protein identification

Modeling Interaction and Dispersion Effects in the Analysis of Gene-by-Environment Interaction

Genotype-by-environment interaction (GxE) studies probe heterogeneity in response to risk factors or interventions. Popular methods for estimation of GxE examine multiplicative interactions between individual genetic and environmental measures. However, risk factors and interventions may modulate the total variance of an epidemiological outcome that itself represents the aggregation of many other etiological components. We expand the traditional GxE model to directly model genetic and environmental moderation of the dispersion of the outcome. We derive a test statistic, [Formula: see text], for inferring whether an interaction identified between individual genetic and environmental measures represents a more general pattern of moderation of the total variance in the phenotype by either the genetic or the environmental measure. We validate our method via extensive simulation, and apply it to investigate genotype-by-birth year interactions for Body Mass Index (BMI) with polygenic scores in the Health and Retirement Study (N = 11,586) and individual genetic variants in the UK Biobank (N = 380,605). We find that changes in the penetrance of a genome-wide polygenic score for BMI across birth year are partly representative of a more general pattern of expanding BMI variation across generations. Three individual variants found to be more strongly associated with BMI among later born individuals, were also associated with the magnitude of variability in BMI itself within any given birth year, suggesting that they may confer general sensitivity of BMI to a range of unmeasured factors beyond those captured by birth year. We introduce an expanded GxE regression model that explicitly models genetic and environmental moderation of the dispersion of the outcome under study. This approach can determine whether GxE interactions identified are specific to the measured predictors or represent a more general pattern of moderation of the total variance in the outcome by the genetic and environmental measures.

Detecting Heterogeneity of Intervention Effects Using Analysis and Meta-analysis of Differences in Variance Between Trial Arms

BACKGROUND

Randomized controlled trials (RCTs) with continuous outcomes usually only examine mean differences in response between trial arms. If the intervention has heterogeneous effects, then outcome variances will also differ between arms. Power of an individual trial to assess heterogeneity is lower than the power to detect the same size of main effect.

METHODS

We describe several methods for assessing differences in variance in trial arms and apply them to a single trial with individual patient data and to meta-analyses using summary data. Where individual data are available, we use regression-based methods to examine the effects of covariates on variation. We present an additional method to meta-analyze differences in variances with summary data.

RESULTS

In the single trial, there was agreement between methods, and the difference in variance was largely due to differences in prevalence of depression at baseline. In two meta-analyses, most individual trials did not show strong evidence of a difference in variance between arms, with wide confidence intervals. However, both meta-analyses showed evidence of greater variance in the control arm, and in one example, this was perhaps because mean outcome in the control arm was higher.

CONCLUSIONS

Using meta-analysis, we overcame low power of individual trials to examine differences in variance using meta-analysis. Evidence of differences in variance should be followed up to identify potential effect modifiers and explore other possible causes such as varying compliance.

Gene-Environment Interactions for Cardiovascular Disease

PURPOSE OF REVIEW

We provide an overview of recent findings with respect to gene-environment (GxE) interactions for cardiovascular disease (CVD) risk and discuss future opportunities for advancing the field.

RECENT FINDINGS

Over the last several years, GxE interactions for CVD have mostly been identified for smoking and coronary artery disease (CAD) or related risk factors. By comparison, there is more limited evidence for GxE interactions between CVD outcomes and other exposures, such as physical activity, air pollution, diet, and sex. The establishment of large consortia and population-based cohorts, in combination with new computational tools and mouse genetics platforms, can potentially overcome some of the limitations that have hindered human GxE interaction studies and reveal additional association signals for CVD-related traits. The identification of novel GxE interactions is likely to provide a better understanding of the pathogenesis and genetic liability of CVD, with significant implications for healthy lifestyles and therapeutic strategies.

The Meaning of "Cause" in Genetics

Causation has multiple distinct meanings in genetics. One reason for this is meaning slippage between two concepts of the gene: Mendelian and molecular. Another reason is that a variety of genetic methods address different kinds of causal relationships. Some genetic studies address causes of traits in individuals, which can only be assessed when single genes follow predictable inheritance patterns that reliably cause a trait. A second sense concerns the causes of trait differences within a population. Whereas some single genes can be said to cause population-level differences, most often these claims concern the effects of many genes. Polygenic traits can be understood using heritability estimates, which estimate the relative influences of genetic and environmental differences to trait differences within a population. Attempts to understand the molecular mechanisms underlying polygenic traits have been developed, although causal inference based on these results remains controversial. Genetic variation has also recently been leveraged as a randomizing factor to identify environmental causes of trait differences. This technique-Mendelian randomization-offers some solutions to traditional epidemiological challenges, although it is limited to the study of environments with known genetic influences.

PRICKLE1 × FOCAD Interaction Revealed by Genome-Wide vQTL Analysis of Human Facial Traits

The human face is a highly complex and variable structure resulting from the intricate coordination of numerous genetic and non-genetic factors. Hundreds of genomic loci impacting quantitative facial features have been identified. While these associations have been shown to influence morphology by altering the mean size and shape of facial measures, their effect on trait variance remains unclear. We conducted a genome-wide association analysis for the variance of 20 quantitative facial measurements in 2,447 European individuals and identified several suggestive variance quantitative trait loci (vQTLs). These vQTLs guided us to conduct an efficient search for gene-by-gene (G × G) interactions, which uncovered an interaction between PRICKLE1 and FOCAD affecting cranial base width. We replicated this G × G interaction signal at the locus level in an additional 5,128 Korean individuals. We used the hypomorphic Prickle1 Beetlejuice (Prickle1 Bj ) mouse line to directly test the function of Prickle1 on the cranial base and observed wider cranial bases in Prickle1 Bj/Bj . Importantly, we observed that the Prickle1 and Focadhesin proteins co-localize in murine cranial base chondrocytes, and this co-localization is abnormal in the Prickle1 Bj/Bj mutants. Taken together, our findings uncovered a novel G × G interaction effect in humans with strong support from both epidemiological and molecular studies. These results highlight the potential of studying measures of phenotypic variability in gene mapping studies of facial morphology.

How does genetic variation modify ND-CNV phenotypes?

Rare copy-number variants (CNVs) associated with neurodevelopmental disorders (NDDs), i.e., ND-CNVs, provide an insight into the neurobiology of NDDs and, potentially, a link between biology and clinical outcomes. However, ND-CNVs are characterised by incomplete penetrance resulting in heterogeneous carrier phenotypes, ranging from non-affected to multimorbid psychiatric, neurological, and physical phenotypes. Recent evidence indicates that other variants in the genome, or 'other hits', may partially explain the variable expressivity of ND-CNVs. These may be other rare variants or the aggregated effects of common variants that modify NDD risk. Here we discuss the recent findings, current questions, and future challenges relating to other hits research in the context of ND-CNVs and their potential for improved clinical diagnostics and therapeutics for ND-CNV carriers.

Sex-specific autosomal genetic effects across 26 human complex traits

Previous studies have shown that men and women have different genetic architectures across many traits. However, except waist-to-hip ratio (WHR) and waist circumference (WC), it remains unknown whether the genetic effects of a certain trait are weaker or stronger on men/women. With ~18 000 Taiwan Biobank subjects, we comprehensively investigate sexual heterogeneity in autosomal genetic effects, for traits regarding cardiovascular health, diabetes, kidney, liver, anthropometric profiles, blood, etc. 'Gene-by-sex interactions' (G $\times$ S) were detected in 18 out of 26 traits, each with an interaction P-value (${{P}}_{{INT}}$) less than $0.05/104={0.00048}$, where 104 is the number of tests conducted in this study. The most significant evidence of G $\times$ S was found in WHR (${{P}}_{{INT}}$ = 3.2 $\times{{10}}^{-{55}}$) and WC (${{P}}_{{INT}}$ = 2.3$\times{{10}}^{-{41}}$). As a novel G$\times$S investigation for other traits, we here find that the autosomal genetic effects are weaker on women than on men, for low-density lipoprotein cholesterol (LDL-C), uric acid (UA) and diabetes-related traits such as fasting glucose and glycated hemoglobin. For LDL-C and UA, the evidence of G$\times$S is especially notable in subjects aged less than 50 years, where estrogen can play a role in attenuating the autosomal genetic effects of these two traits. Men and women have systematically distinct environmental contexts caused by hormonal milieu and their specific society roles, which may trigger diverse gene expressions despite the same DNA materials. As many environmental exposures are difficult to collect and quantify, sex can serve as a good surrogate for these factors.

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.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Relevance of Higher-Order Epistasis in Drug Resistance

We studied five chemically distinct but related 1,3,5-triazine antifolates with regard to their effects on growth of a set of mutants in dihydrofolate reductase. The mutants comprise a combinatorially complete data set of all 16 possible combinations of four amino acid replacements associated with resistance to pyrimethamine in the malaria parasite Plasmodium falciparum. Pyrimethamine was a mainstay medication for malaria for many years, and it is still in use in intermittent treatment during pregnancy or as a partner drug in artemisinin combination therapy. Our goal was to investigate the extent to which the alleles yield similar adaptive topographies and patterns of epistasis across chemically related drugs. We find that the adaptive topographies are indeed similar with the same or closely related alleles being fixed in computer simulations of stepwise evolution. For all but one of the drugs the topography features at least one suboptimal fitness peak. Our data are consistent with earlier results indicating that third order and higher epistatic interactions appear to contribute only modestly to the overall adaptive topography, and they are largely conserved. In regard to drug development, our data suggest that higher-order interactions are likely to be of little value as an advisory tool in the choice of lead compounds.

Evaluating marginal genetic correlation of associated loci for complex diseases and traits between European and East Asian populations

Genome-wide association studies (GWASs) have successfully identified a large amount of single-nucleotide polymorphisms associated with many complex phenotypes in diverse populations. However, a comprehensive understanding of the genetic correlation of associated loci of phenotypes across populations remains lacking and the extent to which associations discovered in one population can be generalized to other populations or can be utilized for trans-ethnic genetic prediction is also unclear. By leveraging summary statistics, we proposed MAGIC to evaluate the trans-ethnic marginal genetic correlation (r_m) of per-allele effect sizes for associated SNPs (P < 5E-8) under the framework of measurement error models. We confirmed the methodological advantage of MAGIC over general approaches through simulations and demonstrated its utility by analyzing 34 GWAS summary statistics of phenotypes from the East Asian (N_max = 254,373) and European (N_max = 1,220,901) populations. Among these phenotypes, r_m was estimated to range from 0.584 (se = 0.140) for breast cancer to 0.949 (se = 0.035) for age of menarche, with an average of 0.835 (se = 0.045). We also uncovered that the trans-ethnic genetic prediction accuracy for phenotypes in the target population would substantially become low when using associated SNPs identified in non-target populations, indicating that associations discovered in the one population cannot be simply generalized to another population and that the accuracy of trans-ethnic phenotype prediction is generally dissatisfactory. Overall, our study provides in-depth insight into trans-ethnic genetic correlation and prediction for complex phenotypes across diverse populations.

Novel EDGE encoding method enhances ability to identify genetic interactions

Assumptions are made about the genetic model of single nucleotide polymorphisms (SNPs) when choosing a traditional genetic encoding: additive, dominant, and recessive. Furthermore, SNPs across the genome are unlikely to demonstrate identical genetic models. However, running SNP-SNP interaction analyses with every combination of encodings raises the multiple testing burden. Here, we present a novel and flexible encoding for genetic interactions, the elastic data-driven genetic encoding (EDGE), in which SNPs are assigned a heterozygous value based on the genetic model they demonstrate in a dataset prior to interaction testing. We assessed the power of EDGE to detect genetic interactions using 29 combinations of simulated genetic models and found it outperformed the traditional encoding methods across 10%, 30%, and 50% minor allele frequencies (MAFs). Further, EDGE maintained a low false-positive rate, while additive and dominant encodings demonstrated inflation. We evaluated EDGE and the traditional encodings with genetic data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes: age-related macular degeneration (AMD), age-related cataract, glaucoma, type 2 diabetes (T2D), and resistant hypertension. A multi-encoding genome-wide association study (GWAS) for each phenotype was performed using the traditional encodings, and the top results of the multi-encoding GWAS were considered for SNP-SNP interaction using the traditional encodings and EDGE. EDGE identified a novel SNP-SNP interaction for age-related cataract that no other method identified: rs7787286 (MAF: 0.041; intergenic region of chromosome 7)–rs4695885 (MAF: 0.34; intergenic region of chromosome 4) with a Bonferroni LRT p of 0.018. A SNP-SNP interaction was found in data from the UK Biobank within 25 kb of these SNPs using the recessive encoding: rs60374751 (MAF: 0.030) and rs6843594 (MAF: 0.34) (Bonferroni LRT p: 0.026). We recommend using EDGE to flexibly detect interactions between SNPs exhibiting diverse action.

Although traditional genetic encodings are widely implemented in genetics research, including in genome-wide association studies (GWAS) and epistasis, each method makes assumptions that may not reflect the underlying etiology. Here, we introduce a novel encoding method that estimates and assigns an individualized data-driven encoding for each single nucleotide polymorphism (SNP): the elastic data-driven genetic encoding (EDGE). With simulations, we demonstrate that this novel method is more accurate and robust than traditional encoding methods in estimating heterozygous genotype values, reducing the type I error, and detecting SNP-SNP interactions. We further applied the traditional encodings and EDGE to biomedical data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes, and EDGE identified a novel interaction for age-related cataract not detected by traditional methods, which replicated in data from the UK Biobank. EDGE provides an alternative approach to understanding and modeling diverse SNP models and is recommended for studying complex genetics in common human phenotypes.

Admixed Populations Improve Power for Variant Discovery and Portability in Genome-Wide Association Studies

Genome-wide association studies (GWAS) are primarily conducted in single-ancestry settings. The low transferability of results has limited our understanding of human genetic architecture across a range of complex traits. In contrast to homogeneous populations, admixed populations provide an opportunity to capture genetic architecture contributed from multiple source populations and thus improve statistical power. Here, we provide a mechanistic simulation framework to investigate the statistical power and transferability of GWAS under directional polygenic selection or varying divergence. We focus on a two-way admixed population and show that GWAS in admixed populations can be enriched for power in discovery by up to 2-fold compared to the ancestral populations under similar sample size. Moreover, higher accuracy of cross-population polygenic score estimates is also observed if variants and weights are trained in the admixed group rather than in the ancestral groups. Common variant associations are also more likely to replicate if first discovered in the admixed group and then transferred to an ancestral population, than the other way around (across 50 iterations with 1,000 causal SNPs, training on 10,000 individuals, testing on 1,000 in each population, p = 3.78e-6, 6.19e-101, ∼0 for F_ST = 0.2, 0.5, 0.8, respectively). While some of these F_ST values may appear extreme, we demonstrate that they are found across the entire phenome in the GWAS catalog. This framework demonstrates that investigation of admixed populations harbors significant advantages over GWAS in single-ancestry cohorts for uncovering the genetic architecture of traits and will improve downstream applications such as personalized medicine across diverse populations.

The Promise of AI in Detection, Diagnosis, and Epidemiology for Combating COVID-19: Beyond the Hype

COVID-19 has created enormous suffering, affecting lives, and causing deaths. The ease with which this type of coronavirus can spread has exposed weaknesses of many healthcare systems around the world. Since its emergence, many governments, research communities, commercial enterprises, and other institutions and stakeholders around the world have been fighting in various ways to curb the spread of the disease. Science and technology have helped in the implementation of policies of many governments that are directed toward mitigating the impacts of the pandemic and in diagnosing and providing care for the disease. Recent technological tools, artificial intelligence (AI) tools in particular, have also been explored to track the spread of the coronavirus, identify patients with high mortality risk and diagnose patients for the disease. In this paper, areas where AI techniques are being used in the detection, diagnosis and epidemiological predictions, forecasting and social control for combating COVID-19 are discussed, highlighting areas of successful applications and underscoring issues that need to be addressed to achieve significant progress in battling COVID-19 and future pandemics. Several AI systems have been developed for diagnosing COVID-19 using medical imaging modalities such as chest CT and X-ray images. These AI systems mainly differ in their choices of the algorithms for image segmentation, classification and disease diagnosis. Other AI-based systems have focused on predicting mortality rate, long-term patient hospitalization and patient outcomes for COVID-19. AI has huge potential in the battle against the COVID-19 pandemic but successful practical deployments of these AI-based tools have so far been limited due to challenges such as limited data accessibility, the need for external evaluation of AI models, the lack of awareness of AI experts of the regulatory landscape governing the deployment of AI tools in healthcare, the need for clinicians and other experts to work with AI experts in a multidisciplinary context and the need to address public concerns over data collection, privacy, and protection. Having a dedicated team with expertise in medical data collection, privacy, access and sharing, using federated learning whereby AI scientists hand over training algorithms to the healthcare institutions to train models locally, and taking full advantage of biomedical data stored in biobanks can alleviate some of problems posed by these challenges. Addressing these challenges will ultimately accelerate the translation of AI research into practical and useful solutions for combating pandemics.

Modeling genome-wide by environment interactions through omnigenic interactome networks

How genes interact with the environment to shape phenotypic variation and evolution is a fundamental question intriguing to biologists from various fields. Existing linear models built on single genes are inadequate to reveal the complexity of genotype-environment (G-E) interactions. Here, we develop a conceptual model for mechanistically dissecting G-E interplay by integrating previously disconnected theories and methods. Under this integration, evolutionary game theory, developmental modularity theory, and a variable selection method allow us to reconstruct environment-induced, maximally informative, sparse, and casual multilayer genetic networks. We design and conduct two mapping experiments by using a desert-adapted tree species to validate the biological application of the model proposed. The model identifies previously uncharacterized molecular mechanisms that mediate trees' response to saline stress. Our model provides a tool to comprehend the genetic architecture of trait variation and evolution and trace the information flow of each gene toward phenotypes within omnigenic networks.

Genetic impacts on DNA methylation: research findings and future perspectives

Multiple recent studies highlight that genetic variants can have strong impacts on a significant proportion of the human DNA methylome. Methylation quantitative trait loci, or meQTLs, allow for the exploration of biological mechanisms that underlie complex human phenotypes, with potential insights for human disease onset and progression. In this review, we summarize recent milestones in characterizing the human genetic basis of DNA methylation variation over the last decade, including heritability findings and genome-wide identification of meQTLs. We also discuss challenges in this field and future areas of research geared to generate insights into molecular processes underlying human complex traits.

The impact of late-career job loss and genetic risk on body mass index: Evidence from variance polygenic scores

Unemployment shocks from the COVID-19 pandemic have reignited concerns over the long-term effects of job loss on population health. Past research has highlighted the corrosive effects of unemployment on health and health behaviors. This study examines whether the effects of job loss on changes in body mass index (BMI) are moderated by genetic predisposition using data from the U.S. Health and Retirement Study (HRS). To improve detection of gene-by-environment (G × E) interplay, we interacted layoffs from business closures-a plausibly exogenous environmental exposure-with whole-genome polygenic scores (PGSs) that capture genetic contributions to both the population mean (mPGS) and variance (vPGS) of BMI. Results show evidence of genetic moderation using a vPGS (as opposed to an mPGS) and indicate genome-wide summary measures of phenotypic plasticity may further our understanding of how environmental stimuli modify the distribution of complex traits in a population.

Exploring a Region on Chromosome 8p23.1 Displaying Positive Selection Signals in Brazilian Admixed Populations: Additional Insights Into Predisposition to Obesity and Related Disorders

We recently reported a deviation of local ancestry on the chromosome (ch) 8p23.1, which led to positive selection signals in a Brazilian population sample. The deviation suggested that the genetic variability of candidate genes located on ch 8p23.1 may have been evolutionarily advantageous in the early stages of the admixture process. In the present work, we aim to extend the previous work by studying additional Brazilian admixed individuals and examining DNA sequencing data from the ch 8p23.1 candidate region. Thus, we inferred the local ancestry of 125 exomes from individuals born in five towns within the Southeast region of Brazil (São Paulo, Campinas, Barretos, and Ribeirão Preto located in the state of São Paulo and Belo Horizonte, the capital of the state of Minas Gerais), and compared to data from two public Brazilian reference genomic databases, BIPMed and ABraOM, and with information from the 1000 Genomes Project phase 3 and gnomAD databases. Our results revealed that ancestry is similar among individuals born in the five Brazilian towns assessed; however, an increased proportion of sub-Saharan African ancestry was observed in individuals from Belo Horizonte. In addition, individuals from the five towns considered, as well as those from the ABRAOM dataset, had the same overrepresentation of Native-American ancestry on the ch 8p23.1 locus that was previously reported for the BIPMed reference sample. Sequencing analysis of ch 8p23.1 revealed the presence of 442 non-synonymous variants, including frameshift, inframe deletion, start loss, stop gain, stop loss, and splicing site variants, which occurred in 24 genes. Among these genes, 13 were associated with obesity, type II diabetes, lipid levels, and waist circumference (PRAG1, MFHAS1, PPP1R3B, TNKS, MSRA, PRSS55, RP1L1, PINX1, MTMR9, FAM167A, BLK, GATA4, and CTSB). These results strengthen the hypothesis that a set of variants located on ch 8p23.1 that result from positive selection during early admixture events may influence obesity-related disease predisposition in admixed individuals of the Brazilian population. Furthermore, we present evidence that the exploration of local ancestry deviation in admixed individuals may provide information with the potential to be translated into health care improvement.

Assessing the co-variability of DNA methylation across peripheral cells and tissues: Implications for the interpretation of findings in epigenetic epidemiology

Most epigenome-wide association studies (EWAS) quantify DNA methylation (DNAm) in peripheral tissues such as whole blood to identify positions in the genome where variation is statistically associated with a trait or exposure. As whole blood comprises a mix of cell types, it is unclear whether trait-associated DNAm variation is specific to an individual cellular population. We collected three peripheral tissues (whole blood, buccal epithelial and nasal epithelial cells) from thirty individuals. Whole blood samples were subsequently processed using fluorescence-activated cell sorting (FACS) to purify five constituent cell-types (monocytes, granulocytes, CD4⁺ T cells, CD8⁺ T cells, and B cells). DNAm was profiled in all eight sample-types from each individual using the Illumina EPIC array. We identified significant differences in both the level and variability of DNAm between different sample types, and DNAm data-derived estimates of age and smoking were found to differ dramatically across sample types from the same individual. We found that for the majority of loci variation in DNAm in individual blood cell types was only weakly predictive of variance in DNAm measured in whole blood, although the proportion of variance explained was greater than that explained by either buccal or nasal epithelial samples. Covariation across sample types was much higher for DNAm sites influenced by genetic factors. Overall, we observe that DNAm variation in whole blood is additively influenced by a combination of the major blood cell types. For a subset of sites, however, variable DNAm detected in whole blood can be attributed to variation in a single blood cell type providing potential mechanistic insight about EWAS findings. Our results suggest that associations between whole blood DNAm and traits or exposures reflect differences in multiple cell types and our data will facilitate the interpretation of findings in epigenetic epidemiology.

As epigenetic variation is cell-type specific, an ongoing challenge in epigenetic epidemiology is how to interpret studies performed using bulk tissue (for example, whole blood) which comprises a mix of different cell types. In this study, we identified major differences in DNA methylation (DNAm) across multiple peripheral tissues and different blood cell types, with each sample type being characterized by a unique signature across multiple genomic loci. We demonstrate how these differences influence commonly used prediction scores derived from DNAm data for age and tobacco smoking, with estimates for the same individual being highly variable across tissues and cell types. Our results enabled us to assess the extent to which variable DNAm in each individual blood cell type relates to variation measured in whole blood. We found that although individual blood cell types predict more of the variation in DNAm in whole blood compared to buccal and nasal epithelial cells, the actual proportion of variance explained is relatively small, except for at sites where DNAm is under genetic control. Our data indicate that for most sites variation in multiple blood cell types additively combines to drive variation in DNAm measured in whole blood. Of note, for a subset of sites, variation in DNAm detected in whole blood can be attributed to a specific blood cell type, potentially facilitating the interpretation of EWAS findings.

Leveraging phenotypic variability to identify genetic interactions in human phenotypes

Although thousands of loci have been associated with human phenotypes, the role of gene-environment (GxE) interactions in determining individual risk of human diseases remains unclear. This is partly because of the severe erosion of statistical power resulting from the massive number of statistical tests required to detect such interactions. Here, we focus on improving the power of GxE tests by developing a statistical framework for assessing quantitative trait loci (QTLs) associated with the trait means and/or trait variances. When applying this framework to body mass index (BMI), we find that GxE discovery and replication rates are significantly higher when prioritizing genetic variants associated with the variance of the phenotype (vQTLs) compared to when assessing all genetic variants. Moreover, we find that vQTLs are enriched for associations with other non-BMI phenotypes having strong environmental influences, such as diabetes or ulcerative colitis. We show that GxE effects first identified in quantitative traits such as BMI can be used for GxE discovery in disease phenotypes such as diabetes. A clear conclusion is that strong GxE interactions mediate the genetic contribution to body weight and diabetes risk.

Estimating the Prevalence and Genetic Risk Mechanisms of ARFID in a Large Autism Cohort

This study is the first genetically-informed investigation of avoidant/restrictive food intake disorder (ARFID), an eating disorder that profoundly impacts quality of life for those affected. ARFID is highly comorbid with autism, and we provide the first estimate of its prevalence in a large and phenotypically diverse autism cohort (a subsample of the SPARK study, N = 5,157 probands). This estimate, 21% (at a balanced accuracy 80%), is at the upper end of previous estimates from studies based on clinical samples, suggesting under-diagnosis and potentially lack of awareness among caretakers and clinicians. Although some studies suggest a decrease of disordered eating symptoms by age 6, our estimates indicate that up to 17% (at a balanced accuracy 87%) of parents of autistic children are also at heightened risk for ARFID, suggesting a lifelong risk for disordered eating. We were also able to provide the first estimates of narrow-sense heritability (h²) for ARFID risk, at 0.45. Genome-wide association revealed a single hit near ZSWIM6, a gene previously implicated in neurodevelopmental conditions. While, the current sample was not well-powered for GWAS, effect size and heritability estimates allowed us to project the sample sizes necessary to more robustly discover ARFID-linked loci via common variants. Further genetic analysis using polygenic risk scores (PRS) affirmed genetic links to autism as well as neuroticism and metabolic syndrome.

How wide is the application of genetic big data in biomedicine

In the era of big data, massive genetic data, as a new industry, has quickly swept almost all industries, especially the pharmaceutical industry. As countries around the world start to build their own gene banks, scientists study the data to explore the origins and migration of humans. Moreover, big data encourage the development of cancer therapy and bring good news to cancer patients. Big datum has been involved in the study of many diseases, and it has been found that analyzing diseases at the gene level can lead to more beneficial treatment options than ordinary treatments. This review will introduce the development of extensive data in medical research from the perspective of big data and tumor, neurological and psychiatric diseases, cardiovascular diseases, other applications and the development direction of big data in medicine.

The Parkinson's Disease Genome-Wide Association Study Locus Browser

BACKGROUND

Parkinson's disease (PD) is a neurodegenerative disease with an often complex component identifiable by genome-wide association studies. The most recent large-scale PD genome-wide association studies have identified more than 90 independent risk variants for PD risk and progression across more than 80 genomic regions. One major challenge in current genomics is the identification of the causal gene(s) and variant(s) at each genome-wide association study locus. The objective of the current study was to create a tool that would display data for relevant PD risk loci and provide guidance with the prioritization of causal genes and potential mechanisms at each locus.

METHODS

We included all significant genome-wide signals from multiple recent PD genome-wide association studies including themost recent PD risk genome-wide association study, age-at-onset genome-wide association study, progression genome-wide association study, and Asian population PD risk genome-wide association study. We gathered data for all genes 1 Mb up and downstream of each variant to allow users to assess which gene(s) are most associated with the variant of interest based on a set of self-ranked criteria. Multiple databases were queried for each gene to collect additional causal data.

RESULTS

We created a PD genome-wide association study browser tool (https://pdgenetics.shinyapps.io/GWASBrowser/) to assist the PD research community with the prioritization of genes for follow-up functional studies to identify potential therapeutic targets.

CONCLUSIONS

Our PD genome-wide association study browser tool provides users with a useful method of identifying potential causal genes at all known PD risk loci from large-scale PD genome-wide association studies. We plan to update this tool with new relevant data as sample sizes increase and new PD risk loci are discovered. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Large-Scale Analyses Provide No Evidence for Gene-Gene Interactions Influencing Type 2 Diabetes Risk

A growing number of genetic loci have been shown to influence individual predisposition to type 2 diabetes (T2D). Despite longstanding interest in understanding whether nonlinear interactions between these risk variants additionally influence T2D risk, the ability to detect significant gene-gene interaction (GGI) effects has been limited to date. To increase power to detect GGI effects, we combined recent advances in the fine-mapping of causal T2D risk variants with the increased sample size available within UK Biobank (375,736 unrelated European participants, including 16,430 with T2D). In addition to conventional single variant-based analysis, we used a complementary polygenic score-based approach, which included partitioned T2D risk scores that capture biological processes relevant to T2D pathophysiology. Nevertheless, we found no evidence in support of GGI effects influencing T2D risk. The current study was powered to detect interactions between common variants with odds ratios >1.2, so these findings place limits on the contribution of GGIs to the overall heritability of T2D.

Genetics of Obesity in East Asians

Obesity has become a public health problem worldwide. Compared with Europe, people in Asia tend to suffer from type 2 diabetes with a lower body mass index (BMI). Genome-wide association studies (GWASs) have identified over 750 loci associated with obesity. Although the majority of GWAS results were conducted in individuals of European ancestry, a recent GWAS in individuals of Asian ancestry has made a significant contribution to the identification of obesity susceptibility loci. Indeed, owing to the multifactorial character of obesity with a strong environmental component, the revealed loci may have distinct contributions in different ancestral genetic backgrounds and in different environments as presented through diet and exercise among other factors. Uncovering novel, yet unrevealed genes in non-European ancestries may further contribute to explaining the missing heritability for BMI. In this review, we aimed to summarize recent advances in obesity genetics in individuals of Asian ancestry. We therefore compared proposed mechanisms underlying susceptibility loci for obesity associated with individuals of European and Asian ancestries and discussed whether known genetic variants might explain ethnic differences in obesity risk. We further acknowledged that GWAS implemented in individuals of Asian ancestries have not only validated the potential role of previously specified obesity susceptibility loci but also exposed novel ones, which have been missed in the initial genetic studies in individuals of European ancestries. Thus, multi-ethnic studies have a great potential not only to contribute to a better understanding of the complex etiology of human obesity but also potentially of ethnic differences in the prevalence of obesity, which may ultimately pave new avenues in more targeted and personalized obesity treatments.

The interplay between host genetics and the gut microbiome reveals common and distinct microbiome features for complex human diseases

BACKGROUND

Interest in the interplay between host genetics and the gut microbiome in complex human diseases is increasing, with prior evidence mainly being derived from animal models. In addition, the shared and distinct microbiome features among complex human diseases remain largely unclear.

RESULTS

This analysis was based on a Chinese population with 1475 participants. We estimated the SNP-based heritability, which suggested that Desulfovibrionaceae and Odoribacter had significant heritability estimates (0.456 and 0.476, respectively). We performed a microbiome genome-wide association study to identify host genetic variants associated with the gut microbiome. We then conducted bidirectional Mendelian randomization analyses to examine the potential causal associations between the gut microbiome and complex human diseases. We found that Saccharibacteria could potentially decrease the concentration of serum creatinine and increase the estimated glomerular filtration rate. On the other hand, atrial fibrillation, chronic kidney disease and prostate cancer, as predicted by host genetics, had potential causal effects on the abundance of some specific gut microbiota. For example, atrial fibrillation increased the abundance of Burkholderiales and Alcaligenaceae and decreased the abundance of Lachnobacterium, Bacteroides coprophilus, Barnesiellaceae, an undefined genus in the family Veillonellaceae and Mitsuokella. Further disease-microbiome feature analysis suggested that systemic lupus erythematosus and chronic myeloid leukaemia shared common gut microbiome features.

CONCLUSIONS

These results suggest that different complex human diseases share common and distinct gut microbiome features, which may help reshape our understanding of disease aetiology in humans. Video Abstract.

Inferring Gene-by-Environment Interactions with a Bayesian Whole-Genome Regression Model

The contribution of gene-by-environment (GxE) interactions for many human traits and diseases is poorly characterized. We propose a Bayesian whole-genome regression model for joint modeling of main genetic effects and GxE interactions in large-scale datasets, such as the UK Biobank, where many environmental variables have been measured. The method is called LEMMA (Linear Environment Mixed Model Analysis) and estimates a linear combination of environmental variables, called an environmental score (ES), that interacts with genetic markers throughout the genome. The ES provides a readily interpretable way to examine the combined effect of many environmental variables. The ES can be used both to estimate the proportion of phenotypic variance attributable to GxE effects and to test for GxE effects at genetic variants across the genome. GxE effects can induce heteroskedasticity in quantitative traits, and LEMMA accounts for this by using robust standard error estimates when testing for GxE effects. When applied to body mass index, systolic blood pressure, diastolic blood pressure, and pulse pressure in the UK Biobank, we estimate that 9.3%, 3.9%, 1.6%, and 12.5%, respectively, of phenotypic variance is explained by GxE interactions and that low-frequency variants explain most of this variance. We also identify three loci that interact with the estimated environmental scores (−log10p>7.3).

Interactions between Polygenic Scores and Environments: Methodological and Conceptual Challenges

Interest in the study of gene-environment interaction has recently grown due to the sudden availability of molecular genetic data-in particular, polygenic scores-in many long-running longitudinal studies. Identifying and estimating statistical interactions comes with several analytic and inferential challenges; these challenges are heightened when used to integrate observational genomic and social science data. We articulate some of these key challenges, provide new perspectives on the study of gene-environment interactions, and end by offering some practical guidance for conducting research in this area. Given the sudden availability of well-powered polygenic scores, we anticipate a substantial increase in research testing for interaction between such scores and environments. The issues we discuss, if not properly addressed, may impact the enduring scientific value of gene-environment interaction studies.

Schizophrenia: Developmental Variability Interacts with Risk Factors to Cause the Disorder: Nonspecific Variability-Enhancing Factors Combine with Specific Risk Factors to Cause Schizophrenia

A new etiological model is proposed for schizophrenia that combines variability-enhancing nonspecific factors acting during development with more specific risk factors. This model is better suited than the current etiological models of schizophrenia, based on the risk factors paradigm, for predicting and/or explaining several important findings about schizophrenia: high co-morbidity rates, low specificity of many risk factors, and persistence in the population of the associated genetic polymorphisms. Compared with similar models, e.g., de-canalization, common psychopathology factor, sexual-selection, or differential sensitivity to the environment, this proposal is more general and integrative. Recently developed research methods have proven the existence of genetic and environmental factors that enhance developmental variability. Applying such methods to newly collected or already available data can allow for testing the hypotheses upon which this model is built. If validated, this model may change the understanding of the etiology of schizophrenia, the research models, and preventionbrk paradigms.

Genetic background and diet affect brown adipose gene coexpression networks associated with metabolic phenotypes

Adipose is a dynamic endocrine organ that is critical for regulating metabolism and is highly responsive to nutritional environment. Brown adipose tissue is an exciting potential therapeutic target; however, there are no systematic studies of gene-by-environment interactions affecting function of this organ. We leveraged a weighted gene coexpression network analysis to identify transcriptional networks in brown adipose tissue from LG/J and SM/J inbred mice fed high- or low-fat diets and correlate these networks with metabolic phenotypes. We identified eight primary gene network modules associated with variation in obesity and diabetes-related traits. Four modules were enriched for metabolically relevant processes such as immune and cytokine response, cell division, peroxisome functions, and organic molecule metabolic processes. The relative expression of genes in these modules is highly dependent on both genetic background and dietary environment. Genes in the immune/cytokine response and cell division modules are particularly highly expressed in high fat-fed SM/J mice, which show unique brown adipose-dependent remission of diabetes. The interconnectivity of genes in these modules is also heavily dependent on diet and strain, with most genes showing both higher expression and coexpression under the same context. We highlight several genes of interest, Col28a1, Cyp26b1, Bmp8b, and Ngef, that have distinct expression patterns among strain-by-diet contexts and fall under metabolic quantitative trait loci previously mapped in an F16 generation of an advanced intercross between LG/J and SM/J. Each of these genes have some connection to obesity and diabetes-related traits, but have not been studied in brown adipose tissue. Our results provide important insights into the relationship between brown adipose and systemic metabolism by being the first gene-by-environment study of brown adipose transcriptional networks.