Assessing Gene-Environment Interactions for Common and Rare Variants with Binary Traits Using Gene-Trait Similarity Regression

Accounting for gene–environment (G×E) interactions in complex trait association studies can facilitate our understanding of genetic heterogeneity under different environmental exposures, improve the ability to discover susceptible genes that exhibit little marginal effect, provide insight into the biological mechanisms of complex diseases, help to identify high-risk subgroups in the population, and uncover hidden heritability. However, significant G×E interactions can be difficult to find. The sample sizes required for sufficient power to detect association are much larger than those needed for genetic main effects, and interactions are sensitive to misspecification of the main-effects model. These issues are exacerbated when working with binary phenotypes and rare variants, which bear less information on association. In this work, we present a similarity-based regression method for evaluating G×E interactions for rare variants with binary traits. The proposed model aggregates the genetic and G×E information across markers, using genetic similarity, thus increasing the ability to detect G×E signals. The model has a random effects interpretation, which leads to robustness against main-effect misspecifications when evaluating G×E interactions. We construct score tests to examine G×E interactions and a computationally efficient EM algorithm to estimate the nuisance variance components. Using simulations and data applications, we show that the proposed method is a flexible and powerful tool to study the G×E effect in common or rare variant studies with binary traits.     

The adaptive potential of a plant pathogenic fungus, <Emphasis Type="Italic">Rhizoctonia solani</Emphasis> AG-3, under heat and fungicide stress

The ability to improve fitness via adaptive evolution may be affected by environmental change. We tested this hypothesis in an in vitro experiment with the plant pathogen Rhizoctonia solani Anastomosis Group 3 (AG-3), assessing genetic and environmental variances under two temperatures (optimal and higher than optimal) and three fungicide concentrations (no fungicide, low and high concentration of a copper-based fungicide). We measured the mean daily growth rate, the coefficient of variation for genotypic (I _G) and environmental variance (I _E) in growth, and broad-sense heritability in growth. Both higher temperature and increased fungicide concentration caused a decline in growth, confirming their potential as stressors for the pathogen. All types of standardized variances in growth—I _G, phenotypic variance, and I _E as a trend—increased with elevated stress. However, heritability was not significantly higher under enhanced stress because the increase in I _G was counterbalanced by somewhat increased I _E. The results illustrate that predictions for adaptation under environmental stress may depend on the type of short-term evolvability measure. Because mycelial growth is linked to fitness, I _G reflects short-term evolvability better than heritability, and it indicates that the evolutionary potential of R. solani is positively affected by stress.     

Dominance and environmental variances: their effect on heritabilities estimated from twin data

A method for partitioning genetic variance estimated from twin data into additive and dominance variances was presented using Falconer's variance component model. The effects of dominance and environmental variances on a number of heritability estimates were also reviewed. A heritability estimate, based on the analysis of variance and the genetic variance estimates presented by HASEMAN and ELSTON and CHRISTIAN et al. which utilizes all available information from twin data, was proposed and discussed. This estimate seems to be the least affected by fluctuations in the magnitudes of dominance and environmental variances.     

Genetic and Environmental Contributions to BMI in Adolescent and Young Adult Women

The objective of this study was to determine the genetic and environmental contributions to variation in BMI over time in European‐American (EA) and African‐American (AA) adolescent and young adult women. Self‐reported BMI (kg/m²) data from 2,816 EA (1,306 twin pairs, 56.5% monozygotic (MZ)) and 404 AA (178 twin pairs, 42.7% MZ) women at baseline (T1; median age 15 years) and 3,225 EA (1,511 twin pairs, 55.3% MZ) and 539 AA (252 pairs, 43.3% MZ) women at follow‐up (T2; median age 22 years) from a Midwestern US, population‐based twin registry were used to construct biometrical genetic models. For EA women, the majority of the variance in BMI was attributable to additive genetic effects at both time points (82% for each), with the remaining variance attributable to nonshared environment. Genetic and nonshared environment correlations between adolescent and young adult BMI were 0.87 and 0.23, respectively. Among AA women, nonadditive genetic effects comprised 68% of the variance at T1 and 73% at T2, and were highly correlated (r_D = 0.94). The proportions of variance attributable to nonshared environment at T1 (29%) and T2 (25%) were more modestly correlated (r_E = 0.31). The remaining variance in AA women could be attributed to additive genetic effects. Additive vs. nonadditive genetic effects contribute differentially to BMI in AA vs. EA adolescent and young adult women. Additional research is needed to better characterize the environmental and genetic factors related to BMI in persons of different races to aid understanding of the complex determinants of body weight in individuals.     

