Tools for efficient epistasis detection in genome-wide association study

<p>Abstract</p> <p>Background</p> <p>Genome-wide association study (GWAS) aims to find genetic factors underlying complex phenotypic traits, for which epistasis or gene-gene interaction detection is often preferred over single-locus approach. However, the computational burden has been a major hurdle to apply epistasis test in the genome-wide scale due to a large number of single nucleotide polymorphism (SNP) pairs to be tested.</p> <p>Results</p> <p>We have developed a set of three efficient programs, FastANOVA, COE and TEAM, that support epistasis test in a variety of problem settings in GWAS. These programs utilize permutation test to properly control error rate such as family-wise error rate (FWER) and false discovery rate (FDR). They guarantee to find the optimal solutions, and significantly speed up the process of epistasis detection in GWAS.</p> <p>Conclusions</p> <p>A web server with user interface and source codes are available at the website <url>http://www.csbio.unc.edu/epistasis/</url>. The source codes are also available at SourceForge <url>http://sourceforge.net/projects/epistasis/</url>.</p

Characterization of ftsZ Mutations that Render Bacillus subtilis Resistant to MinC

Background: Cell division in Bacillus subtilis occurs precisely at midcell. Positional control of cell division is exerted by two mechanisms: nucleoid occlusion, through Noc, which prevents division through nucleoids, and the Min system, where the combined action of the MinC, D and J proteins prevents formation of the FtsZ ring at cell poles or recently completed division sites. Methodology/Principal Findings: We used a genetic screen to identify mutations in ftsZ that confer resistance to the lethal overexpression of the MinC/MinD division inhibitor. The FtsZ mutants were purified and found to polymerize to a similar or lesser extent as wild type FtsZ, and all mutants displayed reduced GTP hydrolysis activity indicative of a reduced polymerization turnover. We found that even though the mutations conferred in vivo resistance to MinC/D, the purified FtsZ mutants did not display strong resistance to MinC in vitro. Conclusions/Significance: Our results show that in B. subtilis, overproduction of MinC can be countered by mutations that alter FtsZ polymerization dynamics. Even though it would be very likely that the FtsZ mutants found depend on other Z-ring stabilizing proteins such as ZapA, FtsA or SepF, we found this not to be the case. This indicates that the cell division process in B. subtilis is extremely robust.

On the Use of Variance per Genotype as a Tool to Identify Quantitative Trait Interaction Effects: A Report from the Women's Genome Health Study

Testing for genetic effects on mean values of a quantitative trait has been a very successful strategy. However, most studies to date have not explored genetic effects on the variance of quantitative traits as a relevant consequence of genetic variation. In this report, we demonstrate that, under plausible scenarios of genetic interaction, the variance of a quantitative trait is expected to differ among the three possible genotypes of a biallelic SNP. Leveraging this observation with Levene's test of equality of variance, we propose a novel method to prioritize SNPs for subsequent gene–gene and gene–environment testing. This method has the advantageous characteristic that the interacting covariate need not be known or measured for a SNP to be prioritized. Using simulations, we show that this method has increased power over exhaustive search under certain conditions. We further investigate the utility of variance per genotype by examining data from the Women's Genome Health Study. Using this dataset, we identify new interactions between the LEPR SNP rs12753193 and body mass index in the prediction of C-reactive protein levels, between the ICAM1 SNP rs1799969 and smoking in the prediction of soluble ICAM-1 levels, and between the PNPLA3 SNP rs738409 and body mass index in the prediction of soluble ICAM-1 levels. These results demonstrate the utility of our approach and provide novel genetic insight into the relationship among obesity, smoking, and inflammation

Estrogen receptor alpha gene polymorphism and endometrial cancer risk – a case-control study

Background: Estrogen is an established endometrial carcinogen. One of the most important mediators of estrogenic action is the estrogen receptor alpha. We have investigated whether polymorphic variation in the estrogen receptor alpha gene (ESR1) is associated with endometrial cancer risk. Methods: In 702 cases with invasive endometrial cancer and 1563 controls, we genotyped five markers in ESR1 and used logistic regression models to estimate odds ratios (OR) and 95 percent confidence intervals (CI). Results: We found an association between rs2234670, rs2234693, as well as rs9340799, markers in strong linkage disequilibrium (LD), and endometrial cancer risk. The association with rs9340799 was the strongest, OR 0.75 (CI 0.60–0.93) for heterozygous and OR 0.53 (CI 0.37–0.77) for homozygous rare compared to those homozygous for the most common allele. Haplotype models did not fit better to the data than single marker models. Conclusion: We found that intronic variation in ESR1 was associated with endometrial cancer risk

Genetic Structure of the Polymorphic Metrosideros (Myrtaceae) Complex in the Hawaiian Islands Using Nuclear Microsatellite Data

Five species of Metrosideros (Myrtaceae) are recognized in the Hawaiian Islands, including the widespread M. polymorpha, and are characterized by a multitude of distinctive, yet overlapping, habit, ecological, and morphological forms. It remains unclear, despite several previous studies, whether the morphological variation within Hawaiian Metrosideros is due to hybridization, genetic polymorphism, phenotypic plasticity, or some combination of these processes. The Hawaiian Metrosideros complex has become a model system to study ecology and evolution; however this is the first study to use microsatellite data for addressing inter-island patterns of variation from across the Hawaiian Islands.Ten nuclear microsatellite loci were genotyped from 143 individuals of Metrosideros. We took advantage of the bi-parental inheritance and rapid mutation rate of these data to examine the validity of the current taxonomy and to investigate whether Metrosideros plants from the same island are more genetically similar than plants that are morphologically similar. The Bayesian algorithm of the program structure was used to define genetic groups within Hawaiian Metrosideros and the closely related taxon M. collina from the Marquesas and Austral Islands. Several standard and nested AMOVAs were conducted to test whether the genetic diversity is structured geographically or taxonomically.The results suggest that Hawaiian Metrosideros have dynamic gene flow, with genetic and morphological diversity structured not simply by geography or taxonomy, but as a result of parallel evolution on islands following rampant island-island dispersal, in addition to ancient chloroplast capture. Results also suggest that the current taxonomy requires major revisions in order to reflect the genetic structure revealed in the microsatellite data

