Experimental evolution and genome sequencing reveal variation in levels of clonal interference in large populations of bacteriophage φX174

Background

The first step of GPI anchor biosynthesis is catalyzed by PIG-A, an enzyme that transfers N -acetylglucosamine from UDP- N -acetylglucosamine to phosphatidylinositol. This protein is present in all eukaryotic organisms ranging from protozoa to higher mammals, as part of a larger complex of five to six 'accessory' proteins whose individual roles in the glycosyltransferase reaction are as yet unclear. The PIG-A gene has been shown to be an essential gene in various eukaryotes. In humans, mutations in the protein have been associated with paroxysomal noctural hemoglobuinuria. The corresponding PIG-A gene has also been recently identified in the genome of many archaeabacteria although genes of the accessory proteins have not been discovered in them. The present study explores the evolution of PIG-A and the phylogenetic relationship between this protein and other glycosyltransferases.

Results

In this paper we show that out of the twelve conserved motifs identified by us eleven are exclusively present in PIG-A and, therefore, can be used as markers to identify PIG-A from newly sequenced genomes. Three of these motifs are absent in the primitive eukaryote, G. lamblia. Sequence analyses show that seven of these conserved motifs are present in prokaryote and archaeal counterparts in rudimentary forms and can be used to differentiate PIG-A proteins from glycosyltransferases. Using partial least square regression analysis and data involving presence or absence of motifs in a range of PIG-A and glycosyltransferases we show that (i) PIG-A may have evolved from prokaryotic glycosyltransferases and lipopolysaccharide synthases, members of the GT4 family of glycosyltransferases and (ii) it is possible to uniquely classify PIG-A proteins versus glycosyltransferases.

Conclusion

Besides identifying unique motifs and showing that PIG-A protein from G. lamblia and some putative PIG-A proteins from archaebacteria are evolutionarily closer to glycosyltransferases, these studies provide a new method for identification and classification of PIG-A proteins.     

Genome-wide Comparative Analysis of Annexin Superfamily in Plants

Most annexins are calcium-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signaling during plant growth and development. They have previously been identified and characterized in Arabidopsis and rice, and constitute a multigene family in plants. In this study, we performed a comparative analysis of annexin gene families in the sequenced genomes of Viridiplantae ranging from unicellular green algae to multicellular plants, and identified 149 genes. Phylogenetic studies of these deduced annexins classified them into nine different arbitrary groups. The occurrence and distribution of bona fide type II calcium binding sites within the four annexin domains were found to be different in each of these groups. Analysis of chromosomal distribution of annexin genes in rice, Arabidopsis and poplar revealed their localization on various chromosomes with some members also found on duplicated chromosomal segments leading to gene family expansion. Analysis of gene structure suggests sequential or differential loss of introns during the evolution of land plant annexin genes. Intron positions and phases are well conserved in annexin genes from representative genomes ranging from Physcomitrella to higher plants. The occurrence of alternative motifs such as K/R/HGD was found to be overlapping or at the mutated regions of the type II calcium binding sites indicating potential functional divergence in certain plant annexins. This study provides a basis for further functional analysis and characterization of annexin multigene families in the plant lineage.     

Analysis of a Streptococcus pneumoniae gene encoding signal peptidase I and overproduction of the enzyme

The spi gene of Streptococcus pneumoniae was cloned and its nucleotide sequence was determined. It encodes a protein of 204 amino acids that is homologous to bacterial signal peptidase I proteins. The S. pneumoniae protein contains all of the conserved amino acid sequence motifs previously identified in this enzyme from both prokaryotic and eukaryotic sources. Sequence comparisons revealed several additional motifs characteristic of the enzyme. The cloned S. pneumoniae gene complemented an Escherichia coli mutant defective in its leader peptidase gene. Expression of the spi gene in S. pneumoniae appeared to be essential for viability. The cloned gene was shown to produce a polypeptide of approximately 20 kDa. Overproduction of the S. pneumoniae spi gene in an E. coli expression system gave a native protein product, soluble in the presence of a non-ionic detergent, which should be amenable to structural determination.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

Exploring the active site of tripeptidyl-peptidase II through studies of pH dependence of reaction kinetics

