Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L.

CONTEXT: Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited. OBJECTIVE: To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts. METHOD: We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed. RESULTS: We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH. CONCLUSIONS: This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Statistical Power and the Classical Twin Design.

Dr Nick Martin has made enormous contributions to the field of behavior genetics over the past 50 years. Of his many seminal papers that have had a profound impact, we focus on his early work on the power of twin studies. He was among the first to recognize the importance of sample size calculation before conducting a study to ensure sufficient power to detect the effects of interest. The elegant approach he developed, based on the noncentral chi-squared distribution, has been adopted by subsequent researchers for other genetic study designs, and today remains a standard tool for power calculations in structural equation modeling and other areas of statistical analysis. The present brief article discusses the main aspects of his seminal paper, and how it led to subsequent developments, by him and others, as the field of behavior genetics evolved into the present era.

Nicholas (Nick) G. Martin and the Extended Twin Model.

The extended twin model is a unique design in the genetic epidemiology toolbox that allows to simultaneously estimate multiple causes of variation such as genetic and cultural transmission, genotype-environment covariance and assortative mating, among others. Nick Martin has played a key role in the conception of the model, the collection of substantially large data sets to test the model, the application of the model to a range of phenotypes, the publication of the results including cross-cultural comparisons, the evaluation of bias and power of the design and the further elaborations of the model, such as the children-of-twins design.

A Brief History of the Collaboration between Dr Nathan Gillespie and Professor Nick Martin Including Personal Reflections.

This article describes Dr Nathan Gillespie's PhD training and supervision under Professor Nick Martin and their ongoing collaborations. Drs Gillespie and Martin have collaborated on numerous biometrical genetic analyses applied to cross-sectional and longitudinal twin data, combined molecular and phenotypic modeling, as well as genomewide meta-analyses of psychoactive substance use and misuse. Dr Gillespie remains an active collaborator with Professor Martin, including ongoing data collection, analysis and publications related to the Brisbane Longitudinal Twin Study.

Twin studies of subjective sleep quality and sleep duration, and their behavioral correlates: Systematic review and meta-analysis of heritability estimates.

Twin studies have shown that a substantial proportion of the variance for sleep variables is due to genetic factors. However, there is still considerable heterogeneity among research reports. Our main objectives were to: 1) Review the twin literature regarding sleep quality and duration, as well as their behavioural correlates; 2) Estimate the mean heritability of subjective sleep quality and sleep duration; 3) Assess heterogeneity among studies on these topics; and 4) Search for moderator variables. Two parallel meta-analyses were carried out for sleep quality and sleep duration. Seventeen articles were included in the meta-analysis. Mean MZ correlations were consistently higher than DZ correlations. A mean heritability of 0.31 (95% CI: 0.20, 0.41) was found for subjective sleep quality (range: 0-0.43) and 0.38 (95% CI: 0.16, 0.56) for sleep duration (range: 0-1). Heterogeneity indexes were significant for both sleep quality (I2 = 98.77, p < .001) and sleep duration (I2 = 99.73, p < .001). The high heterogeneity warrants further research considering possible moderators that may affect heritability.

The Chinese National Twin Registry: a 'gold mine' for scientific research.

The Chinese National Twin Registry (CNTR) currently includes data from 61 566 twin pair from 11 provinces or cities in China. Of these, 31 705, 15 060 and 13 531 pairs are monozygotic, same-sex dizygotic and opposite-sex dizygotic pairs, respectively, determined by opposite sex or intrapair similarity. Since its establishment in 2001, the CNTR has provided an important resource for analysing genetic and environmental influences on chronic diseases especially cardiovascular diseases. Recently, the CNTR has focused on collecting biologic specimens from disease-concordant or disease-discordant twin pairs or from twin pairs reared apart. More than 8000 pairs of these twins have been registered, and blood samples have been collected from more than 1500 pairs. In this review, we summarize the main findings from univariate and multivariate genetic effects analyses, gene-environment interaction studies, omics studies exploring DNA methylation and metabolomic markers associated with phenotypes. There remains further scope for CNTR research and data mining. The plan for future development of the CNTR is described. The CNTR welcomes worldwide collaboration.

Type I Error Rates and Parameter Bias in Multivariate Behavioral Genetic Models.

