Allele-specific gene expression in a wild nonhuman primate population

Natural populations hold enormous potential for evolutionary genetic studies, especially when phenotypic, genetic and environmental data are all available on the same individuals. However, untangling the genotype-phenotype relationship in natural populations remains a major challenge. Here, we describe results of an investigation of one class of phenotype, allele-specific gene expression (ASGE), in the well-studied natural population of baboons of the Amboseli basin, Kenya. ASGE measurements identify cases in which one allele of a gene is overexpressed relative to the alternative allele of the same gene, within individuals, thus providing a control for background genetic and environmental effects. Here, we characterize the incidence of ASGE in the Amboseli baboon population, focusing on the genetic and environmental contributions to ASGE in a set of eleven genes involved in immunity and defence. Within this set, we identify evidence for common ASGE in four genes. We also present examples of two relationships between cis-regulatory genetic variants and the ASGE phenotype. Finally, we identify one case in which this relationship is influenced by a novel gene-environment interaction. Specifically, the dominance rank of an individual's mother during its early life (an aspect of that individual's social environment) influences the expression of the gene CCL5 via an interaction with cis-regulatory genetic variation. These results illustrate how environmental and ecological data can be integrated into evolutionary genetic studies of functional variation in natural populations. They also highlight the potential importance of early life environmental variation in shaping the genetic architecture of complex traits in wild mammals.     

Structural insights into the role of the WW2 domain on tandem WW–PPxY motif interactions of oxidoreductase WWOX

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW–peptide interactions is not always intuitive. The WW domain–containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1–WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.     

Powerful testing via hierarchical linkage disequilibrium in haplotype association studies

Marginal tests based on individual SNPs are routinely used in genetic association studies. Studies have shown that haplotype‐based methods may provide more power in disease mapping than methods based on single markers when, for example, multiple disease‐susceptibility variants occur within the same gene. A limitation of haplotype‐based methods is that the number of parameters increases exponentially with the number of SNPs, inducing a commensurate increase in the degrees of freedom and weakening the power to detect associations. To address this limitation, we introduce a hierarchical linkage disequilibrium model for disease mapping, based on a reparametrization of the multinomial haplotype distribution, where every parameter corresponds to the cumulant of each possible subset of a set of loci. This hierarchy present in the parameters enables us to employ flexible testing strategies over a range of parameter sets: from standard single SNP analyses through the full haplotype distribution tests, reducing degrees of freedom and increasing the power to detect associations. We show via extensive simulations that our approach maintains the type I error at nominal level and has increased power under many realistic scenarios, as compared to single SNP and standard haplotype‐based studies. To evaluate the performance of our proposed methodology in real data, we analyze genome‐wide data from the Wellcome Trust Case‐Control Consortium.     

Role of the Cys 2-Cys 10 disulfide bond for the structure, stability, and folding kinetics of ribonuclease T1.

The Cys 2-Cys 10 disulfide bond in ribonuclease T1 was broken by substituting Cys 2 and Cys 10 by Ser and Asn, respectively, as present in ribonuclease F1. This C2S/C10N variant resembles the wild-type protein in structure and in catalytic activity. Minor structural changes were observed by 2-dimensional NMR in the local environment of the substituted amino acids only. The thermodynamic stability of ribonuclease T1 is strongly reduced by breaking the Cys 2-Cys 10 bond, and the free energy of denaturation is decreased by about 10 kJ/mol. The folding mechanism is not affected, and the trans to cis isomerizations of Pro 39 and Pro 55 are still the rate-limiting steps of the folding process. The differences in the time courses of unfolding and refolding are correlated with the decrease in stability: the folding kinetics of the wild-type protein and the C2S/C10N variant become indistinguishable when they are compared under conditions of identical stability. Apparently, the Cys 2-Cys 10 disulfide bond is important for the stability but not for the folding mechanism of ribonuclease T1. The breaking of this bond has the same effect on stability and folding kinetics as adding 1 M guanidinium chloride to the wild-type protein.     

Grandisol as an attractant for the sugar beet pest Bothynoderes affinis and data on other Lixinae species