Tripeptidyl-peptidase II (TPP II) is a subtilisin-like serine protease which forms a large enzyme complex (> 4 MDa). It is considered a potential drug target due to its involvement in specific physiological processes. However, information is scarce concerning the kinetic characteristics of TPP II and its active site features, which are important for design of efficient inhibitors. To amend this, we probed the active site by determining the pH dependence of TPP II catalysis. Access to pure enzyme is a prerequisite for kinetic investigations and herein we introduce the first efficient purification system for heterologously expressed mammalian TPP II. The pH dependence of kinetic parameters for hydrolysis of two different chromogenic substrates, Ala-Ala-Phe-pNA and Ala-Ala-Ala-pNA, was determined for murine, human and Drosophila melanogaster TPP II as well as mutant variants thereof. The investigation demonstrated that TPP II, in contrast to subtilisin, has a bell-shaped pH dependence of k_cat^app/K_M probably due to deprotonation of the N-terminal amino group of the substrate at higher pH. Since both the K_M and k_cat^app are lower for cleavage of AAA-pNA than for AAF-pNA we propose that the former can bind non-productively to the active site of the enzyme, a phenomenon previously observed with some substrates for subtilisin. Two mutant variants, H267A and D387G, showed bell-shaped pH-dependence of k_cat^app, possibly due to an impaired protonation of the leaving group. This work reveals previously unknown differences between TPP II orthologues and subtilisin as well as features that might be conserved within the entire family of subtilisin-like serine peptidases.     

cumA Multicopper Oxidase Genes from Diverse Mn(II)-Oxidizing and Non-Mn(II)-Oxidizing Pseudomonas Strains

A multicopper oxidase gene, cumA, required for Mn(II) oxidation was recently identified in Pseudomonas putida strain GB-1. In the present study, degenerate primers based on the putative copper-binding regions of the cumA gene product were used to PCR amplify cumA gene sequences from a variety of Pseudomonas strains, including both Mn(II)-oxidizing and non-Mn(II)-oxidizing strains. The presence of highly conserved cumA gene sequences in several apparently non-Mn(II)-oxidizing Pseudomonas strains suggests that this gene may not be expressed, may not be sufficient alone to confer the ability to oxidize Mn(II), or may have an alternative function in these organisms. Phylogenetic analysis of both CumA and 16S rRNA sequences revealed similar topologies between the respective trees, including the presence of several distinct phylogenetic clusters. Overall, our results indicate that both the cumA gene and the capacity to oxidize Mn(II) occur in phylogenetically diverse Pseudomonas strains.     

G-Quadruplex DNA Sequences Are Evolutionarily Conserved and Associated with Distinct Genomic Features in Saccharomyces cerevisiae

G-quadruplex DNA is a four-stranded DNA structure formed by non-Watson-Crick base pairing between stacked sets of four guanines. Many possible functions have been proposed for this structure, but its in vivo role in the cell is still largely unresolved. We carried out a genome-wide survey of the evolutionary conservation of regions with the potential to form G-quadruplex DNA structures (G4 DNA motifs) across seven yeast species. We found that G4 DNA motifs were significantly more conserved than expected by chance, and the nucleotide-level conservation patterns suggested that the motif conservation was the result of the formation of G4 DNA structures. We characterized the association of conserved and non-conserved G4 DNA motifs in Saccharomyces cerevisiae with more than 40 known genome features and gene classes. Our comprehensive, integrated evolutionary and functional analysis confirmed the previously observed associations of G4 DNA motifs with promoter regions and the rDNA, and it identified several previously unrecognized associations of G4 DNA motifs with genomic features, such as mitotic and meiotic double-strand break sites (DSBs). Conserved G4 DNA motifs maintained strong associations with promoters and the rDNA, but not with DSBs. We also performed the first analysis of G4 DNA motifs in the mitochondria, and surprisingly found a tenfold higher concentration of the motifs in the AT-rich yeast mitochondrial DNA than in nuclear DNA. The evolutionary conservation of the G4 DNA motif and its association with specific genome features supports the hypothesis that G4 DNA has in vivo functions that are under evolutionary constraint.     

Identification of the nature of reading frame transitions observed in prokaryotic genomes