Genetic architecture of lipid traits changes over time and differs by race: Princeton Lipid Follow-up Study

Dyslipidemia is a major risk factor for CVD. Previous studies on lipid heritability have largely focused on white populations assessed after the obesity epidemic. Given secular trends and racial differences in lipid levels, this study explored whether lipid heritability is consistent across time and between races. African American and white nuclear families had fasting lipids measured in the 1970s and 22–30 years later. Heritability was estimated, and bivariate analyses between visits were conducted by race using variance components analysis. A total of 1,454 individuals (age 14.1/40.6 for offspring/parents at baseline; 39.6/66.5 at follow-up) in 373 families (286 white, 87 African American) were included. Lipid trait heritabilities were typically stronger during the 1970s than the 2000s. At baseline, additive genetic variation for LDL was significantly lower in African Americans than whites (P = 0.015). Shared genetic contribution to lipid variability over time was significant in both whites (all P < 0.0001) and African Americans (P ≤ 0.05 for total, LDL, and HDL cholesterol). African American families demonstrated shared environmental contributions to lipid variation over time (all P ≤ 0.05). Lower heritability, lower LDL genetic variance, and durable environmental effects across the obesity epidemic in African American families suggest race-specific approaches are needed to clarify the genetic etiology of lipids.     

Glycosylation of Immunoglobulin G: Role of Genetic and Epigenetic Influences

Objective

To determine the extent to which genetic and epigenetic factors contribute to variations in glycosylation of immunoglobulin G (IgG) in humans.

Methods

76  N-glycan traits in circulating IgG were analyzed by UPLC in 220 monozygotic and 310 dizygotic twin pairs from TwinsUK. A classical twin study design was used to derive the additive genetic, common and unique environmental components defining the variance in these traits. Epigenome-wide association analysis was performed using the Illumina 27k chip.

Results

51 of the 76 glycan traits studied have an additive genetic component (heritability, h ²)≥  0.5. In contrast, 12 glycan traits had a low genetic contribution (h²<0.35). We then tested for association between methylation levels and glycan levels (P<2 x10^-6). Among glycan traits with low heritability probe cg08392591 maps to a CpG island 5’ from the ANKRD11 gene, a p53 activator on chromosome 16. Probe cg26991199 maps to the SRSF10 gene involved in regulation of RNA splicing and particularly in regulation of splicing of mRNA precursors upon heat shock. Among those with high heritability we found cg13782134 (mapping to the NRN1L gene) and cg16029957 mapping near the QPCT gene to be array-wide significant. The proportion of array-wide epigenetic associations was significantly larger (P<0.005) among glycans with low heritability (42%) than in those with high heritability (6.2%).

Conclusions

Glycome analyses might provide a useful integration of genetic and non-genetic factors to further our understanding of the role of glycosylation in both normal physiology and disease.     

Integrated analysis of direct and proxy genome wide association studies highlights polygenicity of Alzheimer’s disease outside of the APOE region

Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer’s disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer’s Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7–11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.     

Genetic Epidemiology of BMI and Body Mass Change From Adolescence to Young Adulthood

The complex interplay between genes and environment affecting body mass gain over lifecycle periods of risk is not well understood. We use longitudinal sibling cohort data to examine the role of shared household environment, additive genetic, and shared genetic effects on BMI and BMI change. In the National Longitudinal Study of Adolescent Health, siblings and twin pairs sharing households for ≥10 years as adolescents (N = 5,524; mean = 16.5 ± 1.7 years) were followed into young adulthood (N = 4,368; mean = 22.4 ± 1.8 years). Using a variance component approach, we quantified genetic and household effects on BMI in siblings and nonsiblings sharing household environments over time. Adjusting for race, age, sex, and age‐by‐sex interaction, we detected a heritability of 0.43 ± 0.05 for BMI change. Significant household effects were noted during the young adulthood period only (0.11 ± 0.06). We find evidence for shared genetic effects between BMI and BMI change during adolescence (genetic correlation (ρ_G) = 0.61 ± 0.03) and young adulthood (ρ_G = 0.23 ± 0.06). Our findings support a complex etiology of BMI and BMI change.     

