Gene Re-annotation in Genome of the Extremophile Pyrobaculum Aerophilum by Using Bioinformatics Methods

Abstract

In this paper, we re-annotated the genome of Pyrobaculum aerophilum str. IM2, particularly for hypothetical ORFs. The annotation process includes three parts. Firstly and most importantly, 23 new genes, which were missed in the original annotation, are found by combining similarity search and the ab initio gene finding approaches. Among these new genes, five have significant similarities with function-known genes and the rest have significant similarities with hypothetical ORFs contained in other genomes. Secondly, the coding potentials of the 1645 hypothetical ORFs are re-predicted by using 33 Z curve variables combined with Fisher linear discrimination method. With the accuracy being 99.68%, 25 originally annotated hypothetical ORFs are recognized as non-coding by our method. Thirdly, 80 hypothetical ORFs are assigned with potential functions by using similarity search with BLAST program. Re-annotation of the genome will benefit related researches on this hyperthermophilic crenarchaeon. Also, the re-annotation procedure could be taken as a reference for other archaeal genomes. Details of the revised annotation are freely available at http://cobi.uestc.edu.cn/resource/paero/     

A class of highly polymorphic tetranucleotide repeats for canine genetic mapping

We have identified and characterized a new class of polymorphic markers for the canine genome from a simple tetranucleotide repeat sequence, (GAAA)_n. Genetic markers derived from this repeat are highly polymorphic compared with other canine microsatellites, yet are stable enough to be useful for following Mendelian inheritance in multigeneration pedigrees. We show further that (GAAA)_n repeats are distributed throughout the canine genome and occur with sufficient frequency to be useful in the development of a framework map of the canine genome. Received: 25 October 1995 / Accepted: 17 January 1996     

An evidence of laccases in archaea

Laccases (benzenediol:oxygen oxidoreductase, EC 1.10.3.2) are a diverse group of multicopper oxidases that catalyze the oxidation of a variety of aromatic compounds. Here we present evidence for distribution of laccases among archaea and their probable functions. Putative laccase genes have been found in different archaeal groups that might have branched off early during evolution, e.g. Haloarcula marismortui ATCC 43049, Natronomonas pharaonis DSM2160, Pyrobaculum aerophilum IM2, Candidatus Nitrosopumilus maritimus SCM1, Halorubrum lacusprofundi ATCC 49239. Most of the archaeal multicopper oxidases reported here are of Type 1 and Type 2 whereas type 3 copper-binding domain could be found in Pyrobaculum aerophilum IM2 and Halorubrum lacusprofundi ATCC49239. An analysis of the genome sequence database revealed the presence of novel types of two-domain laccases in archaea. ed using this method. CyMVin the positive samples of Phalaenopsis sp. and Arachnis sp. was confirmed by DNA sequencing and cp gene homeology blast. The results showed that CyMV extracted from the leaves of orchid in Hangzhou, Zhejiang Province, China, could be derived from Kunming city (KM), Yunnan Province, China. This method characterized by high sensitivity, specificity, and precision is suitable for early diagnosis and quantitative detection of CyMV.     

Random sheared fosmid library as a new genomic tool to accelerate complete finishing of rice (Oryza sativa spp. Nipponbare) genome sequence: sequencing of gap-specific fosmid clones uncovers new euchromatic portions of the genome

The International Rice Genome Sequencing Project has recently announced the high-quality finished sequence that covers nearly 95% of the japonica rice genome representing 370 Mbp. Nevertheless, the current physical map of japonica rice contains 62 physical gaps corresponding to approximately 5% of the genome, that have not been identified/represented in the comprehensive array of publicly available BAC, PAC and other genomic library resources. Without finishing these gaps, it is impossible to identify the complete complement of genes encoded by rice genome and will also leave us ignorant of some 5% of the genome and its unknown functions. In this article, we report the construction and characterization of a tenfold redundant, 40 kbp insert fosmid library generated by random mechanical shearing. We demonstrated its utility in refining the physical map of rice by identifying and in silico mapping 22 gap-specific fosmid clones with particular emphasis on chromosomes 1, 2, 6, 7, 8, 9 and 10. Further sequencing of 12 of the gap-specific fosmid clones uncovered unique rice genome sequence that was not previously reported in the finished IRGSP sequence and emphasizes the need to complete finishing of the rice genome. Fosmid library reported here is publicly available from our web site     