Riverbed sediments buffer phosphorus concentrations downstream of sewage treatment works across the River Wensum catchment, UK

Purpose: Wastewater effluent discharged into rivers from sewage treatment works (STWs) represents one of the most important point sources of soluble reactive phosphorus (SRP) pollution and is a major driver of freshwater eutrophication. In this study, we assess the ability of riverbed sediments to act as a self-regulating buffering system to reduce SRP dissolved in the water column downstream of STW outflows. Materials and methods: River water and riverbed sediment samples were collected from 10 tributary outlets across the River Wensum catchment, Norfolk, UK, at monthly intervals between July and October 2016, such that 40 sediment and 40 water samples were collected in total. Of these locations, five were located downstream of STWs and five were on tributaries without STWs. Dissolved SRP concentrations were analysed and the Equilibrium Phosphorus Concentration (EPC0) of each sediment sample was measured to determine whether riverbed sediments were acting as net sources or sinks of SRP. Results and discussion: The mean SRP concentration downstream of STWs (382 µg P L-1) was double that of sites without a STW (185 µg P L-1), whilst the mean EPC0 for effluent impacted sites (105 µg P L-1) was 70% higher than that recorded at unaffected sites (62 µg P L-1). Regardless of STW influence, riverbed sediments across all 10 sites almost always acted as net sinks for SRP from the overlying water column. This was particularly true at sites downstream of STWs which displayed enhanced potential to buffer the river against increases in SRP released in sewage effluent. Conclusions: Despite EPC0 values revealing riverbed sediments were consistently acting as sinks for SRP, elevated SRP concentrations downstream of STWs clearly demonstrate the sediments have insufficient SRP sorption capacity to completely buffer the river against effluent discharge. Consequently, SRP concentrations across the catchment continue to exceed recommended standards for good chemical status, thus emphasising the need for enhanced mitigation efforts at STWs to minimise riverine phosphorus loading

Comparison of Strategies to Detect Epistasis from eQTL Data

Genome-wide association studies have been instrumental in identifying genetic variants associated with complex traits such as human disease or gene expression phenotypes. It has been proposed that extending existing analysis methods by considering interactions between pairs of loci may uncover additional genetic effects. However, the large number of possible two-marker tests presents significant computational and statistical challenges. Although several strategies to detect epistasis effects have been proposed and tested for specific phenotypes, so far there has been no systematic attempt to compare their performance using real data. We made use of thousands of gene expression traits from linkage and eQTL studies, to compare the performance of different strategies. We found that using information from marginal associations between markers and phenotypes to detect epistatic effects yielded a lower false discovery rate (FDR) than a strategy solely using biological annotation in yeast, whereas results from human data were inconclusive. For future studies whose aim is to discover epistatic effects, we recommend incorporating information about marginal associations between SNPs and phenotypes instead of relying solely on biological annotation. Improved methods to discover epistatic effects will result in a more complete understanding of complex genetic effects

Design considerations in a sib-pair study of linkage for susceptibility loci in cancer

<p>Abstract</p> <p>Background</p> <p>Modern approaches to identifying new genes associated with disease allow very fine analysis of associaton and can be performed in population based case-control studies. However, the sibpair design is still valuable because it requires few assumptions other than acceptably high penetrance to identify genetic loci.</p> <p>Methods</p> <p>We conducted simulation studies to assess the impact of design factors on relative efficiency for a linkage study of colorectal cancer. We considered two test statistics, one comparing the mean IBD probability in affected pairs to its null value of 0.5, and one comparing the mean IBD probabilities between affected and discordant pairs. We varied numbers of parents available, numbers of affected and unaffected siblings, reconstructing the genotype of an unavailable affected sibling by a spouse and offspring, and elimination of sibships where the proband carries a mutation at another locus.</p> <p>Results</p> <p>Power and efficiency were most affected by the number of affected sibs, the number of sib pairs genotyped, and the risk attributable to linked and unlinked loci. Genotyping unaffected siblings added little power for low penetrance models, but improved validity of tests when there was genetic heterogeneity and for multipoint testing. The efficiency of the concordant-only test was nearly always better than the concordant-discordant test. Replacement of an unavailable affected sibling by a spouse and offspring recovered some linkage information, particularly if several offspring were available. In multipoint analysis, the concordant-only test was showed a small anticonservative bias at 5 cM, while the multipoint concordant-discordant test was generally the most powerful test, and was not biased away from the null at 5 cM.</p> <p>Conclusion</p> <p>Genotyping parents and unaffected siblings is useful for detecting genotyping errors and if allele frequencies are uncertain. If adequate allele frequency data are available, we suggest a single-point affecteds-only analysis for an initial scan, followed by a multipoint analysis of affected and unaffected members of all available sibships with additional markers around initial hits.</p

Reciprocal Sign Epistasis between Frequently Experimentally Evolved Adaptive Mutations Causes a Rugged Fitness Landscape

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors