Selection-mutation balance models with epistatic selection

We present an application of birth-and-death processes on configuration spaces to a generalized mutation4 selection balance model. The model describes the aging of population as a process of accumulation of mu5 tations in a genotype. A rigorous treatment demands that mutations correspond to points in abstract spaces. 6 Our model describes an infinite-population, infinite-sites model in continuum. The dynamical equation which 7 describes the system, is of Kimura-Maruyama type. The problem can be posed in terms of evolution of states 8 (differential equation) or, equivalently, represented in terms of Feynman-Kac formula. The questions of interest 9 are the existence of a solution, its asymptotic behavior, and properties of the limiting state. In the non-epistatic 10 case the problem was posed and solved in [Steinsaltz D., Evans S.N., Wachter K.W., Adv. Appl. Math., 2005, 11 35(1)]. In our model we consider a topological space X as the space of positions of mutations and the influence of epistatic potentials

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Detection of unusual mutation within the VP1 region of different re-isolates of poliovirus Sabin vaccine

In the present study, a genomic analysis of full VP1 sequence region of 15 clinical re-isolates (14 healthy vaccinees and one bone marrow tumor patient) was conducted, aiming to the identification of mutations and to the assessment of their impact on virus fitness, providing also insights relevant with the natural evolution of Sabin strains. Clinical re-isolates were analyzed by RT-PCR, sequencing and computational analysis. Some re-isolates were characterized by an unusual mutational pattern in which non-synonymous mutations outnumbered the synonymous ones. Furthermore, the majority of amino-acid substitutions were located in the capsid exterior, specifically in N-Ags, near N-Ags and in the north rim of the canyon. Also mutations, which are well-known determinants of attenuation, were identified. The results of this study propose that some re-isolates are characterized by an evolutionary pattern in which non-synonymous mutations with a direct phenotypic impact on viral fitness are fixed in viral genomes, in spite of synonymous ones with no phenotypic impact on viral fitness. Results of the present retrospective characterization of Sabin clinical re-isolates, based on the full VP1 sequence, suggest that vaccine-derived viruses may make their way through narrow breaches and may evolve into transmissible pathogens even in adequately immunized populations. For this reason increased poliovirus laboratory surveillance should be permanent and full VP1 sequence analysis should be conducted even in isolates originating from healthy vaccinees.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Dual role of the beta(2)-adrenergic receptor C terminus for the binding of beta-arrestin and receptor internalization

Homologous desensitization of beta(2)-adrenergic and other G-protein-coupled receptors is a two-step process. After phosphorylation of agonist-occupied receptors by G-protein-coupled receptor kinases, they bind beta-arrestins, which triggers desensitization and internalization of the receptors. Because it is not known which regions of the receptor are recognized by beta-arrestins, we have investigated beta-arrestin interaction and internalization of a set of mutants of the human beta(2)-adrenergic receptor. Mutation of the four serine/threonine residues between residues 355 and 364 led to the loss of agonist-induced receptor-beta-arrestin2 interaction as revealed by fluorescence resonance energy transfer (FRET), translocation of beta-arrestin2 to the plasma membrane, and receptor internalization. Mutation of all seven serine/threonine residues distal to residue 381 did not affect agonist-induced receptor internalization and beta-arrestin2 translocation. A beta(2)-adrenergic receptor truncated distal to residue 381 interacted normally with beta-arrestin2, whereas its ability to internalize in an agonist-dependent manner was compromised. A similar impairment of internalization was observed when only the last eight residues of the C terminus were deleted. Our experiments show that the C terminus distal to residue 381 does not affect the initial interaction between receptor and beta-arrestin, but its last eight amino acids facilitate receptor internalization in concert with beta-arrestin2.

An intronic growth hormone receptor mutation causing activation of a pseudoexon is associated with a broad spectrum of growth hormone insensitivity phenotypes

Context: Inherited GH insensitivity (GHI) is usually caused by mutations in the GH receptor (GHR). Patients present with short stature associated with high GH and low IGF-I levels and may have midfacial hypoplasia ( typical Laron syndrome facial features). We previously described four mildly affected GHI patients with an intronic mutation in the GHR gene (A.(1) -> G.(1) substitution in intron 6), resulting in the activation of a pseudoexon (6 Psi) and inclusion of 36 amino acids. Objective: The study aimed to analyze the clinical and genetic characteristics of additional GHI patients with the pseudoexon (6 Psi) mutation. Design/Patients: Auxological, biochemical, genetic, and haplotype data from seven patients with severe short stature and biochemical evidence of GHI were assessed. Main Outcome Measures: We assessed genotype-phenotype relationship. Results: One patient belongs to the same extended family, previously reported. She has normal facial features, and her IGF-I levels are in the low-normal range for age. The six unrelated patients, four of whom have typical Laron syndrome facial features, have heights ranging from -3.3 to -6.0 SD and IGF-I levels that vary from normal to undetectable. We hypothesize that the marked difference in biochemical and clinical phenotypes might be caused by variations in the splicing efficiency of the pseudoexon. Conclusions: Activation of the pseudoexon in the GHR gene can lead to a variety of GHI phenotypes. Therefore, screening for the presence of this mutation should be performed in all GHI patients without mutations in the coding exons.