Genetic Dissection of the Drosophila melanogaster Female Head Transcriptome Reveals Widespread Allelic Heterogeneity

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.     

Dominant Genetic Variation and Missing Heritability for Human Complex Traits: Insights from Twin versus Genome-wide Common SNP Models

In order to further illuminate the potential role of dominant genetic variation in the “missing heritability” debate, we investigated the additive (narrow-sense heritability, h²) and dominant (δ²) genetic variance for 18 human complex traits. Within the same study base (10,682 Swedish twins), we calculated and compared the estimates from classic twin-based structural equation model with SNP-based genomic-relatedness-matrix restricted maximum likelihood [GREML(d)] method. Contributions of δ² were evident for 14 traits in twin models (average δ²_twin = 0.25, range 0.14–0.49), two of which also displayed significant δ² in the GREMLd analyses (triglycerides δ²_SNP = 0.28 and waist circumference δ²_SNP = 0.19). On average, the proportion of h²_SNP/h²_twin was 70% for ADE-fitted traits (for which the best-fitting model included additive and dominant genetic and unique environmental components) and 31% for AE-fitted traits (for which the best-fitting model included additive genetic and unique environmental components). Independent evidence for contribution from shared environment, also in ADE-fitted traits, was obtained from self-reported within-pair contact frequency and age at separation. We conclude that despite the fact that additive genetics appear to constitute the bulk of genetic influences for most complex traits, dominant genetic variation might often be masked by shared environment in twin and family studies and might therefore have a more prominent role than what family-based estimates often suggest. The risk of erroneously attributing all inherited genetic influences (additive and dominant) to the h² in too-small twin studies might also lead to exaggerated “missing heritability” (the proportion of h² that remains unexplained by SNPs).     

Mapping the Genetic Variation of Regional Brain Volumes as Explained by All Common SNPs from the ADNI Study

Typically twin studies are used to investigate the aggregate effects of genetic and environmental influences on brain phenotypic measures. Although some phenotypic measures are highly heritable in twin studies, SNPs (single nucleotide polymorphisms) identified by genome-wide association studies (GWAS) account for only a small fraction of the heritability of these measures. We mapped the genetic variation (the proportion of phenotypic variance explained by variation among SNPs) of volumes of pre-defined regions across the whole brain, as explained by 512,905 SNPs genotyped on 747 adult participants from the Alzheimer''s Disease Neuroimaging Initiative (ADNI). We found that 85% of the variance of intracranial volume (ICV) (p = 0.04) was explained by considering all SNPs simultaneously, and after adjusting for ICV, total grey matter (GM) and white matter (WM) volumes had genetic variation estimates near zero (p = 0.5). We found varying estimates of genetic variation across 93 non-overlapping regions, with asymmetry in estimates between the left and right cerebral hemispheres. Several regions reported in previous studies to be related to Alzheimer''s disease progression were estimated to have a large proportion of volumetric variance explained by the SNPs.     

Genetic Influences Are Virtually Absent for Trust

Over the past decades, numerous twin studies have revealed moderate to high heritability estimates for individual differences in a wide range of human traits, including cognitive ability, psychiatric disorders, and personality traits. Even factors that are generally believed to be environmental in nature have been shown to be under genetic control, albeit modest. Is such heritability also present in social traits that are conceptualized as causes and consequences of social interactions or in other ways strongly shaped by behavior of other people? Here we examine a population-based sample of 1,012 twins and relatives. We show that the genetic influence on generalized trust in other people (trust-in-others: h² = 5%, ns), and beliefs regarding other people’s trust in the self (trust-in-self: h² = 13%, ns), is virtually absent. As test-retest reliability for both scales were found to be moderate or high (r = .76 and r = .53, respectively) in an independent sample, we conclude that all variance in trust is likely to be accounted for by non-shared environmental influences. We show that, relative to cognitive abilities, psychiatric disorders, and classic personality variables, genetic influences are smaller for trust, and propose that experiences with or observations of the behavior of other people shape trust more strongly than other traits.     

Polygenic and socioeconomic risk for high body mass index: 69 years of follow-up across life

