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.

Inclusão da epistasia em modelo de avaliação genética de bovinos de corte compostos

Dados de bovinos compostos foram analisados para avaliar o efeito da epistasia nos modelos de avaliação genética. As características analisadas foram os pesos aos 205 (P205) e 390 dias (P390) e perímetro escrotal aos 390 dias (PE390). As análises foram realizadas pela metodologia de máxima verossimilhança considerando-se dois modelos: o modelo 1 incluiu como covariáveis os efeitos aditivos diretos e maternos, e os não aditivos das heterozigoses para os efeitos diretos e para o materno total, e o modelo 2 considerou também o efeito direto de epistasia. Para comparação dos modelos, foram utilizados o critério de informação de Akaike (AIC) e o critério de informação Bayesiano de Schwartz (BIC), e o teste de razão de verossimilhança. A inclusão da epistasia no modelo de avaliação genética pouco alterou as estimativas de componentes de (co)variâncias genéticas aditivas e, consequentemente, as herdabilidades. O teste de verossimilhança e o critério de Akaike sugeriram que o modelo 2, que inclui a epistasia, apresentou maior aderência aos dados para todas as características analisadas. O critério BIC indicou este modelo como o melhor apenas para P205. Para análise genética dessa população, o modelo que considerou o efeito de epistasia foi o mais adequado.

Quantitative Trait Loci Mapping and The Genetic Basis of Heterosis in Maize and Rice

Despite its importance to agriculture, the genetic basis of heterosis is still not well understood. The main competing hypotheses include dominance, overdominance, and epistasis. NC design III is an experimental design that. has been used for estimating the average degree of dominance of quantitative trait 106 (QTL) and also for studying heterosis. In this study, we first develop a multiple-interval mapping (MIM) model for design III that provides a platform to estimate the number, genomic positions, augmented additive and dominance effects, and epistatic interactions of QTL. The model can be used for parents with any generation of selling. We apply the method to two data sets, one for maize and one for rice. Our results show that heterosis in maize is mainly due to dominant gene action, although overdominance of individual QTL could not completely be ruled out due to the mapping resolution and limitations of NC design III. For rice, the estimated QTL dominant effects could not explain the observed heterosis. There is evidence that additive X additive epistatic effects of QTL could be the main cause for the heterosis in rice. The difference in the genetic basis of heterosis seems to be related to open or self pollination of the two species. The MIM model for NC design III is implemented in Windows QTL Cartographer, a freely distributed software.

The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

The protein phosphatase 7 regulates phytochrome signaling in Arabidopsis.

The psi2 mutant of Arabidopsis displays amplification of the responses controlled by the red/far red light photoreceptors phytochrome A (phyA) and phytochrome B (phyB) but no apparent defect in blue light perception. We found that loss-of-function alleles of the protein phosphatase 7 (AtPP7) are responsible for the light hypersensitivity in psi2 demonstrating that AtPP7 controls the levels of phytochrome signaling. Plants expressing reduced levels of AtPP7 mRNA display reduced blue-light induced cryptochrome signaling but no noticeable deficiency in phytochrome signaling. Our genetic analysis suggests that phytochrome signaling is enhanced in the AtPP7 loss of function alleles, including in blue light, which masks the reduced cryptochrome signaling. AtPP7 has been found to interact both in yeast and in planta assays with nucleotide-diphosphate kinase 2 (NDPK2), a positive regulator of phytochrome signals. Analysis of ndpk2-psi2 double mutants suggests that NDPK2 plays a critical role in the AtPP7 regulation of the phytochrome pathway and identifies NDPK2 as an upstream element involved in the modulation of the salicylic acid (SA)-dependent defense pathway by light. Thus, cryptochrome- and phytochrome-specific light signals synchronously control their relative contribution to the regulation of plant development. Interestingly, PP7 and NDPK are also components of animal light signaling systems.

PKS2: A link between phototropin signalling and auxin transport - a study on how plants sense and respond to light.

The blue light photoreceptors phototropins (phot1 and phot2 in Arabidopsis thaliana (L.)) carry out various light responses of great adaptive value that optimize plant growth. These processes include phototropism (the bending of an organ induced by unequal light distribution), chloroplast movements, stomatal opening, leaf flattening and solar tracking. The biochemical pathways controlling these important blue light responses are just starting to be elucidated. The PHYTOCHROME KINASE SUBSTRATE (PKS1-4) proteins - the subject of this research - have recently been identified as novel phototropism signalling components. PKS1 (the founding member of this family) interacts in a same complex in vivo with phot1 and the important phot1 signalling element NON-PHOTOTROPIC HYPOCOTYL 3 (NPH3). This suggested that the PKS may act as early components of phot signalling. This work further investigates the role of this protein family during phototropin signalling Genetic experiments clearly showed that the PKS do not control chloroplast movements or stomatal opening. However, PKS2 plays a critical role with NPH3 during leaf flattening and solar tracking. Epistasis data indicated that both proteins act in phot1 and phot2 pathways, which is consistent with their in vivo interaction with both phototropins. Because phototropism, leaf flattening and solar tracking are developmental processes regulated by the hormone auxin, the role of PKS2 and NPH3 during auxin homeostasis was also investigated. Interestingly, PKS2 loss-of-function restores leaf flattening in the auxin transporter mutant aux1. Moreover, PKS2 and NPH3 are found in a same complex with AUX1 in vivo. Taken together, these results suggest that PKS2 may act with NPH3 as a connecting point between phot signalling and auxin transport. Further experiments were performed to explore the molecular mode of action of PKS2 and NPH3 in this process. The significance of these results is discussed.

Reduced Prostasin (CAP1/PRSS8) Activity Eliminates HAI-1 and HAI-2 Deficiency-Associated Developmental Defects by Preventing Matriptase Activation.

Loss of either hepatocyte growth factor activator inhibitor (HAI)-1 or -2 is associated with embryonic lethality in mice, which can be rescued by the simultaneous inactivation of the membrane-anchored serine protease, matriptase, thereby demonstrating that a matriptase-dependent proteolytic pathway is a critical developmental target for both protease inhibitors. Here, we performed a genetic epistasis analysis to identify additional components of this pathway by generating mice with combined deficiency in either HAI-1 or HAI-2, along with genes encoding developmentally co-expressed candidate matriptase targets, and screening for the rescue of embryonic development. Hypomorphic mutations in Prss8, encoding the GPI-anchored serine protease, prostasin (CAP1, PRSS8), restored placentation and normal development of HAI-1-deficient embryos and prevented early embryonic lethality, mid-gestation lethality due to placental labyrinth failure, and neural tube defects in HAI-2-deficient embryos. Inactivation of genes encoding c-Met, protease-activated receptor-2 (PAR-2), or the epithelial sodium channel (ENaC) alpha subunit all failed to rescue embryonic lethality, suggesting that deregulated matriptase-prostasin activity causes developmental failure independent of aberrant c-Met and PAR-2 signaling or impaired epithelial sodium transport. Furthermore, phenotypic analysis of PAR-1 and matriptase double-deficient embryos suggests that the protease may not be critical for focal proteolytic activation of PAR-2 during neural tube closure. Paradoxically, although matriptase auto-activates and is a well-established upstream epidermal activator of prostasin, biochemical analysis of matriptase- and prostasin-deficient placental tissues revealed a requirement of prostasin for conversion of the matriptase zymogen to active matriptase, whereas prostasin zymogen activation was matriptase-independent.

Variation in the intensity of inbreeding depression among successive life-cycle stages and generations in gynodioecious Silene vulgaris (Caryophyllaceae).

Inbreeding depression is one of the hypotheses explaining the maintenance of females within gynodioecious plant populations. However, the measurement of fitness components in selfed and outcrossed progeny depends on life-cycle stage and the history of inbreeding. Comparative data indicate that strong inbreeding depression is more likely to occur at later life-cycle stages. We used hermaphrodite individuals of Silene vulgaris originating from three populations located in different valleys in the Swiss Alps to investigate the effect of two generations of self- and cross-fertilization on fitness components among successive stages of the life cycle in a glasshouse experiment. We detected significant inbreeding depression for most life-cycle stages including: the number of viable and aborted seeds per fruit, probability of germination, above ground biomass, probability of flowering, number of flowers per plant, flower size and pollen viability. Overall, the intensity of inbreeding depression increased among successive stages of the life cycle and cumulative inbreeding depression was significantly stronger in the first generation (delta approximately 0.5) compared with the second generation (delta approximately 0.35). We found no evidence for synergistic epistasis in our experiment. Our finding of more intense inbreeding depression during later stages of the life cycle may help to explain the maintenance of females in gynodioecious populations of S. vulgaris because purging of genetic load is less likely to occur.

The fou2 gain-of-function allele and the wild-type allele of Two Pore Channel 1 contribute to different extents or by different mechanisms to defense gene expression in Arabidopsis.

The fatty acid oxygenation up-regulated 2 (fou2) mutant in Arabidopsis thaliana creates a gain-of-function allele in a non-selective cation channel encoded by the Two Pore Channel 1 (TPC1) gene. This mutant genetically implicates cation fluxes in the control of the positive feedback loop whereby jasmonic acid (JA) stimulates its own synthesis. In this study we observed extensive transcriptome reprogramming in healthy fou2 leaves closely resembling that induced by treatment with methyl jasmonate, biotic stresses and the potassium starvation response. Proteomic analysis of fou2 leaves identified increased levels of seven biotic stress- and JA-inducible proteins. In agreement with these analyses, epistasis studies performed by crossing fou2 with aos indicated that elevated levels of JA in fou2 are the major determinant of the mutant phenotype. In addition, generation of fou2 aba1-5, fou2 etr1-1 and fou2 npr1-1 double mutants showed that the fou2 phenotype was only weakly affected by ABA levels and unaffected by mutations in NPR1 and ETR1. The results now suggest possible mechanisms whereby fou2 could induce JA synthesis/signaling early in the wound response. In contrast to fou2, transcriptome analysis of a loss-of-function allele of TPC1, tpc1-2, revealed no differential expression of JA biosynthesis genes in resting leaves. However, the analysis disclosed reduced mRNA levels of the pathogenesis-related genes PDF1.2a and THI2.1 in healthy and diseased tpc1-2 leaves. The results suggest that wild-type TPC1 contributes to their expression by mechanisms somewhat different from those affecting their expression in fou2.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.