Differential recognition of Staphylococcus aureus quorum-sensing signals depends on both extracellular loops 1 and 2 of the transmembrane sensor AgrC.

Virulence in Staphylococcus aureus is regulated via agr-dependent quorum sensing in which an autoinducing peptide (AIP) activates AgrC, a histidine protein kinase. AIPs are usually thiolactones containing seven to nine amino acid residues in which the thiol of the central cysteine is linked to the alpha-carboxyl of the C-terminal amino acid residue. The staphylococcal agr locus has diverged such that the AIPs of the four different S. aureus agr groups self-activate but cross-inhibit. Consequently, although the agr system is conserved among the staphylococci, it has undergone significant evolutionary divergence whereby to retain functionality, any changes in the AIP-encoding gene (agrD) that modifies AIP structure must be accompanied by corresponding changes in the AgrC receptor. Since AIP-1 and AIP-4 only differ by a single amino acid, we compared the transmembrane topology of AgrC1 and AgrC4 to identify amino acid residues involved in AIP recognition. As only two of the three predicted extracellular loops exhibited amino acid differences, site-specific mutagenesis was used to exchange the key AgrC1 and AgrC4 amino acid residues in each loop either singly or in combination. A novel lux-based agrP3 reporter gene fusion was constructed to evaluate the response of the mutated AgrC receptors. The data obtained revealed that while differential recognition of AIP-1 and AIP-4 depends primarily on three amino acid residues in loop 2, loop 1 is essential for receptor activation by the cognate AIP. Furthermore, a single mutation in the AgrC1 loop 2 resulted in conversion of (Ala5)AIP-1 from a potent antagonist to an activator, essentially resulting in the forced evolution of a new AIP group. Taken together, our data indicate that loop 2 constitutes the predicted hydrophobic pocket that binds the AIP thiolactone ring while the exocyclic amino acid tail interacts with loop 1 to facilitate receptor activation.

Regulation of expression of the rat orthologue of mouse double minute 2 (MDM2) by H2O2-induced oxidative stress in neonatal rat cardiac myocytes.