For many multivariate twin models, the numerical Type I error rates are lower than theoretically expected rates using a likelihood ratio test (LRT), which implies that the significance threshold for statistical hypothesis tests is more conservative than most twin researchers realize. This makes the numerical Type II error rates higher than theoretically expected. Furthermore, the discrepancy between the observed and expected error rates increases as more variables are included in the analysis and can have profound implications for hypothesis testing and statistical inference. In two simulation studies, we examine the Type I error rates for the Cholesky decomposition and Correlated Factors models. Both show markedly lower than nominal Type I error rates under the null hypothesis, a discrepancy that increases with the number of variables in the model. In addition, we observe slightly biased parameter estimates for the Cholesky decomposition and Correlated Factors models. By contrast, if the variance-covariance matrices for variance components are estimated directly (without constraints), the numerical Type I error rates are consistent with theoretical expectations and there is no bias in the parameter estimates regardless of the number of variables analyzed. We call this the direct symmetric approach. It appears that each model-implied boundary, whether explicit or implicit, increases the discrepancy between the numerical and theoretical Type I error rates by truncating the sampling distributions of the variance components and inducing bias in the parameters. The direct symmetric approach has several advantages over other multivariate twin models as it corrects the Type I error rate and parameter bias issues, is easy to implement in current software, and has fewer optimization problems. Implications for past and future research, and potential limitations associated with direct estimation of genetic and environmental covariance matrices are discussed.

Limited common origins of multiple adult health-related behaviors: Evidence from U.S. twins.

Health-related behaviors are significant contributors to morbidity and mortality in the United States, yet evidence on the underlying causes of the vast within-population variation in behaviors is mixed. While many potential causes of health-related behaviors have been identified-such as schooling, genetics, and environments-little is known on how much of the variation across multiple behaviors is due to a common set of causes. We use three separate datasets on U.S. twins to investigate the degree to which multiple health-related behaviors correlate and can be explained by a common set of factors. We find that aside from smoking and drinking, most behaviors are not strongly correlated among individuals. Based on the results of both within-identical-twins regressions and multivariate behavioral genetics models, we find some evidence that schooling may be related to smoking but not to the covariation between multiple behaviors. Similarly, we find that a large fraction of the variance in each of the behaviors is consistent with genetic factors; however, we do not find strong evidence that a single common set of genes explains variation in multiple behaviors. We find, however, that a large portion of the correlation between smoking and heavy drinking is consistent with common, mostly childhood, environments. This suggests that the initiation and patterns of these two behaviors might arise from a common childhood origin. Research and policy to identify and modify this source may provide a strong way to reduce the population health burden of smoking and heavy drinking.

The latent class twin method.

The twin method refers to the use of data from same-sex identical and fraternal twins to estimate the genetic and environmental contributions to a trait or outcome. The standard twin method is the variance component twin method that estimates heritability, the fraction of variance attributed to additive genetic inheritance. The latent class twin method estimates two quantities that are easier to interpret than heritability: the genetic prevalence, which is the fraction of persons in the genetic susceptibility latent class, and the heritability fraction, which is the fraction of persons in the genetic susceptibility latent class with the trait or outcome. We extend the latent class twin method in three important ways. First, we incorporate an additive genetic model to broaden the sensitivity analysis beyond the original autosomal dominant and recessive genetic models. Second, we specify a separate survival model to simplify computations and improve convergence. Third, we show how to easily adjust for covariates by extending the method of propensity scores from a treatment difference to zygosity. Applying the latent class twin method to data on breast cancer among Nordic twins, we estimated a genetic prevalence of 1%, a result with important implications for breast cancer prevention research.

Marginal genetic effects estimation in family and twin studies using random-effects models.

Random-effects models are often used in family-based genetic association studies to properly capture the within families relationships. In such models, the regression parameters have a conditional on the random effects interpretation and they measure, e.g., genetic effects for each family. Estimating parameters that can be used to make inferences at the population level is often more relevant than the family-specific effects, but not straightforward. This is mainly for two reasons: First the analysis of family data often requires high-dimensional random-effects vectors to properly model the familial relationships, for instance when members with a different degree of relationship are considered, such as trios, mix of monozygotic and dizygotic twins, etc. The second complication is the biased sampling design, such as the multiple cases families design, which is often employed to enrich the sample with genetic information. For these reasons deriving parameters with the desired marginal interpretation can be challenging. In this work we consider the marginalized mixed-effects models, we discuss challenges in applying them in ascertained family data and propose penalized maximum likelihood methodology to stabilize the parameter estimation by using external information on the disease prevalence or heritability. The performance of our methodology is evaluated via simulation and is illustrated on data from Rheumatoid Arthritis patients, where we estimate the marginal effect of HLA-DRB1*13 and shared epitope alleles across three different study designs and combine them using meta-analysis.

