Signaling pathway genes for blood pressure, folate and cholesterol levels among hypertensives : an epistasis analysis

Irregular atrial pressure, defective folate and cholesterol metabolism contribute to the pathogenesis of hypertension. However, little is known about the combined roles of the methylenetetrahydrofolate reductase (MTHFR), apolipoprotein-E (ApoE) and angiotensin-converting enzyme (ACE) genes, which are involved in metabolism and homeostasis. The objective of this study is to investigate the association of the MTHFR 677 C>T and 1298A>C, ACE insertion–deletion (I/D) and ApoE genetic polymorphisms with hypertension and to further explore the epistasis interactions that are involved in these mechanisms. A total of 594 subjects, including 348 normotensive and 246 hypertensive ischemic stroke subjects were recruited. The MTHFR 677 C>T and 1298A>C, ACE I/D and ApoEpolymorphisms were genotyped and the epistasis interaction were analyzed. The MTHFR 677 C>T and ApoE polymorphisms demonstrated significant associations with susceptibility to hypertension in multiple logistic regression models, multifactor dimensionality reduction and a classification and regression tree. In addition, the logistic regression model demonstrated that significant interactions between the ApoE E3E3, E2E4, E2E2 and MTHFR 677 C>T polymorphisms existed. In conclusion, the results of this epistasis study indicated significant association between the ApoE and MTHFR polymorphisms and hypertension.

Major histocompatibility complex associations of ankylosing spondylitis are complex and involve further epistasis with ERAP1

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthritis for which HLA-B*27 is the major genetic risk factor, although its role in the aetiology of AS remains elusive. To better understand the genetic basis of the MHC susceptibility loci, we genotyped 7,264 MHC SNPs in 22,647 AS cases and controls of European descent. We impute SNPs, classical HLA alleles and amino-acid residues within HLA proteins, and tested these for association to AS status. Here we show that in addition to effects due to HLA-B*27 alleles, several other HLA-B alleles also affect susceptibility. After controlling for the associated haplotypes in HLA-B, we observe independent associations with variants in the HLA-A, HLA-DPB1 and HLA-DRB1 loci. We also demonstrate that the ERAP1 SNP rs30187 association is not restricted only to carriers of HLA-B*27 but also found in HLA-B*40:01 carriers independently of HLA-B*27 genotype.

Prevalent positive epistasis in Escherichia coli and Saccharomyces cerevisiae metabolic networks

Epistasis refers to the interaction between genes. Although high-throughput epistasis data from model organisms are being generated and used to construct genetic networks(1-3), the extent to which genetic epistasis reflects biologically meaningful interactions remains unclear(4-6). We have addressed this question through in silico mapping of positive and negative epistatic interactions amongst biochemical reactions within the metabolic networks of Escherichia coli and Saccharomyces cerevisiae using flux balance analysis. We found that negative epistasis occurs mainly between nonessential reactions with overlapping functions, whereas positive epistasis usually involves essential reactions, is highly abundant and, unexpectedly, often occurs between reactions without overlapping functions. We offer mechanistic explanations of these findings and experimentally validate them for 61 S. cerevisiae gene pairs.

Epistasis among Drosophila persimilis factors conferring hybrid male sterility with D. pseudoobscura bogotana.

The Bateson-Dobzhansky-Muller model posits that hybrid incompatibilities result from genetic changes that accumulate during population divergence. Indeed, much effort in recent years has been devoted to identifying genes associated with hybrid incompatibilities, often with limited success, suggesting that hybrid sterility and inviability are frequently caused by complex interactions between multiple loci and not by single or a small number of gene pairs. Our previous study showed that the nature of epistasis between sterility-conferring QTL in the Drosophila persimilis-D. pseudoobscura bogotana species pair is highly specific. Here, we further dissect one of the three QTL underlying hybrid male sterility between these species and provide evidence for multiple factors within this QTL. This result indicates that the number of loci thought to contribute to hybrid dysfunction may have been underestimated, and we discuss how linkage and complex epistasis may be characteristic of the genetics of hybrid incompatibilities. We further pinpoint the location of one locus that confers hybrid male sterility when homozygous, dubbed "mule-like", to roughly 250 kilobases.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Analysis of epistasis in linked regions in the Irish study of high-density schizophrenia families