The Mdm2 ubiquitin ligase is an important regulator of p53 abundance and p53-dependent apoptosis. Mdm2 expression is frequently regulated by a p53 Mdm2 autoregulatory loop whereby p53 stimulates Mdm2 expression and hence its own degradation. Although extensively studied in cell lines, relatively little is known about Mdm2 expression in heart where oxidative stress (exacerbated during ischemia-reperfusion) is an important pro-apoptotic stimulus. We demonstrate that Mdm2 transcript and protein expression are induced by oxidative stress (0.2 mm H(2)O(2)) in neonatal rat cardiac myocytes. In other cells, constitutive Mdm2 expression is regulated by the P1 promoter (5' to exon 1), with inducible expression regulated by the P2 promoter (in intron 1). In myocytes, H(2)O(2) increased Mdm2 expression from the P2 promoter, which contains two p53-response elements (REs), one AP-1 RE, and two Ets REs. H(2)O(2) did not detectably increase expression of p53 mRNA or protein but did increase expression of several AP-1 transcription factors. H(2)O(2) increased binding of AP-1 proteins (c-Jun, JunB, JunD, c-Fos, FosB, and Fra-1) to an Mdm2 AP-1 oligodeoxynucleotide probe, and chromatin immunoprecipitation assays showed it increased binding of c-Jun or JunB to the P2 AP-1 RE. Finally, antisense oligonucleotide-mediated reduction of H(2)O(2)-induced Mdm2 expression increased caspase 3 activation. Thus, increased Mdm2 expression is associated with transactivation at the P2 AP-1 RE (rather than the p53 or Ets REs), and Mdm2 induction potentially represents a cardioprotective response to oxidative stress.

Inhibition of synaptic transmission and G protein modulation by synthetic CaV2.2 Ca2+ channel peptides

Abstract: Modulation of presynaptic voltage-dependent Ca+ channels is a major means of controlling neurotransmitter release. The CaV 2.2 Ca2+ channel subunit contains several inhibitory interaction sites for Gβγ subunits, including the amino terminal (NT) and I–II loop. The NT and I–II loop have also been proposed to undergo a G protein-gated inhibitory interaction, whilst the NT itself has also been proposed to suppress CaV 2 channel activity. Here, we investigate the effects of an amino terminal (CaV 2.2[45–55]) ‘NT peptide’ and a I–II loop alpha interaction domain (CaV 2.2[377–393]) ‘AID peptide’ on synaptic transmission, Ca2+ channel activity and G protein modulation in superior cervical ganglion neurones (SCGNs). Presynaptic injection of NT or AID peptide into SCGN synapses inhibited synaptic transmission and also attenuated noradrenaline-induced G protein modulation. In isolated SCGNs, NT and AID peptides reduced whole-cell Ca2+ current amplitude, modified voltage dependence of Ca2+ channel activation and attenuated noradrenaline-induced G protein modulation. Co-application of NT and AID peptide negated inhibitory actions. Together, these data favour direct peptide interaction with presynaptic Ca2+ channels, with effects on current amplitude and gating representing likely mechanisms responsible for inhibition of synaptic transmission. Mutations to residues reported as determinants of Ca2+ channel function within the NT peptide negated inhibitory effects on synaptic transmission, Ca2+ current amplitude and gating and G protein modulation. A mutation within the proposed QXXER motif for G protein modulation did not abolish inhibitory effects of the AID peptide. This study suggests that the CaV 2.2 amino terminal and I–II loop contribute molecular determinants for Ca2+ channel function; the data favour a direct interaction of peptides with Ca2+ channels to inhibit synaptic transmission and attenuate G protein modulation. Non-technical summary: Nerve cells (neurones) in the body communicate with each other by releasing chemicals (neurotransmitters) which act on proteins called receptors. An important group of receptors (called G protein coupled receptors, GPCRs) regulate the release of neurotransmitters by an action on the ion channels that let calcium into the cell. Here, we show for the first time that small peptides based on specific regions of calcium ion channels involved in GPCR signalling can themselves inhibit nerve cell communication. We show that these peptides act directly on calcium channels to make them more difficult to open and thus reduce calcium influx into native neurones. These peptides also reduce GPCR-mediated signalling. This work is important in increasing our knowledge about modulation of the calcium ion channel protein; such knowledge may help in the development of drugs to prevent signalling in pathways such as those involved in pain perception.

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Use of a LightCycler gyrA mutation assay for rapid identification of mutations conferring decreased susceptibility to ciprofloxacin in multiresistant Salmonella enterica serotype typhimurium DT104 isolates

A LightCycler-based PCR-hybridization gyrA mutation assay (GAMA) was developed to rapidly detect gyrA point mutations in multiresistant (MR) Salmonella enterica serotype Typhimurium DT104 with decreased susceptibility to ciprofloxacin (MIC, 0.25 to 1.0 mg/liter). Ninety-two isolates (49 human, 43 animal) were tested with three individual oligonucleotide probes directed against an Asp-87-to-Asn (GAC --> AAC) mutation, an Asp-87-to-Gly (GAC --> GGC) mutation, and a Ser-83-to-Phe (TCC --> TTC) mutation. Strains homologous to the probes could be distinguished from strains that had different mutations by their probe-target melting temperatures. Thirty-seven human and 30 animal isolates had an Asp-87-to-Asn substitution, 6 human and 6 animal isolates had a Ser-83-to-Phe substitution, and 5 human and 2 animal isolates had an Asp-87-to-Gly substitution. The remaining six strains all had mismatches with the three probes and therefore different gyrA mutations. The sequencing of gyrA from these six isolates showed that one human strain and two animal strains had an Asp-87-to-Tyr (GAC --> TAC) substitution and two animal strains had a Ser-83-to-Tyr (TCC --> TAC) substitution. One animal strain had no gyrA mutation, suggesting that this isolate had a different mechanism of resistance. Fifty-eight of the strains tested were indistinguishable by several different typing methods including antibiograms, pulsed-field gel gel electrophoresis, and plasmid profiling, although they could be further subdivided according to gyrA mutation. This study confirmed that MR DT104 with decreased susceptibility to ciprofloxacin from humans and food animals in England and Wales may have arisen independently against a background of clonal spread of MR DT104.

c-Cbl mediates ubiquitination, degradation, and down-regulation of human protease-activated receptor 2

Mechanisms that arrest G-protein-coupled receptor (GPCR) signaling prevent uncontrolled stimulation that could cause disease. Although uncoupling from heterotrimeric G-proteins, which transiently arrests signaling, is well described, little is known about the mechanisms that permanently arrest signaling. Here we reported on the mechanisms that terminate signaling by protease-activated receptor 2 (PAR(2)), which mediated the proinflammatory and nociceptive actions of proteases. Given its irreversible mechanism of proteolytic activation, PAR(2) is a model to study the permanent arrest of GPCR signaling. By immunoprecipitation and immunoblotting, we observed that activated PAR(2) was mono-ubiquitinated. Immunofluorescence indicated that activated PAR(2) translocated from the plasma membrane to early endosomes and lysosomes where it was degraded, as determined by immunoblotting. Mutant PAR(2) lacking intracellular lysine residues (PAR(2)Delta14K/R) was expressed at the plasma membrane and signaled normally but was not ubiquitinated. Activated PAR(2) Delta14K/R internalized but was retained in early endosomes and avoided lysosomal degradation. Activation of wild type PAR(2) stimulated tyrosine phosphorylation of the ubiquitin-protein isopeptide ligase c-Cbl and promoted its interaction with PAR(2) at the plasma membrane and in endosomes in an Src-dependent manner. Dominant negative c-Cbl lacking the ring finger domain inhibited PAR(2) ubiquitination and induced retention in early endosomes, thereby impeding lysosomal degradation. Although wild type PAR(2) was degraded, and recovery of agonist responses required synthesis of new receptors, lysine mutation and dominant negative c-Cbl impeded receptor ubiquitination and degradation and allowed PAR(2) to recycle and continue to signal. Thus, c-Cbl mediated ubiquitination and lysosomal degradation of PAR(2) to irrevocably terminate signaling by this and perhaps other GPCRs.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Protease-activated receptor 2 (PAR2) protein and transient receptor potential vanilloid 4 (TRPV4) protein coupling is required for sustained inflammatory signaling

G protein-coupled receptors of nociceptive neurons can sensitize transient receptor potential (TRP) ion channels, which amplify neurogenic inflammation and pain. Protease-activated receptor 2 (PAR(2)), a receptor for inflammatory proteases, is a major mediator of neurogenic inflammation and pain. We investigated the signaling mechanisms by which PAR(2) regulates TRPV4 and determined the importance of tyrosine phosphorylation in this process. Human TRPV4 was expressed in HEK293 cells under control of a tetracycline-inducible promoter, allowing controlled and graded channel expression. In cells lacking TRPV4, the PAR(2) agonist stimulated a transient increase in [Ca(2+)](i). TRPV4 expression led to a markedly sustained increase in [Ca(2+)](i). Removal of extracellular Ca(2+) and treatment with the TRPV4 antagonists Ruthenium Red or HC067047 prevented the sustained response. Inhibitors of phospholipase A(2) and cytochrome P450 epoxygenase attenuated the sustained response, suggesting that PAR(2) generates arachidonic acid-derived lipid mediators, such as 5',6'-EET, that activate TRPV4. Src inhibitor 1 suppressed PAR(2)-induced activation of TRPV4, indicating the importance of tyrosine phosphorylation. The TRPV4 tyrosine mutants Y110F, Y805F, and Y110F/Y805F were expressed normally at the cell surface. However, PAR(2) was unable to activate TRPV4 with the Y110F mutation. TRPV4 antagonism suppressed PAR(2) signaling to primary nociceptive neurons, and TRPV4 deletion attenuated PAR(2)-stimulated neurogenic inflammation. Thus, PAR(2) activation generates a signal that induces sustained activation of TRPV4, which requires a key tyrosine residue (TRPV4-Tyr-110). This mechanism partly mediates the proinflammatory actions of PAR(2).

In silico identification and three-dimensional modelling of the missense mutation in ADAMTS2 in a sheep flock with dermatosparaxis

Background Dermatosparaxis (Ehlers–Danlos syndrome in humans) is characterized by extreme fragility of the skin. It is due to the lack of mature collagen caused by a failure in the enzymatic processing of procollagen I. We investigated the condition in a commercial sheep flock. Hypothesis/Objectives Mutations in the ADAM metallopeptidase with thrombospondin type 1 motif, 2 (ADAMTS2) locus, are involved in the development of dermatosparaxis in humans, cattle and the dorper sheep breed; consequently, this locus was investigated in the flock. Animals A single affected lamb, its dam, the dam of a second affected lamb and the rams in the flock were studied. Methods DNA was purified from blood, PCR primers were used to detect parts of the ADAMS2 gene and nucleotide sequencing was performed using Sanger's procedure. Skin samples were examined using standard histology procedures. Results A missense mutation was identified in the catalytic domain of ADAMTS2. The mutation is predicted to cause the substitution in the mature ADAMTS2 of a valine molecule by a methionine molecule (V15M) affecting the catalytic domain of the enzyme. Both the ‘sorting intolerant from tolerant’ (SIFT) and the PolyPhen-2 methodologies predicted a damaging effect for the mutation. Three-dimensional modelling suggested that this mutation may alter the stability of the protein folding or distort the structure, causing the protein to malfunction. Conclusions and clinical importance Detection of the mutation responsible for the pathology allowed us to remove the heterozygote ram, thus preventing additional cases in the flock.