Despite the importance of secondary dormancy for plant life cycle timing and survival, there is insufficient knowledge about the (epigenetic) regulation of this trait at the molecular level. Our aim was to determine the role of (epi)genetic processes in the regulation of secondary seed dormancy using natural genotypes of the widely distributed Capsella bursa-pastoris. Seeds of nine ecotypes were exposed to control conditions or histone deacetylase inhibitors [trichostatin A (TSA), valproic acid] during imbibition to study the effects of hyper-acetylation on secondary seed dormancy induction and germination. Valproic acid increased secondary dormancy and both compounds caused a delay of t50 for germination (radicle emergence) but not of t50 for testa rupture, demonstrating that they reduced speed of germination. Transcriptome analysis of one accession exposed to valproic acid versus water showed mixed regulation of ABA, negative regulation of GAs, BRs and auxins, as well as up-regulation of SNL genes, which might explain the observed delay in germination and increase in secondary dormancy. In addition, two accessions differing in secondary dormancy depth (deep vs non-deep) were studied using RNA-seq to reveal the potential regulatory processes underlying this trait. Phytohormone synthesis or signalling was generally up-regulated for ABA (e.g. NCED6, NCED2, ABCG40, ABI3) and down-regulated for GAs (GA20ox1, GA20ox2, bHLH93), ethylene (ACO1, ERF4-LIKE, ERF105, ERF109-LIKE), BRs (BIA1, CYP708A2-LIKE, probable WRKY46, BAK1, BEN1, BES1, BRI1) and auxin (GH3.3, GH3.6, ABCB19, TGG4, AUX1, PIN6, WAT1). Epigenetic candidates for variation in secondary dormancy depth include SNL genes, histone deacetylases and associated genes (HDA14, HDA6-LIKE, HDA-LIKE, ING2, JMJ30), as well as sequences linked to histone acetyltransferases (bZIP11, ARID1A-LIKE), or to gene silencing through histone methylation (SUVH7, SUVH9, CLF). Together, these results show that phytohormones and epigenetic regulation play an important role in controlling differences in secondary dormancy depth between accessions.

The classical twin design relies on a number of strong number of assumptions in order to yield unbiased estimates of heritability. This includes the equal environments assumption — that monozygotic and dizygotic twins experience similar degrees of environmental similarity — an assumption that is likely to be violated in practice for many traits of interest. An alternative method of estimating heritability that does not suffer from many of these limitations is to model trait similarity between sibling pairs as a function of their empirical genome-wide identity by descent sharing, estimated from genetic markers. In this review, I recount the story behind Nick Martin’s and my development of this method, our first attempts at applying it in a human population and more recent studies using the original and related methods to estimate trait heritability.

Monozygotic (MZ) twins form an important system for the study of biological plasticity in humans. While MZ twins are generally considered to be genetically identical, a number of studies have emerged that have demonstrated copy-number differences within a twin pair, particularly in those discordant for disease. The rate of autosomal copy-number variation (CNV) discordance within MZ twin pairs was investigated using a population sample of 376 twin pairs genotyped on Illumina Human610-Quad arrays. After CNV calling using both QuantiSNP and PennCNV followed by manual annotation, only a single CNV difference was observed within the MZ twin pairs, being a 130 KB duplication of chromosome 5. Five other potential discordant CNV were called by the software, but excluded based on manual annotation of the regions. It is concluded that large CNV discordance is rare within MZ twin pairs, indicating that any CNV difference found within phenotypically discordant MZ twin pairs has a high probability of containing the causal gene(s) involved.

Genome-wide association analysis on monozygotic twin-pairs offers a route to discovery of gene–environment interactions through testing for variability loci associated with sensitivity to individual environment/lifestyle. We present a genome-wide scan of loci associated with intra-pair differences in serum lipid and apolipoprotein levels. We report data for 1,720 monozygotic female twin-pairs from GenomEUtwin project with 2.5 million SNPs, imputed or genotyped, and measured serum lipid fractions for both twins. We found one locus associated with intra-pair differences in high-density lipoprotein cholesterol, rs2483058 in an intron of SRGAP2, where twins carrying the C allele are more sensitive to environmental factors (P = 3.98 × 10−8). We followed up the association in further genotyped monozygotic twins (N = 1,261), which showed a moderate association for the variant (P = 0.200, same direction of an effect). In addition, we report a new association on the level of apolipoprotein A-II (P = 4.03 × 10−8).

Genome wide association studies (GWAS) have largely succeeded family-based linkage studies in livestock and human populations as the preferred method to map loci for complex or quantitative traits. However, the type of results produced by the two analyses contrast sharply due to differences in linkage disequilibrium (LD) imposed by the design of studies. In this paper, we demonstrate that association and linkage studies are in agreement provided that (i) the effects from both studies are estimated appropriately as random effects, (ii) all markers are fitted simultaneously and (iii) appropriate adjustments are made for the differences in LD between the study designs. We demonstrate with real data that linkage results can be predicted by the sum of association effects. Our association study captured most of the linkage information because we could predict the linkage results with moderate accuracy. We suggest that the ability of common single nucleotide polymorphism (SNP) to capture the genetic variance in a population will depend on the effective population size of the study organism. The results provide further evidence for many loci of small effect underlying complex traits. The analysis suggests a more informed method for GWAS is to fit statistical models where all SNPs are analysed simultaneously and as random effects.

When using maximum likelihood methods to estimate genetic and environmental components of (co)variance, it is common to test hypotheses using likelihood ratio tests, since such tests have desirable asymptotic properties. In particular, the standard likelihood ratio test statistic is assumed asymptotically to follow a χ2 distribution with degrees of freedom equal to the number of parameters tested. Using the relationship between least squares and maximum likelihood estimators for balanced designs, it is shown why the asymptotic distribution of the likelihood ratio test for variance components does not follow a χ2 distribution with degrees of freedom equal to the number of parameters tested when the null hypothesis is true. Instead, the distribution of the likelihood ratio test is a mixture of χ2 distributions with different degrees of freedom. Implications for testing variance components in twin designs and for quantitative trait loci mapping are discussed. The appropriate distribution of the likelihood ratio test statistic should be used in hypothesis testing and model selection.