Purification and characterization of the tungsten enzyme aldehyde:ferredoxin oxidoreductase from the hyperthermophilic denitrifier Pyrobaculum aerophilum

A tungsten-containing aldehyde:ferredoxin oxidoreductase (AOR) has been purified to homogeneity from Pyrobaculum aerophilum. The N-terminal sequence of the isolated enzyme matches a single open reading frame in the genome. Metal analysis and electron paramagnetic resonance (EPR) spectroscopy indicate that the P. aerophilum AOR contains one tungsten center and one [4Fe-4S]^2+/1+ cluster per 68-kDa monomer. Native AOR is a homodimer. EPR spectroscopy of the purified enzyme that has been reduced with the substrate crotonaldehyde revealed a W(V) species with g_zyx values of 1.952, 1.918, 1.872. The substrate-reduced AOR also contains a [4Fe-4S]¹⁺ cluster with S=3/2 and zero field splitting parameters D=7.5 cm^–1 and E/D=0.22. Molybdenum was absent from the enzyme preparation. The P. aerophilum AOR lacks the amino acid sequence motif indicative for binding of mononuclear iron that is typically found in other AORs. Furthermore, the P. aerophilum AOR utilizes a 7Fe ferredoxin as the putative physiological redox partner, instead of a 4Fe ferredoxin as in Pyrococcus furiosus. This 7Fe ferredoxin has been purified from P. aerophilum, and the amino acid sequence has been identified using mass spectrometry. Direct electrochemistry of the ferredoxin showed two one-electron transitions, at –306 and –445 mV. In the presence of 55 M ferredoxin the AOR activity is 17% of the activity obtained with 1 mM benzyl viologen as an electron acceptor.     

Pa-AGOG, the founding member of a new family of archaeal 8-oxoguanine DNA-glycosylases

Oxidative damage represents a major threat to genomic stability, as the major product of DNA oxidation, 8-oxoguanine (GO), frequently mispairs with adenine during replication. In order to prevent these mutagenic events, organisms have evolved GO-DNA glycosylases that remove this oxidized base from DNA. We were interested to find out how GO is processed in the hyperthermophilic archaeon Pyrobaculum aerophilum, which lives at temperatures around 100°C. To this end, we searched its genome for open reading frames (ORFs) bearing the principal hallmark of GO-DNA glycosylases: a helix–hairpin–helix motif and a glycine/proline-rich sequence followed by an absolutely conserved aspartate (HhH-GPD motif). Interestingly, although the P.aerophilum genome encodes three such ORFs, none of these encodes the potent GO-processing activity detected in P.aerophilum extracts. Fractionation of the extracts, followed by analysis of the active fractions by denaturing polyacrylamide gel electrophoresis, showed that the GO-processing enzyme has a molecular size of ~30 kDa. Mass spectrometric analysis of proteins in this size range identified several peptides originating from P.aerophilum ORF PAE2237. We now show that PAE2237 encodes AGOG (Archaeal GO-Glycosylase), the founding member of a new family of DNA glycosylases, which can remove GO from single- and double-stranded substrates with great efficiency.     

Unique active site formation in a novel galactose 1-phosphate uridylyltransferase from the hyperthermophilic archaeon Pyrobaculum aerophilum

