Variance Heterogeneity in Saccharomyces cerevisiae Expression Data: Trans-Regulation and Epistasis

Here, we describe the results from the first variance heterogeneity Genome Wide Association Study (VGWAS) on yeast expression data. Using this forward genetics approach, we show that the genetic regulation of gene-expression in the budding yeast, Saccharomyces cerevisiae, includes mechanisms that can lead to variance heterogeneity in the expression between genotypes. Additionally, we performed a mean effect association study (GWAS). Comparing the mean and variance heterogeneity analyses, we find that the mean expression level is under genetic regulation from a larger absolute number of loci but that a higher proportion of the variance controlling loci were trans-regulated. Both mean and variance regulating loci cluster in regulatory hotspots that affect a large number of phenotypes; a single variance-controlling locus, mapping close to DIA2, was found to be involved in more than 10% of the significant associations. It has been suggested in the literature that variance-heterogeneity between the genotypes might be due to genetic interactions. We therefore screened the multi-locus genotype-phenotype maps for several traits where multiple associations were found, for indications of epistasis. Several examples of two and three locus genetic interactions were found to involve variance-controlling loci, with reports from the literature corroborating the functional connections between the loci. By using a new analytical approach to re-analyze a powerful existing dataset, we are thus able to both provide novel insights to the genetic mechanisms involved in the regulation of gene-expression in budding yeast and experimentally validate epistasis as an important mechanism underlying genetic variance-heterogeneity between genotypes.     

Statistical epistasis is a generic feature of gene regulatory networks

Functional dependencies between genes are a defining characteristic of gene networks underlying quantitative traits. However, recent studies show that the proportion of the genetic variation that can be attributed to statistical epistasis varies from almost zero to very high. It is thus of fundamental as well as instrumental importance to better understand whether different functional dependency patterns among polymorphic genes give rise to distinct statistical interaction patterns or not. Here we address this issue by combining a quantitative genetic model approach with genotype-phenotype models capable of translating allelic variation and regulatory principles into phenotypic variation at the level of gene expression. We show that gene regulatory networks with and without feedback motifs can exhibit a wide range of possible statistical genetic architectures with regard to both type of effect explaining phenotypic variance and number of apparent loci underlying the observed phenotypic effect. Although all motifs are capable of harboring significant interactions, positive feedback gives rise to higher amounts and more types of statistical epistasis. The results also suggest that the inclusion of statistical interaction terms in genetic models will increase the chance to detect additional QTL as well as functional dependencies between genetic loci over a broad range of regulatory regimes. This article illustrates how statistical genetic methods can fruitfully be combined with nonlinear systems dynamics to elucidate biological issues beyond reach of each methodology in isolation.     

Phenotypic evolution from genetic polymorphisms in a radial network architecture

Background  

The genetic architecture of a quantitative trait influences the phenotypic response to natural or artificial selection. One of the main objectives of genetic mapping studies is to identify the genetic factors underlying complex traits and understand how they contribute to phenotypic expression. Presently, we are good at identifying and locating individual loci with large effects, but there is a void in describing more complex genetic architectures. Although large networks of connected genes have been reported, there is an almost complete lack of information on how polymorphisms in these networks contribute to phenotypic variation and change. To date, most of our understanding comes from theoretical, model-based studies, and it remains difficult to assess how realistic their conclusions are as they lack empirical support.     

The origin of the domestic pig: independent domestication and subsequent introgression

The domestic pig originates from the Eurasian wild boar (Sus scrofa). We have sequenced mitochondrial DNA and nuclear genes from wild and domestic pigs from Asia and Europe. Clear evidence was obtained for domestication to have occurred independently from wild boar subspecies in Europe and Asia. The time since divergence of the ancestral forms was estimated at approximately 500,000 years, well before domestication approximately 9,000 years ago. Historical records indicate that Asian pigs were introduced into Europe during the 18th and early 19th centuries. We found molecular evidence for this introgression and the data indicated a hybrid origin of some major "European" pig breeds. The study is an advance in pig genetics and has important implications for the maintenance and utilization of genetic diversity in this livestock species.     

The Multi-allelic Genetic Architecture of a Variance-Heterogeneity Locus for Molybdenum Concentration in Leaves Acts as a Source of Unexplained Additive Genetic Variance

Genome-wide association (GWA) analyses have generally been used to detect individual loci contributing to the phenotypic diversity in a population by the effects of these loci on the trait mean. More rarely, loci have also been detected based on variance differences between genotypes. Several hypotheses have been proposed to explain the possible genetic mechanisms leading to such variance signals. However, little is known about what causes these signals, or whether this genetic variance-heterogeneity reflects mechanisms of importance in natural populations. Previously, we identified a variance-heterogeneity GWA (vGWA) signal for leaf molybdenum concentrations in Arabidopsis thaliana. Here, fine-mapping of this association reveals that the vGWA emerges from the effects of three independent genetic polymorphisms that all are in strong LD with the markers displaying the genetic variance-heterogeneity. By revealing the genetic architecture underlying this vGWA signal, we uncovered the molecular source of a significant amount of hidden additive genetic variation or “missing heritability”. Two of the three polymorphisms underlying the genetic variance-heterogeneity are promoter variants for Molybdate transporter 1 (MOT1), and the third a variant located ~25 kb downstream of this gene. A fourth independent association was also detected ~600 kb upstream of MOT1. Use of a T-DNA knockout allele highlights Copper Transporter 6; COPT6 (AT2G26975) as a strong candidate gene for this association. Our results show that an extended LD across a complex locus including multiple functional alleles can lead to a variance-heterogeneity between genotypes in natural populations. Further, they provide novel insights into the genetic regulation of ion homeostasis in A. thaliana, and empirically confirm that variance-heterogeneity based GWA methods are a valuable tool to detect novel associations of biological importance in natural populations.     

