Pseudogene rescue: an adaptive mechanism of codon reassignment

Alterations to the genetic code – codon reassignments – have occurred many times in life’s history, despite the fact that genomes are coadapted to their genetic codes and therefore alterations are likely to be maladaptive. A potential mechanism for adaptive codon reassignment, which could trigger either a temporary period of codon ambiguity or a permanent genetic code change, is the reactivation of a pseudogene by a nonsense suppressor mutant transfer RNA. I examine the population genetics of each stage of this process and find that pseudogene rescue is plausible and also readily explains some features of extant variability in genetic codes.

Stops making sense: translational trade-offs and stop codon reassignment

Background Efficient gene expression involves a trade-off between (i) premature termination of protein synthesis; and (ii) readthrough, where the ribosome fails to dissociate at the terminal stop. Sense codons that are similar in sequence to stop codons are more susceptible to nonsense mutation, and are also likely to be more susceptible to transcriptional or translational errors causing premature termination. We therefore expect this trade-off to be influenced by the number of stop codons in the genetic code. Although genetic codes are highly constrained, stop codon number appears to be their most volatile feature. Results In the human genome, codons readily mutable to stops are underrepresented in coding sequences. We construct a simple mathematical model based on the relative likelihoods of premature termination and readthrough. When readthrough occurs, the resultant protein has a tail of amino acid residues incorrectly added to the C-terminus. Our results depend strongly on the number of stop codons in the genetic code. When the code has more stop codons, premature termination is relatively more likely, particularly for longer genes. When the code has fewer stop codons, the length of the tail added by readthrough will, on average, be longer, and thus more deleterious. Comparative analysis of taxa with a range of stop codon numbers suggests that genomes whose code includes more stop codons have shorter coding sequences. Conclusions We suggest that the differing trade-offs presented by alternative genetic codes may result in differences in genome structure. More speculatively, multiple stop codons may mitigate readthrough, counteracting the disadvantage of a higher rate of nonsense mutation. This could help explain the puzzling overrepresentation of stop codons in the canonical genetic code and most variants.

Codon and mRNA sequence optimization of microdystrophin transgenes improves expression and physiological outcome in dystrophic mdx mice following AAV2/8 gene transfer

Duchenne muscular dystrophy is a fatal muscle-wasting disorder. Lack of dystrophin compromises the integrity of the sarcolemma and results in myofibers that are highly prone to contraction-induced injury. Recombinant adenoassociated virus (rAAV)-mediated dystrophin gene transfer strategies to muscle for the treatment of Duchenne muscular dystrophy (DMD) have been limited by the small cloning capacity of rAAV vectors and high titers necessary to achieve efficient systemic gene transfer. In this study, we assess the impact of codon optimization on microdystrophin (ΔAB/R3-R18/ΔCT) expression and function in the mdx mouse and compare the function of two different configurations of codon-optimized microdystrophin genes (ΔAB/R3-R18/ΔCT and ΔR4-R23/ΔCT) under the control of a muscle-restrictive promoter (Spc5-12). Codon optimization of microdystrophin significantly increases levels of microdystrophin mRNA and protein after intramuscular and systemic administration of plasmid DNA or rAAV2/8. Physiological assessment demonstrates that codon optimization of ΔAB/R3-R18/ΔCT results in significant improvement in specific force, but does not improve resistance to eccentric contractions compared with noncodon-optimized ΔAB/ R3-R18/ΔCT. However, codon-optimized microdystrophin ΔR4-R23/ΔCT completely restored specific force generation and provided substantial protection from contraction-induced injury. These results demonstrate that codon optimization of microdystrophin under the control of a muscle-specific promoter can significantly improve expression levels such that reduced titers of rAAV vectors will be required for efficient systemic administration.

Phylogeny, function, and evolution of the cupins, a structurally conserved, functionally diverse superfamily of proteins

The cupin superfamily is a group of functionally diverse proteins that are found in all three kingdoms of life, Archaea, Eubacteria, and Eukaryota. These proteins have a characteristic signature domain comprising two histidine- containing motifs separated by an intermotif region of variable length. This domain consists of six beta strands within a conserved beta barrel structure. Most cupins, such as microbial phosphomannose isomerases (PMIs), AraC- type transcriptional regulators, and cereal oxalate oxidases (OXOs), contain only a single domain, whereas others, such as seed storage proteins and oxalate decarboxylases (OXDCs), are bi-cupins with two pairs of motifs. Although some cupins have known functions and have been characterized at the biochemical level, the majority are known only from gene cloning or sequencing projects. In this study, phylogenetic analyses were conducted on the conserved domain to investigate the evolution and structure/function relationships of cupins, with an emphasis on single- domain plant germin-like proteins (GLPs). An unrooted phylogeny of cupins from a wide spectrum of evolutionary lineages identified three main clusters, microbial PMIs, OXDCs, and plant GLPs. The sister group to the plant GLPs in the global analysis was then used to root a phylogeny of all available plant GLPs. The resulting phylogeny contained three main clades, classifying the GLPs into distinct subfamilies. It is suggested that these subfamilies correlate with functional categories, one of which contains the bifunctional barley germin that has both OXO and superoxide dismutase (SOD) activity. It is proposed that GLPs function primarily as SODs, enzymes that protect plants from the effects of oxidative stress. Closer inspection of the DNA sequence encoding the intermotif region in plant GLPs showed global conservation of thymine in the second codon position, a character associated with hydrophobic residues. Since many of these proteins are multimeric and enzymatically inactive in their monomeric state, this conservation of hydrophobicity is thought to be associated with the need to maintain the various monomer- monomer interactions. The type of structure-based predictive analysis presented in this paper is an important approach for understanding gene function and evolution in an era when genomes from a wide range of organisms are being sequenced at a rapid rate.

Reduced incorporation of the influenza B virus BM2 protein in virus particles decreases infectivity

BM2 is the fourth integral membrane protein encoded by the influenza B virus genome. It is synthesized late in infection and transported to the plasma membrane from where it is subsequently incorporated into progeny virus particles. It has recently been reported that BM2 has ion channel activity and may be the functional homologue of the influenza A virus M2 protein acting as an ion channel involved in viral entry. Using a reverse genetic approach it was not possible to recover virus which lacked BM2. A recombinant influenza B virus was generated in which the BM2 AUG initiation codon was mutated to GUG. This decreased the efficiency of translation of BM2 protein such that progeny virions contained only 1/8 the amount of BM2 seen in wild-type virus. The reduction in BM2 incorporation resulted in a reduction in infectivity although there was no concomitant decrease in the numbers of virions released from the infected cells. These data imply that the incorporation of sufficient BM2 protein into influenza B virions is required for infectivity of the virus particles. (C) 2004 Elsevier Inc. All rights reserved.

A block in virus-like particle maturation following assembly of murine leukaemia virus in insect cells

Expression of the murine leukaemia virus (MLV) major Gag antigen p65(Gag) using the baculovirus expression system leads to efficient assembly and release of virus-like particles (VLP) representative of immature MLV. Expression of P180(Gag-Pol), facilitated normally in mammalian cells by readthrough of the p65(Gag) termination codon, also occurs efficiently in insect cells to provide a source of the MLV protease and a pattern of p65(Gag) processing similar to that observed in mammalian cells. VLP release from P180(Gag-Pol) expressing cells however remains essentially immature with disproportionate levels of the uncleaved p65(Gag) precursor when compared to the intracellular Gag profile. Changing the p65(Gag) termination codon altered the level of p65(Gag) and p180(Gag-Pol) within expressing cells but did not alter the pattern of released VLP, which remained immature. Coexpression of p65(Gag) with a fixed readthrough p180(Gag-Pol) also led to only immature VLP release despite high intracellular protease levels. Our data suggest a mechanism that preferentially selects uncleaved p65(Gag) for the assembly of MLV in this heterologous expression system and implies that, in addition to their relative levels, active sorting of the correct p65(Gag) and p180(Gag-Pol) ratios may occur in producer cells. (C) 2003 Elsevier Inc. All rights reserved.

HAWAIIAN SKIRT: an F-box gene that regulates organ fusion and growth in arabidopsis

A fast neutron-mutagenized population of Arabidopsis ( Arabidopsis thaliana) Columbia-0 wild-type plants was screened for floral phenotypes and a novel mutant, termed hawaiian skirt ( hws), was identified that failed to shed its reproductive organs. The mutation is the consequence of a 28 bp deletion that introduces a premature amber termination codon into the open reading frame of a putative F-box protein ( At3g61590). The most striking anatomical characteristic of hws plants is seen in flowers where individual sepals are fused along the lower part of their margins. Crossing of the abscission marker, Pro(PGAZAT):beta-glucuronidase, into the mutant reveals that while floral organs are retained it is not the consequence of a failure of abscission zone cells to differentiate. Anatomical analysis indicates that the fusion of sepal margins precludes shedding even though abscission, albeit delayed, does occur. Spatial and temporal characterization, using Pro(HWS):beta-glucuronidase or Pro(HWS):green fluorescent protein fusions, has identified HWS expression to be restricted to the stele and lateral root cap, cotyledonary margins, tip of the stigma, pollen, abscission zones, and developing seeds. Comparative phenotypic analyses performed on the hws mutant, Columbia-0 wild type, and Pro(35S):HWS ectopically expressing lines has revealed that loss of HWS results in greater growth of both aerial and below-ground organs while overexpressing the gene brings about a converse effect. These observations are consistent with HWS playing an important role in regulating plant growth and development.

The Glu298Asp single nucleotide polymorphism in the endothelial nitric oxide synthase gene differentially effects the vascular response to acute consumption of fruit and vegetable puree-based drinks

Scope Diets low in fruits and vegetables (FV) are responsible for 2.7 million deaths from cardiovascular diseases (CVD) and certain cancers annually. Many FV and their juices contain flavonoids, some of which increase endothelial nitric oxide synthase (eNOS) activity. A single nucleotide polymorphism in the eNOS gene, where thymine (T) replaces guanine (G) at position 894 predicting substitution of glutamate for aspartate at codon 298 (Glu298Asp), has been associated with increased CVD risk due to effects on nitric oxide synthesis and subsequently vascular reactivity. Individuals can be homozygous for guanine (GG), thymine (TT) or heterozygous (GT). Methods and results We investigated the effects of acute ingestion of a FV-puree-based-drink (FVPD) on vasodilation and antioxidant status in subjects retrospectively genotyped for this polymorphism. Healthy volunteers (n = 24; 11 GG, 11 GT, 2 TT) aged 30–70 were recruited to a randomized, controlled, crossover, acute study. We showed that acute consumption of 400 mL FVPD differentially affected individuals depending on their genotype. There was a significant genotype interaction for endothelium-dependent vasodilation measured by laser Doppler imaging with iontophoresis (P < 0.05) and ex vivo low-density lipoproteins (LDL) oxidation (P = 0.002). GG subjects had increased endothelium-dependent vasodilation 180 min (P = 0.028) and reduced ex vivo LDL oxidation (P = 0.013) after 60 min after FVPD compared with control, no differences were observed in GT subjects. Conclusion eNOS Glu298Asp genotype differentially affects vasodilation and ex vivo LDL oxidation after consumption of FV in the form of a puree-based drink.

Evolutionary dynamics of Staphylococcus aureus during progression from carriage to disease

Whole-genome sequencing offers new insights into the evolution of bacterial pathogens and the etiology of bacterial disease. Staph- ylococcus aureus is a major cause of bacteria-associated mortality and invasive disease and is carried asymptomatically by 27% of adults. Eighty percent of bacteremias match the carried strain. How- ever, the role of evolutionary change in the pathogen during the progression from carriage to disease is incompletely understood. Here we use high-throughput genome sequencing to discover the genetic changes that accompany the transition from nasal carriage to fatal bloodstream infection in an individual colonized with meth- icillin-sensitive S. aureus. We found a single, cohesive population exhibiting a repertoire of 30 single-nucleotide polymorphisms and four insertion/deletion variants. Mutations accumulated at a steady rate over a 13-mo period, except for a cluster of mutations preceding the transition to disease. Although bloodstream bacteria differed by just eight mutations from the original nasally carried bacteria, half of those mutations caused truncation of proteins, including a prema- ture stop codon in an AraC-family transcriptional regulator that has been implicated in pathogenicity. Comparison with evolution in two asymptomatic carriers supported the conclusion that clusters of pro- tein-truncating mutations are highly unusual. Our results demon- strate that bacterial diversity in vivo is limited but nonetheless detectable by whole-genome sequencing, enabling the study of evolutionary dynamics within the host. Regulatory or structural changes that occur during carriage may be functionally important for pathogenesis; therefore identifying those changes is a crucial step in understanding the biological causes of invasive bacterial disease.

Contingency locus in ctsR of Listeria monocytogenes Scott A: a strategy for occurrence of abundant piezotolerant isolates within clonal populations

In a recent study we demonstrated that a high-hydrostatic-pressure-tolerant isolate of Listeria monocytogenes lacks a codon in the class 3 heat shock regulator gene ctsR. This mutation in the region that encodes four consecutive glycines was directly responsible for the observed piezotolerance, increased stress resistance, and reduced virulence. The aim of the present study was to determine whether mutations in ctsR are frequently associated with piezotolerance in L. monocytogenes. Wild-type cultures of L. monocytogenes were therefore exposed to 350 MPa for 20 min, and the piezotolerance of individual surviving isolates was assessed. This rendered 33 isolates with a stable piezotolerant phenotype from a total of 84 survivors. Stable piezotolerant mutants were estimated to be present in the initial wild-type population at frequencies of >10�5. Subsequent sequencing of the ctsR gene of all stable piezotolerant isolates revealed that two-thirds of the strains (i.e., n � 21) had mutations in this gene. The majority of the mutations (16 of 21 strains) consisted of a triplet deletion in the glycine-encoding region of ctsR, identical to what was found in our previous study. Interestingly, 2 of 21 mutants contained a codon insertion in this repeat region. The remaining three stable piezotolerant strains showed a 19-bp insertion in the glycine repeat region, a 16-bp insertion downstream of the glycine repeat area (both leading to frameshifts and a truncated ctsR), and an in-frame 114-bp deletion encoding a drastically shortened carboxy terminus of CtsR. In four instances it was not possible to generate a PCR product. A piezotolerant phenotype could not be linked to mutations in ctsR in 8 of 33 isolates, indicating that other thus-far-unknown mechanisms also lead to stable piezotolerance. The present study highlights the importance of ctsR in piezotolerance and stress tolerance of L. monocytogenes, and it demonstrates that short-sequence repeat regions contribute significantly to the occurrence of a piezotolerant and stress-tolerant subpopulation within L. monocytogenes cultures, thus playing an important role in survival.

The CtsR regulator of Listeria monocytogenes contains a variant glycine repeat region that affects piezotolerance, stress resistance, motility and virulence

A spontaneous high hydrostatic pressure (HHP)-tolerant mutant of Listeria monocytogenes ScottA, named AK01, was isolated previously. This mutant was immotile and showed increased resistance to heat, acid and H2O2 compared with the wild type (wt) (Karatzas, K.A.G. and Bennik, M.H.J. 2002 Appl Environ Microbiol 68: 3183–3189). In this study, we conclusively linked the increased HHP and stress tolerance of strain AK01 to a single codon deletion in ctsR (class three stress gene repressor) in a region encoding a highly conserved glycine repeat. CtsR negatively regulates the expression of the clp genes, including clpP, clpE and the clpC operon (encompassing ctsR itself), which belong to the class III heat shock genes. Allelic replacement of the ctsR gene in the wt background with the mutant ctsR gene, designated ctsRΔGly, rendered mutants with phenotypes and protein expression profiles identical to those of strain AK01. The expression levels of CtsR, ClpC and ClpP proteins were significantly higher in ctsRΔGly mutants than in the wt strain, indicative of the CtsRΔGly protein being inactive. Further evidence that the CtsRΔGly protein lacks its repressor function came from the finding that the Clp proteins in the mutant were not further induced upon heat shock, and that HHP tolerance of a ctsR deletion strain was as high as that of a ctsRΔGly mutant. The high HHP tolerance possibly results from the increased expression of the clp genes in the absence of (active) CtsR repressor. Importantly, the strains expressing CtsRΔGly show significantly attenuated virulence compared with the wt strain; however, no indication of disregulation of PrfA in the mutant strains was found. Our data highlight an important regulatory role of the glycine-rich region of CtsR in stress resistance and virulence.

Prevalence of mutations within the quinolone resistance-determining region of gyrA, gyrB, parC, and parE and association with antibiotic resistance in quinolone-resistant Salmonella enterica

Salmonella enterica isolates (n = 182) were examined for mutations in the quinolone resistance-determining region of gyrA, gyrB, parC, and parE. The frequency, location, and type of GyrA substitution varied with the serovar. Mutations were found in parC that encoded Thr57-Ser, Thr66-Ile, and Ser80-Arg substitutions. Mutations in the gyrB quinolone resistance-determining region were located at codon Tyr420-Cys or Arg437-Len. Novel mutations were also found in parE encoding Glu453-Gly, His461-Tyr, Ala498-Thr, Val512-Gly, and Ser518-Cys. Although it is counterintuitive, isolates with a mutation in both gyrA and parC were more susceptible to ciprofloxacin than were isolates with a mutation in gyrA alone.

Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Proposal for a unified nomenclature for target site mutations associated with resistance to fungicides

Evolved resistance to fungicides is a major problem limiting our ability to control agricultural, medical and veterinary pathogens and is frequently associated with substitutions in the amino acid sequence of the target protein. The convention for describing amino-acid substitutions is to cite the wild type amino acid, the codon number and the new amino acid, using the one letter amino acid code. It has frequently been observed that orthologous amino acid mutations have been selected in different species by fungicides from the same mode of action class, but the amino acids have different numbers. These differences in numbering arise from the different lengths of the proteins in each species. The purpose of the current paper is to propose a system for unifying the labelling of amino acids in fungicide target proteins. To do this we have produced alignments between fungicide target proteins of relevant species fitted to a well-studied “archetype” species. Orthologous amino acids in all species are then assigned numerical “labels” based on the position of the amino acid in the archetype protein.