Genetic influences on body mass index (BMI) appear to markedly differ across life, yet existing research is equivocal and limited by a paucity of life course data. We thus used a birth cohort study to investigate differences in association and explained variance in polygenic risk for high BMI across infancy to old age (2–69 years). A secondary aim was to investigate how the association between BMI and a key purported environmental determinant (childhood socioeconomic position) differed across life, and whether this operated independently and/or multiplicatively of genetic influences. Data were from up to 2677 participants in the MRC National Survey of Health and Development, with measured BMI at 12 timepoints from 2–69 years. We used multiple polygenic indices from GWAS of adult and childhood BMI, and investigated their associations with BMI at each age. For polygenic liability to higher adult BMI, the trajectories of effect size (β) and explained variance (R²) diverged: explained variance peaked in early adulthood and plateaued thereafter, while absolute effect sizes increased throughout adulthood. For polygenic liability to higher childhood BMI, explained variance was largest in adolescence and early adulthood; effect sizes were marginally smaller in absolute terms from adolescence to adulthood. All polygenic indices were related to higher variation in BMI; quantile regression analyses showed that effect sizes were sizably larger at the upper end of the BMI distribution. Socioeconomic and polygenic risk for higher BMI across life appear to operate additively; we found little evidence of interaction. Our findings highlight the likely independent influences of polygenic and socioeconomic factors on BMI across life. Despite sizable associations, the BMI variance explained by each plateaued or declined across adulthood while BMI variance itself increased. This is suggestive of the increasing importance of chance (‘non-shared’) environmental influences on BMI across life.     

Genetics of Callous-Unemotional Behavior in Children

Callous-unemotional behavior (CU) is currently under consideration as a subtyping index for conduct disorder diagnosis. Twin studies routinely estimate the heritability of CU as greater than 50%. It is now possible to estimate genetic influence using DNA alone from samples of unrelated individuals, not relying on the assumptions of the twin method. Here we use this new DNA method (implemented in a software package called Genome-wide Complex Trait Analysis, GCTA) for the first time to estimate genetic influence on CU. We also report the first genome-wide association (GWA) study of CU as a quantitative trait. We compare these DNA results to those from twin analyses using the same measure and the same community sample of 2,930 children rated by their teachers at ages 7, 9 and 12. GCTA estimates of heritability were near zero, even though twin analysis of CU in this sample confirmed the high heritability of CU reported in the literature, and even though GCTA estimates of heritability were substantial for cognitive and anthropological traits in this sample. No significant associations were found in GWA analysis, which, like GCTA, only detects additive effects of common DNA variants. The phrase ‘missing heritability’ was coined to refer to the gap between variance associated with DNA variants identified in GWA studies versus twin study heritability. However, GCTA heritability, not twin study heritability, is the ceiling for GWA studies because both GCTA and GWA are limited to the overall additive effects of common DNA variants, whereas twin studies are not. This GCTA ceiling is very low for CU in our study, despite its high twin study heritability estimate. The gap between GCTA and twin study heritabilities will make it challenging to identify genes responsible for the heritability of CU.     

Epidemiology and Genetic Epidemiology of the Liver Function Test Proteins

Background

The liver function test (LFT) is among the most commonly used clinical investigations to assess hepatic function, severity of liver diseases and the effect of therapies, as well as to detect drug-induced liver injury (DILI).

Aims

To determine the relative contribution of genetic and environmental factors as well as test and quantify the effects of sex, age, BMI and alcohol consumption to variation in liver function test proteins - including alanine amino transaminase (ALT), Albumin, gamma glutamyl transpeptidase (GGT), total bilirubin, total protein, total globulin, aspartate transaminase (AST), and alkaline phosphotase (ALP) - using the classical twin model.

Methods

Blood samples were collected from a total of 5380 twin pairs from the TwinsUK registry. We measured the expression levels of major proteins associated with the LFT, calculated BMI from measured weight and height and questionnaires were completed for alcohol consumption by the twins. The relative contribution of genetic and environmental factors to variation in the LFT proteins was assessed and quantified using a variance components model fitting approach.

Results

Our results show that (1) variation in all the LFTs has a significant heritable basis (h² ranging from 20% to 77%); (2) other than GGT, the LFTs are all affected to some extent by common environmental factors (c² ranging from 24% to 54%); and (3) a small but significant proportion of the variation in the LFTs was due to confounding effects of age, sex, BMI, and alcohol use.

Conclusions

Variation in the LFT proteins is under significant genetic and common environmental control although sex, alcohol use, age and BMI also contribute significantly to inter-individual variation in the LFT proteins. Understanding the underlying genetic contribution of liver function tests may help the interpretation of their results and explain wide variation among individuals.     