A gene encoding galactose 1-phosphate uridylyltransferase (GalT) was identified in the hyperthermophilic archaeon Pyrobaculum aerophilum. The gene was overexpressed in Escherichia coli, after which its product was purified and characterized. The expressed enzyme was highly thermostable and retained about 90% of its activity after incubation for 10 minutes at temperatures up to 90°C. Two different crystal structures of P. aerophilum GalT were determined: the substrate-free enzyme at 2.33 Å and the UDP-bound H140F mutant enzyme at 1.78 Å. The main-chain coordinates of the P. aerophilum GalT monomer were similar to those in the structures of the E. coli and human GalTs, as was the dimeric arrangement. However, there was a striking topological difference between P. aerophilum GalT and the other two enzymes. In the E. coli and human enzymes, the N-terminal chain extends from one subunit into the other and forms part of the substrate-binding pocket in the neighboring subunit. By contrast, the N-terminal chain in P. aerophilum GalT extends to the substrate-binding site in the same subunit. Amino acid sequence alignment showed that a shorter surface loop in the N-terminal region contributes to the unique topology of P. aerophilum GalT. Structural comparison of the substrate-free enzyme with UDP-bound H140F suggests that binding of the glucose moiety of the substrate, but not the UDP moiety, gives rise to a large structural change around the active site. This may in turn provide an appropriate environment for the enzyme reaction.     

Molecular cloning and characterisation of the ribC gene from Bacillus subtilis : a point mutation in ribC results in riboflavin overproduction

A mutation leading to roseoflavin resistance and deregulated riboflavin biosynthesis was mapped in the genome of the riboflavin-overproducing Bacillus subtilis strains RB52 and RB50 at map position 147°. The chromosomal location indicates that the deregulating mutation in RB52 and RB50 is an allele of the previously identified ribC mutation. We cloned the ribC gene and found that it encodes a putative 36-kDa protein. Surprisingly, RibC has significant sequence similarity to flavin kinases and FAD synthases from various other bacterial species. By comparing the deduced amino acid sequence of RibC from the wild-type parent strain of RB50 with the RibC sequence from the riboflavin-overexpressing RB50 mutant we identified a point mutation that resulted in a Gly to Ser exchange in the C-terminal region of the product Received: 3 June 1996 / Accepted: 19 October 1996     

Simple sequence repeats for the genetic analysis of apple

 The development of highly informative markers, such as simple sequence repeats, for tagging genes controlling agronomic characters is essential for apple breeding. Furthermore the use of these markers is fundamental both for variety identification and for the characterisation and management of genetic resources. We have developed 16 reliable simple sequence repeat (SSR) markers that amplify all alleles from a panel of 19 Malus x domestica (Borkh.) cultivars or breeding selections and from Malus floribunda 821. Those markers show a high level of genetic polymorphism, with on average 8.2 alleles per locus and an average heterozygosity of 0.78. Due to this high level of polymorphism, it was possible using two selected SSRs to distinguish all cultivars except Starking and Red Delicious. Ten of the markers we developed have been mapped on a RAPD linkage map, proving their Mendelian segregation as well as their random distribution in the apple genome. Finally, we discuss the importance of using co-dominant markers in outbreeding species. Received: 8 October 1997 / Accepted: 9 December 1997     

Complete Sequence of SpHox8 and Its Linkage in the Single Hox Gene Cluster of Strongylocentrotus purpuratus

SpHox8 is the paralog group 8 Hox gene of Strongylocentrotus purpuratus. This identification follows from an analysis of the sequence of the complete open reading frame of a late-gastrula-stage cDNA clone; from direct linkage to adjacent Hox genes in a fosmid contig; and from blot hybridizations carried out on pulse field gel electrophoretic separations. SpHox8 is a single-copy gene, and there is only one Hox gene cluster per genome in S. purpuratus. Received: 7 November 1996 / Accepted: 4 December 1996     

A thermostable vacuolar-type membrane pyrophosphatase from the archaeon Pyrobaculum aerophilum: implications for the origins of pyrophosphate-energized pumps