Dissection of the genetic architecture of body weight in chicken reveals the impact of epistasis on domestication traits

In this contribution, we study the genetic mechanisms leading to differences in the observed growth patterns of domesticated White Leghorn chickens and their wild ancestor the red jungle fowl. An epistatic QTL analysis for several body-weight measures from hatch to adulthood confirms earlier findings that polymorphisms at >15 loci contribute to body-weight determination in an F(2) intercross between these populations and that many loci are involved in complex genetic interactions. Here, we use a new genetic model to decompose the genetic effects of this multilocus epistatic genetic network. The results show how the functional modeling of genetic effects provides new insights into how genetic interactions in a large set of loci jointly contribute to phenotypic expression. By exploring the functional effects of QTL alleles, we show that some alleles can display temporal shifts in the expression of genetic effects due to their dependencies on the genetic background. Our results demonstrate that the effects of many genes are dependent on genetic interactions with other loci and how their involvement in the domestication process relies on these interactions.     

The influence of hydrological regimes on sex ratios and spatial segregation of the sexes in two dioecious riparian shrub species in northern Sweden

River management practices have altered the hydrological regimes of many rivers and also altered the availability of regeneration niches for riparian species. We investigated the impact of changed hydrological regimes on the sex ratios and the Spatial Segregation of the Sexes (SSS) in the dioecious species Salix myrsinifolia Salisb.–phylicifolia L. and S. lapponum L. by studying the free-flowing Vindel River and the regulated Ume River in northern Sweden. We surveyed sex ratios of these species in 12 river reaches on the Vindel River and in 17 reaches on the Ume River. In addition, we surveyed the sex and location above mean river stage of 1,002 individuals across both river systems to investigate the SSS of both species. Cuttings were collected from male and female individuals of S. myrsinifolia–phylicifolia from both rivers and subjected to four different water table regimes in a greenhouse experiment to investigate growth response between the sexes. We found an M/F sex ratio in both river systems similar to the regional norm of 0.62 for S. myrsinifolia–phylicifolia and of 0.42 for S. lapponum. We found no evidence of SSS in either the free-flowing Vindel River or the regulated Ume River. In the greenhouse experiment, hydrological regime had a significant effect on shoot and root dry weight and on root length. Significantly higher shoot dry weights were found in females than in males and significantly different shoot and root dry weights were found between cuttings taken from the two rivers. We concluded that changed hydrological regimes are likely to alter dimensions of the regeneration niche and therefore to influence sex ratios and SSS at an early successional stage, making it difficult to find clear spatial patterns once these species reach maturity and can be sexed.     

Natural CMT2 Variation Is Associated With Genome-Wide Methylation Changes and Temperature Seasonality

As Arabidopsis thaliana has colonized a wide range of habitats across the world it is an attractive model for studying the genetic mechanisms underlying environmental adaptation. Here, we used public data from two collections of A. thaliana accessions to associate genetic variability at individual loci with differences in climates at the sampling sites. We use a novel method to screen the genome for plastic alleles that tolerate a broader climate range than the major allele. This approach reduces confounding with population structure and increases power compared to standard genome-wide association methods. Sixteen novel loci were found, including an association between Chromomethylase 2 (CMT2) and temperature seasonality where the genome-wide CHH methylation was different for the group of accessions carrying the plastic allele. Cmt2 mutants were shown to be more tolerant to heat-stress, suggesting genetic regulation of epigenetic modifications as a likely mechanism underlying natural adaptation to variable temperatures, potentially through differential allelic plasticity to temperature-stress.     

The Chromosomal Association of the Smc5/6 Complex Depends on Cohesion and Predicts the Level of Sister Chromatid Entanglement

The cohesin complex, which is essential for sister chromatid cohesion and chromosome segregation, also inhibits resolution of sister chromatid intertwinings (SCIs) by the topoisomerase Top2. The cohesin-related Smc5/6 complex (Smc5/6) instead accumulates on chromosomes after Top2 inactivation, known to lead to a buildup of unresolved SCIs. This suggests that cohesin can influence the chromosomal association of Smc5/6 via its role in SCI protection. Using high-resolution ChIP-sequencing, we show that the localization of budding yeast Smc5/6 to duplicated chromosomes indeed depends on sister chromatid cohesion in wild-type and top2-4 cells. Smc5/6 is found to be enriched at cohesin binding sites in the centromere-proximal regions in both cell types, but also along chromosome arms when replication has occurred under Top2-inhibiting conditions. Reactivation of Top2 after replication causes Smc5/6 to dissociate from chromosome arms, supporting the assumption that Smc5/6 associates with a Top2 substrate. It is also demonstrated that the amount of Smc5/6 on chromosomes positively correlates with the level of missegregation in top2-4, and that Smc5/6 promotes segregation of short chromosomes in the mutant. Altogether, this shows that the chromosomal localization of Smc5/6 predicts the presence of the chromatid segregation-inhibiting entities which accumulate in top2-4 mutated cells. These are most likely SCIs, and our results thus indicate that, at least when Top2 is inhibited, Smc5/6 facilitates their resolution.