Genetic and Environmental Influences on Body‐Fat Measures among African‐American Twins

Objective: To investigate the genetic and environmental influences on body‐fat measures including waist circumference (WC), waist‐to‐hip ratio (WHR), and body mass index (BMI) among African‐American men and women. Research Methods and Procedures: Measurements were taken as part of the Carolina African American Twin Study of Aging. This sample currently comprises 146 same‐sex African‐American twins with an average age of 50 years (range, 22 to 88 years). This analysis included 26 monozygotic and 29 dizygotic men and 45 monozygotic and 46 dizygotic women. Maximum likelihood quantitative genetic analysis was used. Results: In men, additive genetic effects accounted for 77% of the variance in WC, 59% in WHR, and 89% in BMI. In women, additive genetic effects accounted for 76% of the variance in WC, 56% in WHR, and 73% in BMI. The remaining variance in both men and women was attributed to unique environmental effects (WC, 21%; WHR, 36%; BMI, 11% in men and WC, 22%; WHR, 38%; BMI, 27% in women) and age (WC, 2%; WHR, 5% in men and WC, 2%; WHR, 6% in women). When BMI was controlled in the analysis of WC and WHR, it accounted for a portion of the genetic and environmental variance in WHR and over one‐half of the genetic and environmental variance in WC. Discussion: There are both genetic and environmental influences on WC, WHR, and BMI, and independent of BMI, there are genetic and environmental effects on WC and WHR among both genders. The results from this African‐American twin sample are similar to findings among white twin samples.     

The Rate of Nonallelic Homologous Recombination in Males Is Highly Variable,Correlated between Monozygotic Twins and Independent of Age

Nonallelic homologous recombination (NAHR) between highly similar duplicated sequences generates chromosomal deletions, duplications and inversions, which can cause diverse genetic disorders. Little is known about interindividual variation in NAHR rates and the factors that influence this. We estimated the rate of deletion at the CMT1A-REP NAHR hotspot in sperm DNA from 34 male donors, including 16 monozygotic (MZ) co-twins (8 twin pairs) aged 24 to 67 years old. The average NAHR rate was 3.5×10⁻⁵ with a seven-fold variation across individuals. Despite good statistical power to detect even a subtle correlation, we observed no relationship between age of unrelated individuals and the rate of NAHR in their sperm, likely reflecting the meiotic-specific origin of these events. We then estimated the heritability of deletion rate by calculating the intraclass correlation (ICC) within MZ co-twins, revealing a significant correlation between MZ co-twins (ICC = 0.784, p = 0.0039), with MZ co-twins being significantly more correlated than unrelated pairs. We showed that this heritability cannot be explained by variation in PRDM9, a known regulator of NAHR, or variation within the NAHR hotspot itself. We also did not detect any correlation between Body Mass Index (BMI), smoking status or alcohol intake and rate of NAHR. Our results suggest that other, as yet unidentified, genetic or environmental factors play a significant role in the regulation of NAHR and are responsible for the extensive variation in the population for the probability of fathering a child with a genomic disorder resulting from a pathogenic deletion.     

Analysis of Genotype × Environment Interaction Effects for Starch Pasting Viscosity Characteristics in Indica Rice

Milled rice (Oryza sativa L.) is composed of approximately 90% starch. The properties of starch have considerable effects on cooked rice palatability and consumer acceptability. Starch pasting viscosity parameters serve as important indices in the estimation of eating, cooking, and processing qualities of rice. In the present study, four cytoplasmic male-sterile (CMS) lines and eight restorer (R) lines have been used in an incomplete diallel cross to analyze seed effects, cytoplasmic effects, maternal gene effects, and their genotype × environment (GE) effects on the following starch pasting viscosity parameters: breakdown (BD), consistency (CS), and setback (SB). The results demonstrated that the total main genetic variances (V_G) accounted for over 64% of the total genetic variance (V_G + V_GE) for the three traits, indicating that these traits were mainly controlled by the main genetic effects in addition to the GE interaction effects. The estimated total narrow-sense heritability were 67.8%, 79.5%, and 79.5% for BD, CS, and SB, respectively. The general heritability (h²_G) accounted for over 75% of the total heritability (h²_G + h²_GE), indicating that early selection would be effective for those traits and the selection efficiencies were relatively stable in different environments.