Vacuolar-type H⁺-translocating pyrophosphatases (V-PPases) have been considered to be restricted to plants, a few species of phototrophic proteobacteria and protists. Here, we describe PVP, a thermostable, sequence-divergent V-PPase from the facultatively aerobic hyperthermophilic archaeon Pyrobaculum aerophilum. PVP shares only 38% sequence identity with both the prototypical V-PPase from Arabidopsis thaliana and the H⁺-PPi synthase from Rhodospirillum rubrum, yet possesses most of the structural features characteristic of V-PPases. Heterologous expression of PVP in Saccharomyces cerevisiae yields a M_r 64 000 membrane polypeptide that specifically catalyzes Mg²⁺-dependent PPi hydrolysis. The existence of PVP implies that PPi-energized H⁺-translocation is phylogenetically more deeply rooted than previously thought.     

Position of the reduced mycorrhizal colonisation (<Emphasis Type="Italic">Rmc</Emphasis>) locus on the tomato genome map

Our research aims to investigate the molecular communication between land plants and arbuscular mycorrhizal (AM) fungi in the establishment of symbiosis. We have identified a mutation in the facultative AM host tomato, which we named rmc. Plants that are homozygous for rmc no longer host most AM fungi. The mutation also affects the interaction of tomato with root knot nematode and Fusarium wilt. However, the function/s encoded by the intact Rmc locus is/are unknown. To clone and sequence the gene or genes that comprise the Rmc locus, we have initiated a positional cloning project. In this paper, we report the construction of mapping populations and use of molecular markers from the published genome map to identify the location of Rmc on tomato chromosome 8. Nucleotide binding site-leucine rich repeat resistance genes, reported to reside in the same region of that chromosome, provided insufficient differences to develop cleaved amplified polymorphic sequence markers. Therefore, we were unable to map these sequences in relation to rmc. Our results potentiate future work to identify the Rmc function and to determine the genetic basis for the multiple plant-microbe interaction functions that the rmc mutation has defined.     

Molecular linkage map for an intraspecific recombinant inbred population of durum wheat (Triticum turgidum L. var. durum)

Durum wheat (Triticum turgidum L. var. durum) is an economically and nutritionally important cereal crop in the Mediterranean region. To further our understanding of durum genome organization we constructed a durum linkage map using restriction fragment length polymorphisms (RFLPs), simple sequence repeats (SSRs) known as Gatersleben wheat microsatellites (GWMs), amplified fragment length polymorphisms (AFLPs), and seed storage proteins (SSPs: gliadins and glutenins). A population of 110 F₉ recombinant inbred lines (RILs) was derived from an intraspecific cross between two durum cultivars, Jennah Khetifa and Cham 1. The two parents exhibit contrasting traits for resistance to biotic and abiotic stresses and for grain quality. In total, 306 markers have been placed on the linkage map – 138 RFLPs, 26 SSRs, 134 AFLPs, five SSPs, and three known genes (one pyruvate decarboxylase and two lipoxygenases). The map is 3598 cM long, with an average distance between markers of 11.8 cM, and 12.1% of the markers deviated significantly from the expected Mendelian ratio 1:1. The molecular markers were evenly distributed between the A and B genomes. The chromosome with the most markers is 1B (41 markers), followed by 3B and 7B, with 25 markers each. The chromosomes with the fewest markers are 2A (11 markers), 5A (12 markers), and 4B (15 markers). In general, there is a good agreement between the map obtained and the Triticeae linkage consensus maps. This intraspecific map provides a useful tool for marker-assisted selection and map-based breeding for resistance to biotic and abiotic stresses and for improvement of grain quality. Received: 14 February 2000 / Accepted: 28 April 2000     

Construction of an oilseed rape (Brassica napus L.) genetic map with SSR markers

