A gene's eye view of epistasis,selection and speciation

In this mini‐review, I discuss the effects of gene interaction or epistasis from a `gene's eye view.' By a `gene's eye view' of epistasis, I mean that I will consider a single, bi‐allelic locus, A , whose effects on fitness result only from its interactions with alleles of another, unknown locus, X . I will show how changes in the frequencies of alleles at the background locus affect the relationship of alleles at the A ‐locus to fitness. Changing the genetic background changes the fundamental characteristics of the A ‐locus, such as the magnitude and sign of allelic effects on fitness, and, consequently, it changes the strength and pattern of selection. I consider each of the four kinds of pair–wise interactions between two loci and show that some kinds of epistasis are more sensitive than others to population genetic subdivision. Lastly, I show that some kinds of epistasis are more likely than others to affect the process of speciation and contribute to or be responsible for general genetic features of interspecific hybrids, such as Haldane's rule.     

Effect on linkage disequilibrium of selection for a quantitative character with epistasis

Summary Selection for a character controlled by additive genes induces linkage disequilibrium which reduces the additive genetic variance usable for further selective gains. Additive x additive epistasis contributes to selection response through development of linkage disequilibrium between interacting loci. To investigate the relative importance of the two effects of linkage disequilibrium, formulae are presented and results are reported of simulations using models involving additive, additive x additive and dominance components. The results suggest that so long as epistatic effects are not large relative to additive effects, and the proportion of pairs of loci which show epistasis is not very high, the predominant effect of linkage disequilibrium will be to reduce the rate of selection response.     

Comparison between marker-assisted selection and phenotypical selection in a set of Arabidopsis thaliana recombinant inbred lines

 Parents were selected from a well-characterised Arabidopsis recombinant inbred line (RIL) population based on (1) their phenotype for flowering time or (2) marker and QTL information that had been assessed previously. The F₂ offspring obtained from pairs of selected RILs was analysed for these traits, and the results obtained with these two methods of selection were compared. Selection based on marker and QTL information gave approximately the same result as selection based on phenotype. The relative high heritability of flowering time in Arabidopsis facilitated successful phenotypical selection. The difference in selection result that was anticipated to be in favour of the marker-assisted approach was therefore not observed. Received: 29 November 1997 / Accepted: 8 June 1998     

The effect of selection for phenotypical characters on the chromosome constitution in spring triticale

A set of spring triticale selection forms was studied with the use of microsatellite markers specific for D-genome chromosomes and genomic in situ hybridization. 2R/2D substitutions were found in six of twelve forms. Two lines displayed segregation for the substitution: plants with the substitution and without it, and plants heterozygous for the substitution were revealed. Segregation of the lines and the presence in them of telocentric chromosomes of wheat and rye in the monosomic state suggest that the process of formation is still in progress. The data obtained also demonstrate that triticale is under selection for the 2R/2D substitution as controlling the spectrum of adaptively and selectively valuable traits.     

Antigen-independent selection of intracellular stable antibody frameworks

The intracellular expression of highly specific antibody fragments (intrabodies) in eukaryotes has a great potential in functional genomics and therapeutics. However, since the intracellular reducing environment prevents formation of the conserved intrachain disulfide bonds, most antibodies do not fold properly and are therefore inactive inside cells. The few antibodies that have been found to function in an intracellular environment and that have been characterized for their biophysical properties have generally shown a high degree of stability and solubility. Thus, for intracellular expression and application, very stable antibody frameworks are needed that can correctly fold even in the absence of disulfide bonds and that do not aggregate. Here, we present and discuss a novel method, named Quality Control, which allows selection of stable and soluble antibody frameworks in vivo without the requirement or knowledge of antigens. This system is based on the expression of single-chain antibody fragments (scFvs) fused to a selectable marker that can control gene expression and cell growth. The activity of such a selectable marker fused to various scFvs that have been biophysically characterized correlated with the solubility and stability of the scFv moieties. This antigen-independent intrabody selection system was applied to screen scFv libraries for identifying stable and soluble frameworks, which subsequently served as acceptor backbones to construct intrabody libraries by randomization of hypervariable loops.     

On the evolution of epistasis I: diploids under selection

One interpretation of recent literature on the evolution of phenotypic modularity is that evolution should act to decrease the degree of interaction between genes that contribute to different phenotypes. This issue is addressed directly here using a fitness scheme determined by two genetic loci and a third locus which modifies a measure of statistical interaction between the fitnesses due to the first two. The equilibrium structure of such an epistasis-modifying locus is studied. It is shown that under well-specified conditions a modifying allele that increases epistasis succeeds. In other words, genetic interactions tend to become stronger. It is speculated that this occurs because the mean fitness in such models is locally increasing as a function of the degree of epistasis.     

Convergence of multilocus systems under weak epistasis or weak selection

 The convergence of multilocus systems under viability selection with constant fitnesses is investigated. Generations are discrete and nonoverlapping; the monoecious population mates at random. The number of multiallelic loci, the linkage map, dominance, and epistasis are arbitrary. It is proved that if epistasis or selection is sufficiently weak (and satisfies a certain nondegeneracy assumption whose genericity we establish), then there is always convergence to some equilibrium point. In particular, cycling cannot occur. The behavior of the mean fitness and some other aspects of the dynamics are also analyzed. Received: 15 November 1997 / Revised version: 25 May 1998     

Physiological epistasis, ontogenetic conflict and natural selection on physiology and life history

Ontogenetic conflict arises when optima for alleles governingfitness variation differ between juveniles and adults or betweenadult sexes. Loci that govern development of alternative phenotypesin the sexes, hereafter termed morph-determining loci, mediatedevelopment through the endocrine system. Morphotypic selectionis defined to be multivariate selection favoring discrete alternativemorphotypes (e.g., optima). When the optimal combinations ofalleles for alternative morphs differ between the sexes, itgenerates conflicting selection pressure and thus ontogeneticconflict. Selection on morph alleles promotes ontogenetic conflictbecause it perturbs physiological epistasis that governs theexpression of male versus female traits. Expression of physiologicaltraits arises from homeostasis that maintains trait expressionwithin a normal range. The genetic basis of homeostasis is likelyto arise from interactions among several genes (e.g., geneticepistasis) or protein products (e.g., physiological epistasis).For example, endocrine regulation arises from interactions betweengondatropins, which are protein hormones produced by the hypothalamic-pituitaryglands, and steroid hormones, which are produced by the gonads(e.g., HPG axis). The side-blotched lizard system is discussedwith respect to physiological bases of ontogenetic conflict.We also describe a novel molecular marker strategy for uncoveringgenome-wide physiological epistasis in nature. Finally, ontogeneticconflict exerts selection on females to evolve mate selectionor cryptic choice that is reflected in different sires beingchosen for son versus daughter production. We describe how side-blotchedlizard females ameliorate ontogenetic conflict by cryptic choiceof male genotypes to produce sons versus daughters.     

Two-locus epistasis with sexually antagonistic selection: a genetic Parrondo's paradox

An example is provided where, with antagonistic selection and epistatic interaction of alleles at two loci, an autosomal allele can rise in frequency, persist in the population, and even continue to fixation, despite having an apparently lower average fitness than the alternative allele, in a process similar to Parrondo's paradox.     

Synergistic epistasis and alternative hypotheses

Inbreeding generally results in deleterious shifts in mean fitness. If the fitness response to increasing inbreeding coefficient is non-linear, this suggests a contribution of epistasis to inbreeding depression. In a cross-breeding experiment, Salathe & Ebert (2003. J. Evol. Biol. 16: 976-985) tested and found the presence of this non-linearity in Daphnia magna. They argue that epistatic interactions cause this non-linearity. We argue here that their experimental protocol does not allow disentangling the effect of synergistic epistasis from two alternative hypotheses, namely hybrid vigour and statistical non-independence of data.     

Naturally segregating loci exhibit epistasis for fitness

The extent to which gene interaction or epistasis contributes to fitness variation within populations remains poorly understood, despite its importance to a myriad of evolutionary questions. Here, we report a multi-year field study estimating fitness of Mimulus guttatus genetic lines in which pairs of naturally segregating loci exist in an otherwise uniform background. An allele at QTL x5b—a locus originally mapped for its effect on flower size—positively affects survival if combined with one genotype at quantitative trait locus x10a (aa) but has negative effects when combined with the other genotypes (Aa and AA). The viability differences between genotypes parallel phenotypic differences for the time and node at which a plant flowers. Viability is negatively correlated with fecundity across genotypes, indicating antagonistic pleiotropy for fitness components. This trade-off reduces the genetic variance for total fitness relative to the individual fitness components and thus may serve to maintain variation. Additionally, we find that the effects of each locus and their interaction often vary with the environment.     

A fast and powerful W-test for pairwise epistasis testing

Epistasis plays an essential role in the development of complex diseases. Interaction methods face common challenge of seeking a balance between persistent power, model complexity, computation efficiency, and validity of identified bio-markers. We introduce a novel W-test to identify pairwise epistasis effect, which measures the distributional difference between cases and controls through a combined log odds ratio. The test is model-free, fast, and inherits a Chi-squared distribution with data adaptive degrees of freedom. No permutation is needed to obtain the P-values. Simulation studies demonstrated that the W-test is more powerful in low frequency variants environment than alternative methods, which are the Chi-squared test, logistic regression and multifactor-dimensionality reduction (MDR). In two independent real bipolar disorder genome-wide associations (GWAS) datasets, the W-test identified significant interactions pairs that can be replicated, including SLIT3-CENPN, SLIT3-TMEM132D, CNTNAP2-NDST4 and CNTCAP2-RTN4R. The genes in the pairs play central roles in neurotransmission and synapse formation. A majority of the identified loci are undiscoverable by main effect and are low frequency variants. The proposed method offers a powerful alternative tool for mapping the genetic puzzle underlying complex disorders.     

Mathematical model of clonal selection and antibody production. II

Several modifications are proposed to a recent mathematical model (Bell, 1970) of the clonal selection and antibody production which take place when an adult animal is injected with an antigen. In the original model, antigen molecules were assumed univalent, i.e. to have only one combining site per molecule, and the first modification is an allowance for multivalent antigens by permitting an antigen molecule to interact with only one cell at a time. It is found that, for antigens with few (? 10) sites per molecule this modification is not important while for many sites per molecule, the modification will reduce the response as compared to that from the same number of independent antigenic sites. In section 3, it is seen that by restricting the number of cells which can arise in an immune response, much more realistic responses to high antigen doses are obtained. Moreover, the response is made more predictable by assuming that both target and proliferating cells are stimulated by antigen when a fraction of their receptors, between f_min and f_max, is bound to antigen. The parameter f_min is shown to determine the extent to which a molecule which can be recognized by antibodies will also serve as an immunogen giving rise to cellular proliferation and antibody production. In section 4, it is found that if more plasma cells than memory cells result from antigen stimulation, then weak stimulation will lead to a depletion of target plus memory cells and to at least partial immunological paralysis. Optimum antigen doses and times for the induction of such paralysis are examined. In section 5, precipitation of antigen-antibody mixtures is considered and it is shown that the number of doubly bound antibody molecules per antigen site determines whether precipitation will occur. This number is easily computed for heterogeneous bivalent antibodies.     

Re-evaluation of genetic formulae and epistasis

Experimental verification of some quantitative results, following from the new theory for investigation of non-allelic genes' interactions in the diallele analysis by Hayman was conducted. We postulated earlier a special kind of epistasis HC which is the new explication of epistasis j--the interaction between homozygote and heterozygote by Mather and Jinks. Genetic interpretation of the epistasis HC establishes some relations between epistasis and the model of ecological-genetical control of quantitative characters. We found the invariant characteristic of genetic control which is stable in a number of characters, environments, treatments, inspite of redetermination of genetic formulae under conditions of environmental variability.     

Morphological and phenotypical characterization of Bacillus sporothermodurans

AIMS: Enumeration of resistant bacteria in ultra-high temperature (UHT) treated milk; morphological characterization and phenotyping of resistant strains by traditional and nontraditional methods and their identification by molecular biology. METHODS AND RESULTS: Modified standard plate count agar (PCA) and modified brain-heart infusion (BHI) agar were used for colony counts. Physiological culture traits were determined as suggested by Bergey's Manual of Systematic Bacteriology or in modified J-broth or in modified BHI agar. Scanning electron microscope (SEM) was used for microscopic examination. Strain identification was carried out by polymerase chain reaction (PCR). A total of 125 (62.81% of 199) samples were positive and the bacterial load was higher than 10(5) CFU ml(-1) in 46 samples (28.80% of 125). The 16S rRNA sequence of bacterial cultures obtained from UHT-treated milk was similar to that of Bacillus sporothermodurans M215 type strain((T)) and different biotypes were found by analysis of colony appearance, cell morphology and physiological traits. CONCLUSIONS: Bacillus sporothermodurans was the predominant sporigenous micro-organisms in UHT milk. SIGNIFICANCE AND IMPACT OF THE STUDY: BHI agar is more suitable than PCA for quality control of milk after UHT treatment. Modified J-broth medium is useful to determine selected physiological traits of B. sporothermodurans. The strains characterized and identified as B. sporothermodurans were significantly different compared with the type strain.     

Equilibrium correlations in a model for multidimensional epistasis

We investigate a statistical model for multidimensional epistasis. The genotype is devided into subsequences, and within each subsequence mutations which occur in a prescribed order are beneficial. The bit-string model used to represent the genotype, may be cast in the form of a ferromagnetic Ising model with a staggered field. We describe the actual correlations between mutations at different sites, within an equilibrium population at a given tolerance, which we define to be the temperature of the statistical ensemble.     

Concordance between phenotypical features and PCR-REA for the identification of Malassezia spp

The genus Malassezia has been recently revised and nowadays includes 11 species that cannot always be differentiated from each other by physiological and morphological tests. This study was aimed to evaluate the correlation between a molecular method and conventional phenotypic features in the identification of Malassezia spp. To achieve this aim, 92 Argentinean clinical strains isolated between 2001 and 2005 were analyzed along with three reference strains (Malassezia furfur CBS 7019, Malassezia sympodialis CBS 7222 and Malassezia slooffiae CBS 7956). By using PCR and restriction enzyme analysis with three different DNA endonucleases (PCR-REA), the molecular method consistently identified all three reference strains and all 92 clinical isolates as follows: 63 M. sympodialis, 18 M. furfur, 10 Malassezia globosa and one Malassezia obtusa. Phenotypic studies undentified 85 clinical isolates and two of the reference strains (total agreement > 91%). In particular for M. sympodialis, M. furfur and M. globosa, the species more frequently involved in human pathology, the agreement ranged between 84 and 96%. This result suggests that phenotypic studies are suitable for the presumptive identification of important Malassezia species in the clinical medical mycology laboratories where molecular methodologies are not available.