Epistasis may be important in the etiology of schizophrenia. Analysis of epistasis has been important in the positional cloning of a gene involved in the etiology of type II diabetes mellitus. We investigated the importance of epistasis among six linked regions in 268 multiplex pedigrees in the Irish Study of High-Density Schizophrenia Families (ISHDSF) by computing pairwise correlations between nonparametric linkage scores for narrow, intermediate, and broad diagnostic definitions. The linked regions were on chromosomes 2, 4, 5, 6, 8, and 10. No correlation reached our a priori level of statistical significance. Using this statistical approach, we did not find evidence of important epistatic effects among these six regions in the ISHDSF.

Marker-trait association and epistasis for brown rust resistance in sugarcane

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Epistasis among eNOS, MMP-9 and VEGF maternal genotypes in hypertensive disorders of pregnancy

Polymorphisms of the endothelial nitric oxide synthase (eNOS), matrix metalloproteinase-9 (MMP-9) and vascular endothelial growth factor (VEGF) genes were shown to be associated with hypertensive disorders of pregnancy. However, epistasis is suggested to be an important component of the genetic susceptibility to preeclampsia (PE). The aim of this study was to characterize the interactions among these genes in PE and gestational hypertension (GH). Seven clinically relevant polymorphisms of eNOS (T-786C, rs2070744, a variable number of tandem repeats in intron 4 and Glu298Asp, rs1799983), MMP-9 (C-1562T, rs3918242 and -90(CA)(13-25), rs2234681) and VEGF (C-2578A, rs699947 and G-634C, rs2010963) were genotyped by TaqMan allelic discrimination assays or PCR and fragment separation by electrophoresis in 122 patients with PE, 107 patients with GH and a control group of 102 normotensive pregnant (NP) women. A robust multifactor dimensionality reduction analysis was used to characterize gene-gene interactions. Although no significant genotype combinations were observed for the comparison between the GH and NP groups (P>0.05), the combination of MMP-9-1562CC with VEGF-634GG was more frequent in NP women than in women with PE (P<0.05). Moreover, the combination of MMP-9-1562CC with VEGF-634CC or MMP-9-1562CT with VEGF-634CC or-634GG was more frequent in women with PE than in NP women (P<0.05). These results are obscured when single polymorphisms in these genes are considered and suggest that specific genotype combinations of MMP-9 and VEGF contribute to PE susceptibility. Hypertension Research (2012) 35, 917-921; doi:10.1038/hr.2012.60; published online 10 May 2012

The role of epistasis in the manifestation of heterosis: a systems-oriented approach

Heterosis is widely used in breeding, but the genetic basis of this biological phenomenon has not been elucidated. We postulate that additive and dominance genetic effects as well as two-locus interactions estimated in classical QTL analyses are not sufficient for quantifying the contributions of QTL to heterosis. A general theoretical framework for determining the contributions of different types of genetic effects to heterosis was developed. Additive x additive epistatic interactions of individual loci with the entire genetic background were identified as a major component of midparent heterosis. On the basis of these findings we defined a new type of heterotic effect denoted as augmented dominance effect di* that comprises the dominance effect at each QTL minus half the sum of additive x additive interactions with all other QTL. We demonstrate that genotypic expectations of QTL effects obtained from analyses with the design III using testcrosses of recombinant inbred lines and composite-interval mapping precisely equal genotypic expectations of midparent heterosis, thus identifying genomic regions relevant for expression of heterosis. The theory for QTL mapping of multiple traits is extended to the simultaneous mapping of newly defined genetic effects to improve the power of QTL detection and distinguish between dominance and overdominance.