We constructed a Brassica napus genetic map with 240 simple sequence repeats (SSR) primer pairs from private and public origins. SSR, or microsatellites, are highly polymorphic and efficient markers for the analysis of plant genomes. Our selection of primer pairs corresponded to 305 genetic loci that we were able to map. In addition, we also used 52 sequence-characterized amplified region primer pairs corresponding to 58 loci that were developed in our lab. Genotyping was performed on six F2 populations, corresponding to a total of 574 F2 individual plants, obtained according to an unbalanced diallel cross design involving six parental lines. The resulting consensus map presented 19 linkage groups ranging from 46.2 to 276.5 cM, which we were able to name after the B. napus map available at , thus enabling the identification of the A genome linkage groups originating from the B. rapa ancestor and the C genome linkage groups originating from the B. oleracea ancestor in the amphidiploid genome of B. napus. Some homoeologous regions were identified between the A and the C genomes. This map could be used to identify more markers, which would eventually be linked to genes controlling important agronomic characters in rapeseed. Furthermore, considering the good genome coverage we obtained, together with an observed homogenous distribution of the loci across the genome, this map is a powerful tool to be used in marker-assisted breeding. Electronic Supplementary Material Supplementary material is available for this article at     

Gene-specific universal mammalian sequence-tagged sites: Application to the canine genome