Operating Characteristics of Statistical Methods for Detecting Gene-by-Measured Environment Interaction in the Presence of Gene-Environment Correlation under Violations of Distributional Assumptions.

Accurately identifying interactions between genetic vulnerabilities and environmental factors is of critical importance for genetic research on health and behavior. In the previous work of Van Hulle et al. (Behavior Genetics, Vol. 43, 2013, pp. 71-84), we explored the operating characteristics for a set of biometric (e.g., twin) models of Rathouz et al. (Behavior Genetics, Vol. 38, 2008, pp. 301-315), for testing gene-by-measured environment interaction (GxM) in the presence of gene-by-measured environment correlation (rGM) where data followed the assumed distributional structure. Here we explore the effects that violating distributional assumptions have on the operating characteristics of these same models even when structural model assumptions are correct. We simulated N = 2,000 replicates of n = 1,000 twin pairs under a number of conditions. Non-normality was imposed on either the putative moderator or on the ultimate outcome by ordinalizing or censoring the data. We examined the empirical Type I error rates and compared Bayesian information criterion (BIC) values. In general, non-normality in the putative moderator had little impact on the Type I error rates or BIC comparisons. In contrast, non-normality in the outcome was often mistaken for or masked GxM, especially when the outcome data were censored.

VEGAS2: Software for More Flexible Gene-Based Testing.

Gene-based tests such as versatile gene-based association study (VEGAS) are commonly used following per-single nucleotide polymorphism (SNP) GWAS (genome-wide association studies) analysis. Two limitations of VEGAS were that the HapMap2 reference set was used to model the correlation between SNPs and only autosomal genes were considered. HapMap2 has now been superseded by the 1,000 Genomes reference set, and whereas early GWASs frequently ignored the X chromosome, it is now commonly included. Here we have developed VEGAS2, an extension that uses 1,000 Genomes data to model SNP correlations across the autosomes and chromosome X. VEGAS2 allows greater flexibility when defining gene boundaries. VEGAS2 offers both a user-friendly, web-based front end and a command line Linux version. The online version of VEGAS2 can be accessed through https://vegas2.qimrberghofer.edu.au/. The command line version can be downloaded from https://vegas2.qimrberghofer.edu.au/zVEGAS2offline.tgz. The command line version is developed in Perl, R and shell scripting languages; source code is available for further development.

Estimating heritability for cause specific mortality based on twin studies.

There has been considerable interest in studying the magnitude and type of inheritance of specific diseases. This is typically derived from family or twin studies, where the basic idea is to compare the correlation for different pairs that share different amount of genes. We here consider data from the Danish twin registry and discuss how to define heritability for cancer occurrence. The key point is that this should be done taking censoring as well as competing risks due to e.g.  death into account. We describe the dependence between twins on the probability scale and show that various models can be used to achieve sensible estimates of the dependence within monozygotic and dizygotic twin pairs that may vary over time. These dependence measures can subsequently be decomposed into a genetic and environmental component using random effects models. We here present several novel models that in essence describe the association in terms of the concordance probability, i.e., the probability that both twins experience the event, in the competing risks setting. We also discuss how to deal with the left truncation present in the Nordic twin registries, due to sampling only of twin pairs where both twins are alive at the initiation of the registries.

How meaningful are heritability estimates of liability?

It is commonly acknowledged that estimates of heritability from classical twin studies have many potential shortcomings. Despite this, in the post-GWAS era, these heritability estimates have come to be a continual source of interest and controversy. While the heritability estimates of a quantitative trait are subject to a number of biases, in this article we will argue that the standard statistical approach to estimating the heritability of a binary trait relies on some additional untestable assumptions which, if violated, can lead to badly biased estimates. The ACE liability threshold model assumes at its heart that each individual has an underlying liability or propensity to acquire the binary trait (e.g., disease), and that this unobservable liability is multivariate normally distributed. We investigated a number of different scenarios violating this assumption such as the existence of a single causal diallelic gene and the existence of a dichotomous exposure. For each scenario, we found that substantial asymptotic biases can occur, which no increase in sample size can remove. Asymptotic biases as much as four times larger than the true value were observed, and numerous cases also showed large negative biases. Additionally, regions of low bias occurred for specific parameter combinations. Using simulations, we also investigated the situation where all of the assumptions of the ACE liability model are met. We found that commonly used sample sizes can lead to biased heritability estimates. Thus, even if we are willing to accept the meaningfulness of the liability construct, heritability estimates under the ACE liability threshold model may not accurately reflect the heritability of this construct. The points made in this paper should be kept in mind when considering the meaningfulness of a reported heritability estimate for any specific disease.

The heritability of human disease: estimation, uses and abuses.

Relatives provide the basic material for the study of inheritance of human disease. However, the methodologies for the estimation of heritability and the interpretation of the results have been controversial. The debate arises from the plethora of methods used, the validity of the methodological assumptions and the inconsistent and sometimes erroneous genetic interpretations made. We will discuss how to estimate disease heritability, how to interpret it, how biases in heritability estimates arise and how heritability relates to other measures of familial disease aggregation.

Interpreting estimates of heritability--a note on the twin decomposition.

While most outcomes may in part be genetically mediated, quantifying genetic heritability is a different matter. To explore data on twins and decompose the variation is a classical method to determine whether variation in outcomes, e.g. IQ or schooling, originate from genetic endowments or environmental factors. Despite some criticism, the model is still widely used. The critique is generally related to how estimates of heritability may encompass environmental mediation. This aspect is sometimes left implicit by authors even though its relevance for the interpretation is potentially profound. This short note is an appeal for clarity from authors when interpreting the magnitude of heritability estimates. It is demonstrated how disregarding existing theoretical contributions can easily lead to unnecessary misinterpretations and/or controversies. The key arguments are relevant also for estimates based on data of adopted children or from modern molecular genetics research.

Adjusted confidence intervals for a bounded parameter.

It is well known that the regular likelihood ratio test of a bounded parameter is not valid if the boundary value is being tested. This is the case for testing the null value of a scalar variance component. Although an adjusted test of variance component has been suggested to account for the effect of its lower bound of zero, no adjustment of its interval estimate has ever been proposed. If left unadjusted, the confidence interval of the variance may still contain zero when the adjusted test rejects the null hypothesis of a zero variance, leading to conflicting conclusions. In this research, we propose two ways to adjust the confidence interval of a parameter subject to a lower bound, one based on the Wald test and the other on the likelihood ratio test. Both are compatible to the adjusted test and parametrization-invariant. A simulation study and two examples are given in the framework of ACDE models in twin studies.

Sibling comparison designs: bias from non-shared confounders and measurement error.

Twins, full siblings, and half-siblings are increasingly used as comparison groups in matched cohort and matched case-control studies. The "within-pair" estimates acquired through these comparisons are free from confounding from all factors that are shared by the siblings. This has made sibling comparisons popular in studying associations thought likely to suffer confounding from socioeconomic or genetic factors. Despite the wide application of these designs in epidemiology, they have received little scrutiny from a statistical or methodological standpoint. In this paper we show, analytically and through a series of simulations, that the standard interpretation of the models is subject to several limitations that are rarely acknowledged.Although within-pair estimates will not be confounded by factors shared by the siblings, such estimates are more severely biased by non-shared confounders than the unpaired estimate. If siblings are less similar with regard to confounders than to the exposure under study, the within-pair estimate will always be more biased than the ordinary unpaired estimate. Attenuation of associations due to random measurement error in exposure will also be higher in the within-pair estimate, leading within-pair associations to be weaker than corresponding unpaired associations, even in the absence of confounding. Implications for the interpretation of sibling comparison results are discussed.

Visual analysis of geocoded twin data puts nature and nurture on the map.

Twin studies allow us to estimate the relative contributions of nature and nurture to human phenotypes by comparing the resemblance of identical and fraternal twins. Variation in complex traits is a balance of genetic and environmental influences; these influences are typically estimated at a population level. However, what if the balance of nature and nurture varies depending on where we grow up? Here we use statistical and visual analysis of geocoded data from over 6700 families to show that genetic and environmental contributions to 45 childhood cognitive and behavioral phenotypes vary geographically in the United Kingdom. This has implications for detecting environmental exposures that may interact with the genetic influences on complex traits, and for the statistical power of samples recruited for genetic association studies. More broadly, our experience demonstrates the potential for collaborative exploratory visualization to act as a lingua franca for large-scale interdisciplinary research.