Novel genetic risk markers for ulcerative colitis in the IL2/IL21 region are in epistasis with IL23R and suggest a common genetic background for ulcerative colitis and celiac disease

OBJECTIVES: Recently, a genome-wide association study showed that single-nucleotide polymorphisms (SNPs) in the chromosome 4q27 region containing IL2 and IL21 are associated with celiac disease. Given the increased prevalence of inflammatory bowel disease (IBD) among celiac disease patients, we investigated the possible involvement of these SNPs in IBD. METHODS: Five SNPs strongly associated with celiac disease within the KIAA1109/TENR/IL2/IL21 linkage disequilibrium block on chromosome 4q27 and one coding SNP within the IL21 gene were analyzed in a large German IBD cohort. The study population comprised a total of 2,948 Caucasian individuals, including 1,461 IBD patients (ulcerative colitis (UC): n=514, Crohn's disease (CD): n=947) and 1,487 healthy unrelated controls. RESULTS: Three of the five celiac disease risk markers had a protective effect on UC susceptibility, and this effect remained significant after correcting for multiple testing: rs6840978: P=0.0082, P(corr)=0.049, odds ratio (OR) 0.77, 95% confidence interval (CI) 0.63-0.93; rs6822844: P=0.0028, P(corr)=0.017, OR 0.73, 95% CI 0.59-0.90; rs13119723: P=0.0058, P(corr)=0.035, OR 0.75, 95% CI 0.61-0.92. A haplotype consisting of the six SNPs tested was markedly associated with UC susceptibility (P=0.0025, P(corr)=0.015, OR 0.72, 95% CI 0.58-0.89). Moreover, in UC, epistasis was observed between the IL23R SNP rs1004819 and three SNPs in the KIAA1109/TENR/IL2/IL21 block (rs13151961, rs13119723, and rs6822844). CONCLUSIONS: Similar to other autoimmune diseases such as celiac disease, rheumatoid arthritis, type 1 diabetes, Graves' disease, and psoriatic arthritis, genetic variation in the chromosome 4q27 region predisposes to UC, suggesting a common genetic background for these diseases.

Intergenomic epistasis causes asynchronous hatch times in whitefish hybrids, but only when parental ecotypes differ

Support for the theory of ecological speciation requires evidence for ecological divergence between species which directly or indirectly causes reproductive isolation. This study investigates effects of ecological vs. genetic disparity of parental species on the presence of endogenous selection (deformation and mortality rates) and potential sources of exogenous selection (growth rates and hatch timing) on hybrids. Hybrid embryonic development is analysed in a common-garden full-sib cross of three species belonging to two different ecotypes within the Coregonus lavaretus species flock in the central Alpine region of Europe. Although hatch timing was similar across the three species, embryonic growth rates and egg sizes differed between ecotypes. This led to a mismatch between embryonic growth rate and egg size in hybrid crosses that reveals epistasis between the maternal and embryonic genomes and transgressive hatch times that were asynchronous with control crosses. A strong constraint of egg size to embryo size at late development was also evident. We argue that this demonstrates potential for coadaptation of a maternal trait (egg size) with offspring growth rate to be an important source of selection against hybridization between ecotypes with different egg sizes. Implications for the measurement and quantification of early life-history traits affected by this additive relationship, such as hatch day and larval size, are also discussed.

Identification of novel susceptibility loci for inflammatory bowel disease on chromosomes 1p, 3q, and 4q: Evidence for epistasis between 1p and IBD1

The idiopathic inflammatory bowel diseases, Crohn’s disease (CD) and ulcerative colitis (UC), are chronic, frequently disabling diseases of the intestines. Segregation analyses, twin concordance, and ethnic differences in familial risks have established that CD and UC are complex, non-Mendelian, related genetic disorders. We performed a genome-wide screen using 377 autosomal markers, on 297 CD, UC, or mixed relative pairs from 174 families, 37% Ashkenazim. We observed evidence for linkage at 3q for all families (multipoint logarithm of the odds score (MLod) = 2.29, P = 5.7 × 10−4), with greatest significance for non-Ashkenazim Caucasians (MLod = 3.39, P = 3.92 × 10−5), and at chromosome 1p (MLod = 2.65, P = 2.4 × 10−4) for all families. In a limited subset of mixed families (containing one member with CD and another with UC), evidence for linkage was observed on chromosome 4q (MLod = 2.76, P = 1.9 × 10−4), especially among Ashkenazim. There was confirmatory evidence for a CD locus, overlapping IBD1, in the pericentromeric region of chromosome 16 (MLod = 1.69, P = 2.6 × 10−3), particularly among Ashkenazim (MLod = 1.51, P = 7.8 × 10−3); however, positive MLod scores were observed over a very broad region of chromosome 16. Furthermore, evidence for epistasis between IBD1 and chromosome 1p was observed. Thirteen additional loci demonstrated nominal (MLod > 1.0, P < 0.016) evidence for linkage. This screen provides strong evidence that there are several major susceptibility loci contributing to the genetic risk for CD and UC.

Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid

The genetic basis of heterosis was investigated in an elite rice hybrid by using a molecular linkage map with 150 segregating loci covering the entire rice genome. Data for yield and three traits that were components of yield were collected over 2 years from replicated field trials of 250 F2:3 families. Genotypic variations explained from about 50% to more than 80% of the total variation. Interactions between genotypes and years were small compared with the main effects. A total of 32 quantitative trait loci (QTLs) were detected for the four traits; 12 were observed in both years and the remaining 20 were detected in only one year. Overdominance was observed for most of the QTLs for yield and also for a few QTLs for the component traits. Correlations between marker heterozygosity and trait expression were low, indicating that the overall heterozygosity made little contribution to heterosis. Digenic interactions, including additive by additive, additive by dominance, and dominance by dominance, were frequent and widespread in this population. The interactions involved large numbers of marker loci, most of which individually were not detectable on single-locus basis; many interactions among loci were detected in both years. The results provide strong evidence that epistasis plays a major role as the genetic basis of heterosis.

Certain observations concerning the effects of epistasis on complex traits and the evolution of genomes.

Thesis (Ph.D.)--University of Washington, 2016-05

Simulating the effects of dominance and epistasis on selection response in the CIMMYT wheat breeding program using QuCim

Plant breeders use many different breeding methods to develop superior cultivars. However, it is difficult, cumbersome, and expensive to evaluate the performance of a breeding method or to compare the efficiencies of different breeding methods within an ongoing breeding program. To facilitate comparisons, we developed a QU-GENE module called QuCim that can simulate a large number of breeding strategies for self-pollinated species. The wheat breeding strategy Selected Bulk used by CIMMYT's wheat breeding program was defined in QuCim as an example of how this is done. This selection method was simulated in QuCim to investigate the effects of deviations from the additive genetic model, in the form of dominance and epistasis, on selection outcomes. The simulation results indicate that the partial dominance model does not greatly influence genetic advance compared with the pure additive model. Genetic advance in genetic systems with overdominance and epistasis are slower than when gene effects are purely additive or partially dominant. The additive gene effect is an appropriate indicator of the change in gene frequency following selection when epistasis is absent. In the absence of epistasis, the additive variance decreases rapidly with selection. However, after several cycles of selection it remains relatively fixed when epistasis is present. The variance from partial dominance is relatively small and therefore hard to detect by the covariance among half sibs and the covariance among full sibs. The dominance variance from the overdominance model can be identified successfully, but it does not change significantly, which confirms that overdominance cannot be utilized by an inbred breeding program. QuCim is an effective tool to compare selection strategies and to validate some theories in quantitative genetics.