Our goal was to identify evolutionary conserved frame transitions in protein coding regions and to uncover an underlying functional role of these structural aberrations. We used the ab initio frameshift prediction program, GeneTack, to detect reading frame transitions in 206 991 genes (fs-genes) from 1106 complete prokaryotic genomes. We grouped 102 731 fs-genes into 19 430 clusters based on sequence similarity between protein products (fs-proteins) as well as conservation of predicted position of the frameshift and its direction. We identified 4010 pseudogene clusters and 146 clusters of fs-genes apparently using recoding (local deviation from using standard genetic code) due to possessing specific sequence motifs near frameshift positions. Particularly interesting was finding of a novel type of organization of the dnaX gene, where recoding is required for synthesis of the longer subunit, τ. We selected 20 clusters of predicted recoding candidates and designed a series of genetic constructs with a reporter gene or affinity tag whose expression would require a frameshift event. Expression of the constructs in Escherichia coli demonstrated enrichment of the set of candidates with sequences that trigger genuine programmed ribosomal frameshifting; we have experimentally confirmed four new families of programmed frameshifts.     

Analyses of phylogeny, evolution, conserved sequences and genome-wide expression of the ICK/KRP family of plant CDK inhibitors

Background and Aims

The cell cycle is controlled by cyclin-dependent kinases (CDKs), and CDK inhibitors are major regulators of their activities. The ICK/KRP family of CDK inhibitors has been reported in several plants, with seven members in arabidopsis; however, the phylogenetic relationship among members in different species is unknown. Also, there is a need to understand how these genes and proteins are regulated. Furthermore, little information is available on the functional differences among ICK/KRP family members.

Methods

We searched publicly available databases and identified over 120 unique ICK/KRP protein sequences from more than 60 plant species. Phylogenetic analysis was performed using 101 full-length sequences from 40 species and intron–exon organization of ICK/KRP genes in model species. Conserved sequences and motifs were analysed using ICK/KRP protein sequences from arabidopsis (Arabidopsis thaliana), rice (Orysa sativa) and poplar (Populus trichocarpa). In addition, gene expression was examined using microarray data from arabidopsis, rice and poplar, and further analysed by RT-PCR for arabidopsis.

Key Results and Conclusions

Phylogenetic analysis showed that plant ICK/KRP proteins can be grouped into three major classes. Whereas the C-class contains sequences from dicotyledons, monocotyledons and gymnosperms, the A- and B-classes contain only sequences from dicotyledons or monocotyledons, respectively, suggesting that the A- and B-classes might have evolved from the C-class. This classification is also supported by exon–intron organization. Genes in the A- and B- classes have four exons, whereas genes in the C-class have only three exons. Analysis of sequences from arabidopsis, rice and poplar identified conserved sequence motifs, some of which had not been described previously, and putative functional sites. The presence of conserved motifs in different family members is consistent with the classification. In addition, gene expression analysis showed preferential expression of ICK/KRP genes in certain tissues. A model has been proposed for the evolution of this gene family in plants.     

Genome-Wide Phylogenetic Analysis of Stress-Activated Protein Kinase Genes in Rice (OsSAPKs) and Expression Profiling in Response to Xanthomonas oryzae pv. oryzicola Infection

All members of the SnRK2 protein kinase gene family encoded by the rice (Oryza sativa L.) genome are activated by hyperosmotic stress, and have been designated as stress-activated protein kinases (SAPKs). In this study, gene structures, phylogeny, and conserved motifs for the entire OsSAPK gene family in rice have been analyzed. Moreover, expression patterns of OsSAPK in response to infection with Xanthomonas oryzae pv. oryzicola (Xoc) were investigated. A total of ten OsSAPK genes in the japonica rice cultivar 9804 were identified and classified into four groups. All genes had similar exon–intron structures and organization of putative motifs/domains, and shared the same four motifs (motifs 1–4). Group I (OsSAPK1 and OsSAPK2) shared another two motifs (motif 5 and motif 10), while group III (OsSAPK8, OsSAPK9 and OsSAPK10) had seven motifs in common (motifs 1–7). Moreover, we found that four OsSAPKs, including OsSAPK3, OsSAPK5, OsSAPK7 and OsSAPK9, were significantly upregulated in response to infection by Xoc in rice plants carrying the nonhost resistance gene Rxo1. Four of the OsSAPK genes in which expression was upregulated were localized to both the cytoplasm and nucleus, but clustered in different groups, suggesting that they are involved in different resistance signal transduction pathways. These results will provide useful information for the future functional dissection of this gene family.     

Molecular characterization of the lysosomal acid phosphatase fromDrosophila melanogaster

InDrosophila, unlike humans, the lysosomal acid phosphatase (Acph-1) is a non-essential enzyme. It is also one of the most rapidly evolving gene-enzyme systems in the genus. In order to determine which parts of the enzyme are conserved and which parts are apparently under little functional constraint, we cloned the gene fromDrosophila melanogaster via a chromosomal walk. Fragments from the gene were used to recover an apparently full-length cDNA. The cDNA was subcloned into aDrosophila transformation vector where it was under the control of the 5′ promoter sequence of thehsp-70 gene. Three independent transformants were obtained; in each, Acph-1 expression from the cDNA was constitutive and not dependent on heat shock, as determined by densitometric analyses of the allozymic forms of the enzyme. The pattern of expression indicates thehsp-70 and endogenousAcph-1 promoters act together in some, but not all, tissues. The sequence of the cDNA was determined using deletions made with exonuclease III, and primers deduced from the cDNA sequence were used to sequence the genomic clone. Five introns were found, and putative 5′ up-stream regulatory sequences were identified. Amino acid sequence comparisons have revealed several highly conserved motifs betweenDrosophila Acph-1 and vertebrate lysosomal and prostatic acid phosphatases.     

Predicting regulons and their cis-regulatory motifs by comparative genomics

We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caernorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.     

Serine/threonine kinases and E2-ubiquitin conjugating enzymes in Planctomycetes: unexpected findings

The regulation of signal transduction by phosphorylation and ubiquitination is essential to ensure the correct behavior of eukaryotic cells. We searched for protein families involved in such signaling in several eukaryotic species and in a limited set of prokaryotes, where two members of the Planctomycetes phylum were included as they exhibit eukaryote-like features (Gemmata obscuriglobus and Pirellula staleyi). We identified sequences homologous to eukaryotic serine/threonine kinases (STKs) and E2-ubiquitin conjugating enzymes in the two Planctomycetes species. To extend these analyses to the Planctomycetes/Verrucomicrobia/Chlamydia super-phylum, we performed comparative analyses using domains from kinases, phosphatases and GTPases that serve as signaling signatures, and we analyzed their distributions. We found substantial differences in kinome densities with regards to other prokaryote clades and among the groups in the Planctomycetes/Verrucomicrobia/Chlamydia super-phylum. In addition, we identified the presence of classic eukaryotic E2-conjugating ubiquitin proteins in prokaryotes, these having previously believed to exist only in eukaryotes. Our phylogenetic analyses of the STKs signature domains and E2-enzymes suggest the existence of horizontal gene transfer.     

Alternative Transposition Generates New Chimeric Genes and Segmental Duplications at the Maize p1 Locus

The maize Ac/Ds transposon family was the first transposable element system identified and characterized by Barbara McClintock. Ac/Ds transposons belong to the hAT family of class II DNA transposons. We and others have shown that Ac/Ds elements can undergo a process of alternative transposition in which the Ac/Ds transposase acts on the termini of two separate, nearby transposons. Because these termini are present in different elements, alternative transposition can generate a variety of genome alterations such as inversions, duplications, deletions, and translocations. Moreover, Ac/Ds elements transpose preferentially into genic regions, suggesting that structural changes arising from alternative transposition may potentially generate chimeric genes at the rearrangement breakpoints. Here we identified and characterized 11 independent cases of gene fusion induced by Ac alternative transposition. In each case, a functional chimeric gene was created by fusion of two linked, paralogous genes; moreover, each event was associated with duplication of the ∼70-kb segment located between the two paralogs. An extant gene in the maize B73 genome that contains an internal duplication apparently generated by an alternative transposition event was also identified. Our study demonstrates that alternative transposition-induced duplications may be a source for spontaneous creation of diverse genome structures and novel genes in maize.     

In silico analysis of PHB gene family in maize

Prohibitins (PHBs) are highly conserved proteins in species ranging from prokaryotes to eukaryotes. Plant PHBs have been implicated in various cellular processes including development, senescence and stress responses. Although PHBs have been investigated in several plant species including Arabidopsis and tobacco, no systematic gene family analysis has been carried in maize. In the present study, 16 putative PHB genes have been identified. Analysis of the conserved protein motifs and gene structures has revealed high levels of conservation within the phylogenetic subgroups. Published microarray database showed that most maize PHB genes exhibited different expression levels in different tissues and developmental stages. Cis-elements analysis showed that ZmPHB2 and ZmPHB12 may play important roles in plant development. Taken together, we provide a comprehensive bioinformatics analysis of the PHB gene family in maize genome and our data provide an important foundation for further functional study of this gene family in maize.