We are developing a genetic map of the dog based partly upon markers contained within known genes. In order to facilitate the development of these markers, we have used polymerase chain reaction (PCR) primers designed to conserved regions of genes that have been sequenced in at least two species. We have refined the method for designing primers to maximize the number that produce successful amplifications across as many mammalian species as possible. We report the development of primer sets for 11 loci in detail:CFTR, COL10A1, CSFIR, CYP1A1, DCN1, FES, GHR, GLB1, PKLR, PVALB, andRB1. We also report an additional 75 primer sets in the appendices. The PCR products were sequenced to show that the primers amplify the expected canine genes. These primer sets thus define a class of gene-specific sequence-tagged sites (STSs). There are a number of uses for these STSs, including the rapid development of various linkage tools and the rapid testing of genomic and cDNA libraries for the presence of their corresponding genes. Six of the eleven gene targets reported in detail have been proposed to serve as “anchored reference loci” for the development of mammalian genetic maps [O'Brien, S. J.,et al., Nat. Genet. 3:103, 1993]. The primer sets should cover a significant portion of the canine genome for the development of a linkage map. In order to determine how useful these primer sets would be for the other genome projects, we tested the 11 primer sets on the DNA from species representing five mammalian orders. Eighty-four percent of the gene-species combinations amplified successfully. We have named these primer sets “universal mammalian sequence-tagged sites” because they should be useful for many mammalian genome projects.     

A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content

We have developed a new DH mapping population for oilseed rape, named TNDH, using genetically and phenotypically diverse parental lines. We used the population in the construction of a high stringency genetic linkage map, consisting of 277 loci, for use in quantitative genetic analysis. A proportion of the markers had been used previously in the construction of linkage maps for Brassica species, thus permitting the alignment of maps. The map includes 68 newly developed Sequence Tagged Site (STS) markers targeted to the homologues of defined genes of A. thaliana. The use of these markers permits the alignment of our linkage map with the A. thaliana genome sequence. An additional 74 loci (31 newly developed STS markers and 43 loci defined by SSR and RFLP markers that had previously been used in published linkage maps) were added to the map. These markers increased the resolution of alignment of the newly constructed linkage map with existing Brassica linkage maps and the A. thaliana genome sequence. We conducted field trials with the TNDH population at two sites, and over 2 years, and identified reproducible QTL for seed oil content and erucic acid content. The results provide new insights into the genetic control of seed oil and erucic acid content in oilseed rape, and demonstrate the utility of the linkage map and population.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.D. Qiu and C. Morgan authors contributed equally to the work.     

Fosmid Library Construction and Initial Analysis of End Sequences in Zhikong Scallop (<Emphasis Type="Italic">Chlamys farreri</Emphasis>)

Zhikong scallop (Chlamys farreri Jones et Preston, 1904) is one of the most commercially important bivalves in China, but research on its genome is underdeveloped. In this study, we constructed the first Zhikong scallop fosmid library, and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 133,851 clones with an average insert size of about 40 kb, amounting to 4.3 genome equivalents. Fosmid stability assays indicate that Zhikong scallop DNA was stable during propagation in the fosmid system. Library screening with two genes and seven microsatellite markers yielded between two and eight positive clones, and none of those tested was absent from the library. End-sequencing of 480 individual clones generated 828 sequences after trimming, with an average sequence length of 624 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 213 (25.72%) and 44 (5.31%) significant hits (E < e⁻⁵), respectively. Repetitive sequences analysis resulted in 375 repeats, accounting for 15.84% of total length, which were composed of interspersed repetitive sequences, tandem repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for physical mapping and positional cloning, and provide a better understanding of the Zhikong scallop genome.     

Biochemical characterization of DNA-binding proteins from Pyrobaculum aerophilum and Aeropyrum pernix

Several representatives of the Crenarchaeal branch of the Archaea contain highly abundant, small, positively charged proteins exemplified by the Sso7d protein from Sulfolobus solfataricus. These proteins bind to DNA in a non-sequence-specific manner. Using publicly available genomic sequence information, we identified a second class of small Crenarchaeal DNA-binding proteins represented by the Pyrobaculum aerophilum open reading frame 3192–encoded (Pae3192) protein and its paralogs. We investigated the biochemical properties of the Pae3192 protein and an orthologous protein (Ape1322b) from Aeropyrum pernix in side-by-side experiments with the Sso7d protein. We demonstrate that the recombinant Ape1322b, Pae3192 and Sso7d proteins bind to DNA and that the DNA-protein complexes formed are slightly different for each protein. We show that like Sso7d, Pae3192 constrains negative supercoils in DNA. In addition, we show that all three proteins raise the melting temperature of duplex DNA upon binding. Finally, we present the equilibrium affinity constants and kinetic association constants of each protein for single-stranded and double-stranded DNA.     

Retrotransposons and Tandem Repeat Sequences in the Nuclear Genomes of Cryptomonad Algae

The cryptomonads are an enigmatic group of unicellular eukaryotic algae that possess two nuclear genomes, having acquired photosynthesis by the uptake and retention of a eukaryotic algal endosymbiont. The endosymbiont nuclear genome, or nucleomorph, of the cryptomonad Guillardia theta has been completely sequenced: at only 551 kilobases (kb) and with a gene density of ∼1 gene/kb, it is a model of compaction. In contrast, very little is known about the structure and composition of the cryptomonad host nuclear genome. Here we present the results of two small-scale sequencing surveys of fosmid clone libraries from two distantly related cryptomonads, Rhodomonas salina CCMP1319 and Cryptomonas paramecium CCAP977/2A, corresponding to ∼150 and ∼235 kb of sequence, respectively. Very few of the random end sequences determined in this study show similarity to known genes in other eukaryotes, underscoring the considerable evolutionary distance between the cryptomonads and other eukaryotes whose nuclear genomes have been completely sequenced. Using a combination of fosmid clone end-sequencing, Southern hybridizations, and PCR, we demonstrate that Ty3-gypsy long-terminal repeat (LTR) retrotransposons and tandem repeat sequences are a prominent feature of the nuclear genomes of both organisms. The complete sequence of a 30.9-kb genomic fragment from R. salina was found to contain a full-length Ty3-gypsy element with near-identical LTRs and a chromodomain, a protein module suggested to mediate the site-specific integration of the retrotransposon. The discovery of chromodomain-containing retroelements in cryptomonads further expands the known distribution of the so-called chromoviruses across the tree of eukaryotes. [Reviewing Editor: Dr. Debashish Bhattacharya]