As a result of field tests in Bulgaria and Hungary, cis‐2‐isopropenyl‐l‐methylcyclobutane ethanol (racemic grandisol) is reported for the first time as an attractant for Bothynoderes affinis (Schrank) (Coleoptera: Curculionidae, Lixinae), a member of the pest weevil complex of sugar beet. Dose–response experiments in the field using Csalomon TAL (modified pitfall) traps (Plant Protection Institute, CAR HAS, Budapest, Hungary) showed that catches of B. affinis adults increased with increasing attractant dose. In a subsequent experiment studying the effect of trap color (white, blue, yellow, fluorescent yellow, and transparent) all traps with the lure caught more than non‐baited control traps, and the highest number of adults was recorded in transparent and yellow baited traps. Trap color had a significant effect on the number of B. affinis females captured. Transparent TAL traps baited with 1–10 mg grandisol applied on rubber dispensers are recommended for the detection and monitoring of B. affinis. In addition to the target species, 17 other Lixinae species were captured during the field experiments, demonstrating for the first time the possible role of grandisol in the chemical communication systems of some of these species. A second locality of Lixus punctiventris Boheman (Lixinae, Lixini) in Bulgaria is reported. TAL traps baited with grandisol might be a useful tool for surveying Lixinae diversity in different biotopes.     

THE LOCI OF EVOLUTION: HOW PREDICTABLE IS GENETIC EVOLUTION?

Is genetic evolution predictable? Evolutionary developmental biologists have argued that, at least for morphological traits, the answer is a resounding yes. Most mutations causing morphological variation are expected to reside in the cis‐regulatory, rather than the coding, regions of developmental genes. This “cis‐regulatory hypothesis” has recently come under attack. In this review, we first describe and critique the arguments that have been proposed in support of the cis‐regulatory hypothesis. We then test the empirical support for the cis‐regulatory hypothesis with a comprehensive survey of mutations responsible for phenotypic evolution in multicellular organisms. Cis‐regulatory mutations currently represent approximately 22% of 331 identified genetic changes although the number of cis‐regulatory changes published annually is rapidly increasing. Above the species level, cis‐regulatory mutations altering morphology are more common than coding changes. Also, above the species level cis‐regulatory mutations predominate for genes not involved in terminal differentiation. These patterns imply that the simple question “Do coding or cis‐regulatory mutations cause more phenotypic evolution?” hides more interesting phenomena. Evolution in different kinds of populations and over different durations may result in selection of different kinds of mutations. Predicting the genetic basis of evolution requires a comprehensive synthesis of molecular developmental biology and population genetics.     

Brownian dynamics simulation of helix-capping motifs

Helix-capping motifs are believed to play an important role in stabilizing alpha-helices and defining helix start and stop signals. We performed microsecond scale Brownian dynamics simulations to study ten XAAD sequences, with X = (A,E,I,L,N,Q,S,T,V,Y), to examine their propensity to form helix capping motifs and correlate these results with those obtained from analyzing a structural database of proteins. For the widely studied capping box motif S**D, where the asterisk can be any amino acid residue, the simulations suggested that one of the two hydrogen bonds proposed earlier as a stabilizing factor might not be as important. On the other hand, side-chain interactions between the capping residue and the third residue downstream on the polypeptide chain might also play a role in stabilizing this motif. These results are consistent with explicit-solvent molecular dynamics simulations of two capping box motifs found in the proteins BPTI and alpha-dendrotoxin. Principal component analysis of the SAAD trajectory showed that the first three principal components, after those corresponding to translational-rotational motion were removed, accounted for more than half of the conformational fluctuations. The first component separated the conformational space into two parts with the all-helical conformation and the capping box motif lying largely in one part. The second component, on the other hand, could be used to describe conformational transitions between the all-helical form and the capping box motif.     

The C‐terminal calcium‐sensitive disordered motifs regulate isoform‐specific polymerization characteristics of calsequestrin

Calsequestrin (CASQ) exists as two distinct isoforms CASQ1 and CASQ2 in all vertebrates. Although the isoforms exhibit unique functional characteristic, the structural basis for the same is yet to be fully defined. Interestingly, the C‐terminal region of the two isoforms exhibit significant differences both in length and amino acid composition; forming Dn‐motif and DEXn‐motif in CASQ1 and CASQ2, respectively. Here, we investigated if the unique C‐terminal motifs possess Ca²⁺‐sensitivity and affect protein function. Sequence analysis shows that both the Dn‐ and DEXn‐motifs are intrinsically disordered regions (IDRs) of the protein, a feature that is conserved from fish to man. Using purified synthetic peptides, we show that these motifs undergo distinctive Ca²⁺‐mediated folding suggesting that these disordered motifs are Ca²⁺‐sensitivity. We generated chimeric proteins by swapping the C‐terminal portions between CASQ1 and CASQ2. Our studies show that the C‐terminal portions do not play significant role in protein folding. An interesting finding of the current study is that the switching of the C‐terminal portion completely reverses the polymerization kinetics. Collectively, these data suggest that these Ca²⁺‐sensitivity IDRs located at the back‐to‐back dimer interface influence isoform‐specific Ca²⁺‐dependent polymerization properties of CASQ. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 15–22, 2015.