Complete genome sequence of the alkaliphilic bacterium Bacillus halodurans and genomic sequence comparison with Bacillus subtilis

The 4 202 353 bp genome of the alkaliphilic bacterium Bacillus halodurans C-125 contains 4066 predicted protein coding sequences (CDSs), 2141 (52.7%) of which have functional assignments, 1182 (29%) of which are conserved CDSs with unknown function and 743 (18.3%) of which have no match to any protein database. Among the total CDSs, 8.8% match sequences of proteins found only in Bacillus subtilis and 66.7% are widely conserved in comparison with the proteins of various organisms, including B.subtilis. The B.halodurans genome contains 112 transposase genes, indicating that transposases have played an important evolutionary role in horizontal gene transfer and also in internal genetic rearrangement in the genome. Strain C-125 lacks some of the necessary genes for competence, such as comS, srfA and rapC, supporting the fact that competence has not been demonstrated experimentally in C-125. There is no paralog of tupA, encoding teichuronopeptide, which contributes to alkaliphily, in the C-125 genome and an ortholog of tupA cannot be found in the B.subtilis genome. Out of 11 σ factors which belong to the extracytoplasmic function family, 10 are unique to B.halodurans, suggesting that they may have a role in the special mechanism of adaptation to an alkaline environment.     

The complete genomic sequence of Mycoplasma penetrans,an intracellular bacterial pathogen in humans

The complete genomic sequence of an intracellular bacterial pathogen, Mycoplasma penetrans HF-2 strain, was determined. The HF-2 genome consists of a 1 358 633 bp single circular chromosome containing 1038 predicted coding sequences (CDSs), one set of rRNA genes and 30 tRNA genes. Among the 1038 CDSs, 264 predicted proteins are common to the Mycoplasmataceae sequenced thus far and 463 are M.penetrans specific. The genome contains the two-component system but lacks the essential cellular gene, uridine kinase. The relatively large genome of M.penetrans HF-2 among mycoplasma species may be accounted for by both its rich core proteome and the presence of a number of paralog families corresponding to 25.4% of all CDSs. The largest paralog family is the p35 family, which encodes surface lipoproteins including the major antigen, P35. A total of 44 genes for p35 and p35 homologs were identified and 30 of them form one large cluster in the chromosome. The genetic tree of p35 paralogs suggests the occurrence of dynamic chromosomal rearrangement in paralog formation during evolution. Thus, M.penetrans HF-2 may have acquired diverse repertoires of antigenic variation-related genes to allow its persistent infection in humans.     

Life on Arginine for Mycoplasma hominis: Clues from Its Minimal Genome and Comparison with Other Human Urogenital Mycoplasmas

Mycoplasma hominis is an opportunistic human mycoplasma. Two other pathogenic human species, M. genitalium and Ureaplasma parvum, reside within the same natural niche as M. hominis: the urogenital tract. These three species have overlapping, but distinct, pathogenic roles. They have minimal genomes and, thus, reduced metabolic capabilities characterized by distinct energy-generating pathways. Analysis of the M. hominis PG21 genome sequence revealed that it is the second smallest genome among self-replicating free living organisms (665,445 bp, 537 coding sequences (CDSs)). Five clusters of genes were predicted to have undergone horizontal gene transfer (HGT) between M. hominis and the phylogenetically distant U. parvum species. We reconstructed M. hominis metabolic pathways from the predicted genes, with particular emphasis on energy-generating pathways. The Embden–Meyerhoff–Parnas pathway was incomplete, with a single enzyme absent. We identified the three proteins constituting the arginine dihydrolase pathway. This pathway was found essential to promote growth in vivo. The predicted presence of dimethylarginine dimethylaminohydrolase suggested that arginine catabolism is more complex than initially described. This enzyme may have been acquired by HGT from non-mollicute bacteria. Comparison of the three minimal mollicute genomes showed that 247 CDSs were common to all three genomes, whereas 220 CDSs were specific to M. hominis, 172 CDSs were specific to M. genitalium, and 280 CDSs were specific to U. parvum. Within these species-specific genes, two major sets of genes could be identified: one including genes involved in various energy-generating pathways, depending on the energy source used (glucose, urea, or arginine) and another involved in cytadherence and virulence. Therefore, a minimal mycoplasma cell, not including cytadherence and virulence-related genes, could be envisaged containing a core genome (247 genes), plus a set of genes required for providing energy. For M. hominis, this set would include 247+9 genes, resulting in a theoretical minimal genome of 256 genes.     

Complete Genome Sequence of the Plant Pathogen Erwinia amylovora Strain ATCC 49946

Erwinia amylovora causes the economically important disease fire blight that affects rosaceous plants, especially pear and apple. Here we report the complete genome sequence and annotation of strain ATCC 49946. The analysis of the sequence and its comparison with sequenced genomes of closely related enterobacteria revealed signs of pathoadaptation to rosaceous hosts.Erwinia amylovora, a plant-associated member of the Enterobacteriaceae, causes fire blight, a devastating disease of rosaceous plants, especially pear and apple (6). The complete genome of Ea273 (ATCC 49946), a virulent strain isolated from an infected apple tree in New York State, was sequenced. Total DNA was extracted and prepared in pMAQ1 shotgun libraries. The complete shotgun sequence was obtained by using dye terminator chemistry in ABI 3730 automated sequencers and contains 88,457 reads (11.12-fold coverage), yielding a theoretical coverage of the genome of 99.99%. The sequence was assembled, finished, and annotated as described previously (1, 5), using Artemis (4) to collate data and facilitate annotation.The genome of E. amylovora consists of a circular chromosome of 3,805,874 bp and two plasmids, AMYP1 (28,243 bp) and AMYP2 (71,487 bp). Coding regions in the chromosome account for 85.1% of the total sequence, with 3,483 identified coding sequences (CDS). Two hundred fifty-four (7%) of the CDSs do not have any matches in current NCBI databases; 114 (3.3%) correspond to conserved hypothetical proteins. Forty-nine CDSs (1.4%) are similar to genes from mobile elements such as integrases, transposases, and bacteriophages, and 110 CDSs (3.2%) were classified as pseudogenes due to interruptions or truncations of the CDSs. The remaining 2,956 annotated CDSs include among other categories genes involved in biosynthesis of the cellular envelope and modifications of surface proteins (299 CDSs [11%]) and genes involved in signal transduction and regulation (228 CDSs [8%]). Seven rRNA operons and 78 tRNA sequences were identified in the chromosome; two new clusters were identified (AMY1550-1575 and AMY2648-2676) that resemble the T3SS-encoding SSR-1 island of Sodalis glossinidius (2), and four clusters that contain genes for biosynthesis of flagella, which based on their location might be regulated independently.The smaller plasmid, AMYP1, had been reported as pEA29 (3); its sequence is nearly identical to the one reported here. The larger plasmid, AMYP2, renamed pEA72 for consistency in nomenclature, contains 87 predicted CDSs, with two predicted mobile-element-related CDSs and one pseudogene. Among the CDSs with annotated functions are a cluster of genes (AMYP2_49 to AMYP2_62) that encode a putative type IV fimbrial system (pil genes).The genome of E. amylovora is only 3.8 Mb long, whereas most free-living enterobacteria, including plant pathogens, have genomes of 4.5 Mb to 5.5 Mb. Comparison of the genome of Ea273 with the sequenced genomes of 15 closely related enterobacteria identified 21 lineage-specific regions, which might be considered genomic islands. E. amylovora has many more predicted pseudogenes, relative to other enterobacteria with similar lifestyles. Given its size and the preponderance of pseudogenes, genome reduction may have occurred via mutational inactivation and subsequent deletion with the following consequences: E. amylovora has fewer genes involved in anaerobic respiration and fermentation than are found in typical related enterobacteria; this likely result in a reduced capacity to live in anaerobic environments.The genome sequence of E. amylovora has revealed clear signs of pathoadaptation to the rosaceous plant environment. For example, T3SS-related proteins are present that are more similar to proteins of other plant pathogens than to proteins of closely related enterobacteria. These include type III effectors, homologous to those of plant-pathogenic pseudomonads, which confer virulence to E. amylovora in plants, and a sorbitol-metabolizing cluster that may confer a competitive advantage for survival in rosaceous plants. The reduced genome size and erosion or loss of genes involved in anaerobic respiration and nitrate assimilation are remarkable, relative to other plant- and animal-pathogenic members of the Enterobacteriaceae.     

Complete Genome Sequence of Halalkalicoccus jeotgali B3T,an Extremely Halophilic Archaeon

Halalkalicoccus jeotgali B3^T, isolated from salt-fermented seafood from South Korea, is an extremely halophilic archaeon belonging to the family Halobacteriaceae. Here, we present the complete genome sequence of the type strain H. jeotgali B3^T (3,698,650 bp, with a G+C content of 62.5%), which consists of one chromosome and six plasmids. This is the first complete genome sequence of the Halalkalicoccus species.Extremely halophilic archaea (haloarchaea) are adapted to hypersaline environments and grow optimally in NaCl solutions of 2.6 M or higher (12). These haloarchaea are classified within the family Halobacteriaceae in the order Halobacteriales; currently, this family comprises 28 genera (3), and only 11 complete genome sequences in Halobacteriaceae have been reported. In a study of archaeal diversity in salt-fermented small shrimp or shellfish from South Korea, our laboratory isolated and characterized 5 novel, extremely halophilic archaeal strains of Halobacteriaceae. These strains included Natronococcus jeotgali (9), Halalkalicoccus jeotgali (11), Halorubrum cibi (7), Haloterrigena jeotgali (10) and Haladaptatus cibarius (8). We have now sequenced the genome of Halalkalicoccus jeotgali B3^T; genome sequencing had not been completed or initiated for any strain in this genus when our sequencing project was begun. The genus Halalkalicoccus currently contains only two species, Halalkalicoccus tibetensis (13) and H. jeotgali, and these species exhibit 98.6% gene sequence similarity in their 16S rRNA. The genome of H. jeotgali B3^T is the first of this genus to be sequenced.The complete genome sequence of H. jeotgali B3^T was determined by a whole-genome shotgun strategy using Roche 454 GS (FLX Titanium) pyrosequencing (898,168 reads totaling ∼348 Mb; ∼94-fold coverage of the genome) and a fosmid library (514 reads totaling ∼680 kb) at the Genome Resource Center, KRIBB (Korea Research Institute of Bioscience and Biotechnology). Genome sequences from pyrosequencing were processed by Roche''s software according to the manufacturer''s instructions, and sequences from the fosmid library were processed by PESTAS (6). A total of 898,196 reads were assembled using Newbler Assembler 2.3 (454 Life Science), which generated 54 large contigs (>100 bp in size) with bases having quality scores of 40 and above. The gaps between contigs were closed by primer walking and sequencing of PCR products across the gaps. The annotation was done by merging results obtained from the RAST (Rapid Annotation using Subsystem Technology) pipeline (1), Glimmer 3.02 (2), tRNAscan-SE 1.21 (5), and RNAmmer 1.2 (4).The H. jeotgali B3^T genome is 3,698,650 bases long with a 62.5% G+C content. The chromosome consists of a single circular chromosome (2,809,118 bp, with a G+C content of 65.0%) and six plasmids (406,285 bp, 55.3%; 363,534 bp, 54.2%; 44,576 bp, 58.9%; 44,459 bp, 54.9%; 23,727 bp, 47.6%; 6,951 bp, 60.6%). The genome contains 3,860 predicted coding sequences and 52 RNA genes (determined using RAST). The chromosome is predicted to contain 3,101 coding sequences with a coding intensity of 90.0%, including 47 tRNA genes, 1 5S rRNA gene, 1 16S rRNA gene, and 1 23S rRNA gene. The largest plasmid contains 466 coding sequences with a coding intensity of 81.2% and 2 tRNA genes, while the other five plasmids contain 425, 44, 48, 29, and 5 coding sequences with coding intensities of 80.2%, 84.2%, 83.0%, 69.6%, and 22.8%, respectively (determined using Glimmer3). More detailed analysis of this genome and comparative analysis with other haloarchaea will provide further insight into the genomic differences and metabolism of the extremely halophilic archaea.     

Distinct Gene Number-Genome Size Relationships for Eukaryotes and Non-Eukaryotes: Gene Content Estimation for Dinoflagellate Genomes

The ability to predict gene content is highly desirable for characterization of not-yet sequenced genomes like those of dinoflagellates. Using data from completely sequenced and annotated genomes from phylogenetically diverse lineages, we investigated the relationship between gene content and genome size using regression analyses. Distinct relationships between log₁₀-transformed protein-coding gene number (Y′) versus log₁₀-transformed genome size (X′, genome size in kbp) were found for eukaryotes and non-eukaryotes. Eukaryotes best fit a logarithmic model, Y′ = ln(-46.200+22.678X′, whereas non-eukaryotes a linear model, Y′ = 0.045+0.977X′, both with high significance (p<0.001, R²>0.91). Total gene number shows similar trends in both groups to their respective protein coding regressions. The distinct correlations reflect lower and decreasing gene-coding percentages as genome size increases in eukaryotes (82%–1%) compared to higher and relatively stable percentages in prokaryotes and viruses (97%–47%). The eukaryotic regression models project that the smallest dinoflagellate genome (3×10⁶ kbp) contains 38,188 protein-coding (40,086 total) genes and the largest (245×10⁶ kbp) 87,688 protein-coding (92,013 total) genes, corresponding to 1.8% and 0.05% gene-coding percentages. These estimates do not likely represent extraordinarily high functional diversity of the encoded proteome but rather highly redundant genomes as evidenced by high gene copy numbers documented for various dinoflagellate species.     

Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom⁺ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod⁺ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom⁺ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom⁺ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.     

Genome Sequence of a Cellulose-Producing Bacterium,Gluconacetobacter hansenii ATCC 23769

The Gram-negative bacterium Gluconacetobacter hansenii is considered a model organism for studying cellulose synthesis. We have determined the genome sequence of strain ATCC 23769.Plants produce cellulose, an unbranched chain of β-1,4-linked glucose units, as a structural polysaccharide. It is the most abundant polymer on earth, recently receiving much interest due to its potential use as a feedstock for bioethanol. Bacteria also produce cellulose. Among these, Gluconacetobacter hansenii (previously named Acetobacter xylinus) (4) has been extensively characterized and is a model system for cellulose biosynthesis (1, 2, 7). G. hansenii produces extracellular cellulose that is devoid of lignin or hemicellulose, making it an excellent source for pure cellulose. A lack of a completely sequenced genome for this organism has been a limiting factor in identifying other key proteins involved in cellulose synthesis.The whole-genome sequencing of G. hansenii ATCC 23769 was performed using the 454 FLX-Titanium pyrosequencing technology (5). A combinatorial sequencing approach using 489,201 reads obtained from the shotgun library and 195,088 reads from an 8-kb pair end library (3) produced a total of 221,294,116 bp. These reads were assembled using the Newbler assembler, producing 88 large contigs (>500 bp) and a chromosome-sized scaffold of 3,646,142 bp with an average coverage of ×50.5. This scaffold contained exclusively chromosomal DNA and no plasmid sequences. The gaps in the large scaffold were filled by primer walking and subsequent sequencing of the PCR products. The resulting high-quality draft assembly, consisting of a large scaffold with 71 contigs, was annotated using the Prokaryotic Genomes Automatic Annotation Pipeline (PGAAP) service of the National Institute of Biotechnology Information (NCBI).The chromosomal sequence of G. hansenii 23769 contains 3,547,122 bp, with a G+C content of 59%. The genome contains 3,351 genes, of which 3,308 are protein-encoding genes, accounting for 84% of the genome. There are 43 genes for tRNAs and 2 rRNA loci. The genes encoding proteins involved in cellulose synthesis are in an operon consisting of acsAB (GXY_04277), acsC (GXY_04282), and acsD (GXY_04292), as previously shown by Saxena et al. (7). Interestingly, there are two additional copies of acsAB, GXY_08864 and GXY_14452, which share 69% and 72% sequence identity, respectively, with the acsAB genes in the operon; the deduced amino acid sequences are 40% and 46% identical, respectively, with that deduced from acsAB in the operon. There are also two additional copies of acsC, GXY_08869 and GXY_014472, which share 72% and 65% DNA sequence identity, respectively, with the acsC gene in the operon; the deduced amino acid sequences share 28% and 30% amino acid identity, respectively, with that deduced from acsC. acsAB (GXY_08864) and acsC (GXY_08869) are only 17 bp apart, less than the distance (66 bp) between the acsAB and acsC genes in the operon. acsAB (GXY_14452) and acsC (GXY_14472) are separated by 3,299 bp, with three genes in between. However, acsD is present only in the operon, not duplicated elsewhere in the genome. The genome also contains three genes encoding diguanylate cyclase, as previously reported by Tal et al. (8). Diguanylate cyclase catalyzes the formation of cyclic di-GMP, a second messenger in bacteria that functions as an allosteric activator of cellulase synthase AcsAB (6).     

The genome sequence of the tomato-pathogenic actinomycete Clavibacter michiganensis subsp. michiganensis NCPPB382 reveals a large island involved in pathogenicity

Clavibacter michiganensis subsp. michiganensis is a plant-pathogenic actinomycete that causes bacterial wilt and canker of tomato. The nucleotide sequence of the genome of strain NCPPB382 was determined. The chromosome is circular, consists of 3.298 Mb, and has a high G+C content (72.6%). Annotation revealed 3,080 putative protein-encoding sequences; only 26 pseudogenes were detected. Two rrn operons, 45 tRNAs, and three small stable RNA genes were found. The two circular plasmids, pCM1 (27.4 kbp) and pCM2 (70.0 kbp), which carry pathogenicity genes and thus are essential for virulence, have lower G+C contents (66.5 and 67.6%, respectively). In contrast to the genome of the closely related organism Clavibacter michiganensis subsp. sepedonicus, the genome of C. michiganensis subsp. michiganensis lacks complete insertion elements and transposons. The 129-kb chp/tomA region with a low G+C content near the chromosomal origin of replication was shown to be necessary for pathogenicity. This region contains numerous genes encoding proteins involved in uptake and metabolism of sugars and several serine proteases. There is evidence that single genes located in this region, especially genes encoding serine proteases, are required for efficient colonization of the host. Although C. michiganensis subsp. michiganensis grows mainly in the xylem of tomato plants, no evidence for pronounced genome reduction was found. C. michiganensis subsp. michiganensis seems to have as many transporters and regulators as typical soil-inhabiting bacteria. However, the apparent lack of a sulfate reduction pathway, which makes C. michiganensis subsp. michiganensis dependent on reduced sulfur compounds for growth, is probably the reason for the poor survival of C. michiganensis subsp. michiganensis in soil.     

Ganglioside GM3 Levels Are Altered in a Mouse Model of HIBM: GM3 as a Cellular Marker of the Disease

Objective

HIBM (Hereditary Inclusion Body Myopathy) is a recessive hereditary disease characterized by adult-onset, slowly progressive muscle weakness sparing the quadriceps. It is caused by a single missense mutation of each allele of the UDP-N-acetylglucosamine 2-epimerase/N-acetylmannosamine kinase (GNE) gene, a bifunctional enzyme catalyzing the first two steps of sialic acid synthesis in mammals. However, the mechanisms and cellular pathways affected by the GNE mutation and causing the muscle weakness could not be identified so far. Based on recent evidence in literature, we investigated a new hypothesis, i.e. the involvement in the disease of the GM3 ganglioside, a specific glycolipid implicated in muscle cell proliferation and differentiation.

Methods

qRT-PCR analysis of St3gal5 (GM3 synthase) gene expression and HPLC quantification of GM3 ganglioside were conducted on muscle tissue from a mouse model of HIBM harboring the M712T mutation of GNE (Gne^M712T/M712T mouse) vs control mice (Gne^+/+ mouse).

Results

St3gal5 mRNA levels were significantly lower in Gne^M712T/M712T mouse muscles vs Gne^+/+ mouse muscles (64.41%±10% of Gne^+/+ levels). GM3 ganglioside levels showed also a significant decrease in Gne^M712T/M712T mouse muscle compared to Gne^+/+ mouse muscle (18.09%±5.33% of Gne^+/+ levels). Although these Gne^M712T/M712T mice were described to suffer severe glomerular proteinuria, no GM3 alterations were noted in kidneys, highlighting a tissue specific alteration of gangliosides.

Conclusion

The M712T mutation of GNE hampers the muscle ability to synthesize normal levels of GM3. This is the first time that a mutation of GNE can be related to the molecular pathological mechanism of HIBM.     

Genome Sequence of “Candidatus Mycoplasma haemolamae” Strain Purdue,a Red Blood Cell Pathogen of Alpacas (Vicugna pacos) and Llamas (Lama glama)

We report the complete genome sequence of “Candidatus Mycoplasma haemolamae,” an endemic red-cell pathogen of camelids. The single, circular chromosome has 756,845 bp, a 39.3% G+C content, and 925 coding sequences (CDSs). A great proportion (49.1%) of these CDSs are organized into paralogous gene families, which can now be further explored with regard to antigenic variation.     

Consensus-Degenerate Hybrid Oligonucleotide Primers for Amplification of Priming Glycosyltransferase Genes of the Exopolysaccharide Locus in Strains of the Lactobacillus casei Group

A primer design strategy named CODEHOP (consensus-degenerate hybrid oligonucleotide primer) for amplification of distantly related sequences was used to detect the priming glycosyltransferase (GT) gene in strains of the Lactobacillus casei group. Each hybrid primer consisted of a short 3′ degenerate core based on four highly conserved amino acids and a longer 5′ consensus clamp region based on six sequences of the priming GT gene products from exopolysaccharide (EPS)-producing bacteria. The hybrid primers were used to detect the priming GT gene of 44 commercial isolates and reference strains of Lactobacillus rhamnosus, L. casei, Lactobacillus zeae, and Streptococcus thermophilus. The priming GT gene was detected in the genome of both non-EPS-producing (EPS⁻) and EPS-producing (EPS⁺) strains of L. rhamnosus. The sequences of the cloned PCR products were similar to those of the priming GT gene of various gram-negative and gram-positive EPS⁺ bacteria. Specific primers designed from the L. rhamnosus RW-9595M GT gene were used to sequence the end of the priming GT gene in selected EPS⁺ strains of L. rhamnosus. Phylogenetic analysis revealed that Lactobacillus spp. form a distinctive group apart from other lactic acid bacteria for which GT genes have been characterized to date. Moreover, the sequences show a divergence existing among strains of L. rhamnosus with respect to the terminal region of the priming GT gene. Thus, the PCR approach with consensus-degenerate hybrid primers designed with CODEHOP is a practical approach for the detection of similar genes containing conserved motifs in different bacterial genomes.     

Genome Size, Complexity, and Ploidy of the Pathogenic Fungus Histoplasma capsulatum

The genome size, complexity, and ploidy of the dimorphic pathogenic fungus Histoplasma capsulatum was determined by using DNA renaturation kinetics, genomic reconstruction, and flow cytometry. Nuclear DNA was isolated from two strains, G186AS and Downs, and analyzed by renaturation kinetics and genomic reconstruction with three putative single-copy genes (calmodulin, α-tubulin, and β-tubulin). G186AS was found to have a genome of approximately 2.3 × 10⁷ bp with less than 0.5% repetitive sequences. The Downs strain, however, was found to have a genome approximately 40% larger with more than 16 times more repetitive DNA. The Downs genome was determined to be 3.2 × 10⁷ bp with approximately 8% repetitive DNA. To determine ploidy, the DNA mass per cell measured by flow cytometry was compared with the 1n genome estimate to yield a DNA index (DNA per cell/1n genome size). Strain G186AS was found to have a DNA index of 0.96, and Downs had a DNA index of 0.94, indicating that both strains are haploid. Genomic reconstruction and Southern blot data obtained with α- and β-tubulin probes indicated that some genetic duplication has occurred in the Downs strain, which may be aneuploid or partially diploid.     

Identification of stop codon readthrough genes in Saccharomyces cerevisiae

We specifically sought genes within the yeast genome controlled by a non-conventional translation mechanism involving the stop codon. For this reason, we designed a computer program using the yeast database genomic regions, and seeking two adjacent open reading frames separated only by a unique stop codon (called SORFs). Among the 58 SORFs identified, eight displayed a stop codon bypass level ranging from 3 to 25%. For each of the eight sequences, we demonstrated the presence of a poly(A) mRNA. Using isogenic [PSI⁺] and [psi^–] yeast strains, we showed that for two of the sequences the mechanism used is a bona fide readthrough. However, the six remaining sequences were not sensitive to the PSI state, indicating either a translation termination process independent of eRF3 or a new stop codon bypass mechanism. Our results demonstrate that the presence of a stop codon in a large ORF may not always correspond to a sequencing error, or a pseudogene, but can be a recoding signal in a functional gene. This emphasizes that genome annotation should take into account the fact that recoding signals could be more frequently used than previously expected.     

GAD2 on chromosome 10p12 is a candidate gene for human obesity

The gene GAD2 encoding the glutamic acid decarboxylase enzyme (GAD65) is a positional candidate gene for obesity on Chromosome 10p11–12, a susceptibility locus for morbid obesity in four independent ethnic populations. GAD65 catalyzes the formation of γ-aminobutyric acid (GABA), which interacts with neuropeptide Y in the paraventricular nucleus to contribute to stimulate food intake. A case-control study (575 morbidly obese and 646 control subjects) analyzing GAD2 variants identified both a protective haplotype, including the most frequent alleles of single nucleotide polymorphisms (SNPs) +61450 C>A and +83897 T>A (OR = 0.81, 95% CI [0.681–0.972], p = 0.0049) and an at-risk SNP (−243 A>G) for morbid obesity (OR = 1.3, 95% CI [1.053–1.585], p = 0.014). Furthermore, familial-based analyses confirmed the association with the obesity of SNP +61450 C>A and +83897 T>A haplotype (χ² = 7.637, p = 0.02). In the murine insulinoma cell line βTC3, the G at-risk allele of SNP −243 A>G increased six times GAD2 promoter activity (p < 0.0001) and induced a 6-fold higher affinity for nuclear extracts. The −243 A>G SNP was associated with higher hunger scores (p = 0.007) and disinhibition scores (p = 0.028), as assessed by the Stunkard Three-Factor Eating Questionnaire. As GAD2 is highly expressed in pancreatic β cells, we analyzed GAD65 antibody level as a marker of β-cell activity and of insulin secretion. In the control group, −243 A>G, +61450 C>A, and +83897 T>A SNPs were associated with lower GAD65 autoantibody levels (p values of 0.003, 0.047, and 0.006, respectively). SNP +83897 T>A was associated with lower fasting insulin and insulin secretion, as assessed by the HOMA-B% homeostasis model of β-cell function (p = 0.009 and 0.01, respectively). These data support the hypothesis of the orexigenic effect of GABA in humans and of a contribution of genes involved in GABA metabolism in the modulation of food intake and in the development of morbid obesity.     

Haploinsufficiency of Activation-Induced Deaminase for Antibody Diversification and Chromosome Translocations both In Vitro and In Vivo

The humoral immune response critically relies on the secondary diversification of antibodies. This diversification takes places through somatic remodelling of the antibody genes by two molecular mechanisms, Class Switch Recombination (CSR) and Somatic Hypermutation (SHM). The enzyme Activation Induced Cytidine Deaminase (AID) initiates both SHM and CSR by deaminating cytosine residues on the DNA of immunoglobulin genes. While crucial for immunity, AID-catalysed deamination is also the triggering event for the generation of lymphomagenic chromosome translocations. To address whether restricting the levels of AID expression in vivo contributes to the regulation of its function, we analysed mice harbouring a single copy of the AID gene (AID^+/−). AID^+/− mice express roughly 50% of normal AID levels, and display a mild hyperplasia, reminiscent of AID deficient mice and humans. Moreover, we found that AID^+/− cells have an impaired competence for CSR and SHM, which indicates that AID gene dose is limiting for its physiologic function. We next evaluated the impact of AID reduction in AID^+/− mice on the generation of chromosome translocations. Our results show that the frequency of AID-promoted c-myc/IgH translocations is reduced in AID^+/− mice, both in vivo and in vitro. Therefore, AID is haploinsufficient for antibody diversification and chromosome translocations. These findings suggest that limiting the physiologic levels of AID expression can be a regulatory mechanism that ensures an optimal balance between immune proficiency and genome integrity.     

Purification and Characterization of the Soluble Methane Monooxygenase of the Type II Methanotrophic Bacterium Methylocystis sp. Strain WI 14

Methane monooxygenase (MMO) catalyzes the oxidation of methane to methanol as the first step of methane degradation. A soluble NAD(P)H-dependent methane monooxygenase (sMMO) from the type II methanotrophic bacterium WI 14 was purified to homogeneity. Sequencing of the 16S rDNA and comparison with that of other known methanotrophic bacteria confirmed that strain WI 14 is very close to the genus Methylocystis. The sMMO is expressed only during growth under copper limitation (<0.1 μM) and with ammonium or nitrate ions as the nitrogen source. The enzyme exhibits a low substrate specificity and is able to oxidize several alkanes and alkenes, cyclic hydrocarbons, aromatics, and halogenic aromatics. It has three components, hydroxylase, reductase and protein B, which is involved in enzyme regulation and increases sMMO activity about 10-fold. The relative molecular masses of the native components were estimated to be 229, 41, and 18 kDa, respectively. The hydroxylase contains three subunits with relative molecular masses of 57, 43, and 23 kDa, which are present in stoichiometric amounts, suggesting that the native protein has an α₂β₂γ₂ structure. We detected 3.6 mol of iron per mol of hydroxylase by atomic absorption spectrometry. sMMO is strongly inhibited by Hg²⁺ ions (with a total loss of enzyme activity at 0.01 mM Hg²⁺) and Cu²⁺, Zn²⁺, and Ni²⁺ ions (95, 80, and 40% loss of activity at 1 mM ions). The complete sMMO gene sequence has been determined. sMMO genes from strain WI 14 are clustered on the chromosome and show a high degree of homology (at both the nucleotide and amino acid levels) to the corresponding genes from Methylosinus trichosporium OB3b, Methylocystis sp. strain M, and Methylococcus capsulatus (Bath).     

Evidence for horizontal transfer of SsuDAT1I restriction-modification genes to the Streptococcus suis genome

Different strains of Streptococcus suis serotypes 1 and 2 isolated from pigs either contained a restriction-modification (R-M) system or lacked it. The R-M system was an isoschizomer of Streptococcus pneumoniae DpnII, which recognizes nucleotide sequence 5′-GATC-3′. The nucleotide sequencing of the genes encoding the R-M system in S. suis DAT1, designated SsuDAT1I, showed that the SsuDAT1I gene region contained two methyltransferase genes, designated ssuMA and ssuMB, as does the DpnII system. The deduced amino acid sequences of M.SsuMA and M.SsuMB showed 70 and 90% identity to M.DpnII and M.DpnA, respectively. However, the SsuDAT1I system contained two isoschizomeric restriction endonuclease genes, designated ssuRA and ssuRB. The deduced amino acid sequence of R.SsuRA was 49% identical to that of R.DpnII, and R.SsuRB was 72% identical to R.LlaDCHI of Lactococcus lactis subsp. cremoris DCH-4. The four SsuDAT1I genes overlapped and were bounded by purine biosynthetic gene clusters in the following gene order: purF-purM-purN-purH-ssuMA-ssuMB-ssuRA-ssuRB-purD-purE. The G+C content of the SsuDAT1I gene region (34.1%) was lower than that of the pur region (48.9%), suggesting horizontal transfer of the SsuDAT1I system. No transposable element or long-repeat sequence was found in the flanking regions. The SsuDAT1I genes were functional by themselves, as they were individually expressed in Escherichia coli. Comparison of the sequences between strains with and without the R-M system showed that only the region from 53 bp upstream of ssuMA to 5 bp downstream of ssuRB was inserted in the intergenic sequence between purH and purD and that the insertion target site was not the recognition site of SsuDAT1I. No notable substitutions or insertions could be found, and the structures were conserved among all the strains. These results suggest that the SsuDAT1I system could have been integrated into the S. suis chromosome by an illegitimate recombination mechanism.     

Biochemical and Genetic Characterization of an Extracellular Protease from Pseudomonas fluorescens CY091

Pseudomonas fluorescens CY091 cultures produce an extracellular protease with an estimated molecular mass of 50 kDa. Production of this enzyme (designated AprX) was observed in media containing CaCl₂ or SrCl₂ but not in media containing ZnCl₂, MgCl₂, or MnCl₂. The requirement of Ca²⁺ (or Sr²⁺) for enzyme production was concentration dependent, and the optimal concentration for production was determined to be 0.35 mM. Following ammonium sulfate precipitation and ion-exchange chromatography, the AprX in the culture supernatant was purified to near electrophoretic homogeneity. Over 20% of the enzyme activity was retained in the AprX sample which had been heated in boiling water for 10 min, indicating that the enzyme is highly resistant to heat inactivation. The enzyme activity was almost completely inhibited in the presence of 1 mM 1,10-phenanthroline, but only 30% of the activity was inhibited in the presence of 1 mM EGTA. The gene encoding AprX was cloned from the genome of P. fluorescens CY091 by isolating cosmid clones capable of restoring the protease production in a nonproteolytic mutant of strain CY091. The genomic region of strain CY091 containing the aprX gene was located within a 7.3-kb DNA fragment. Analysis of the complete nucleotide sequence of this 7.3-kb fragment revealed the presence of a cluster of genes required for the production of extracellular AprX in P. fluorescens and Escherichia coli. The AprX protein showed 50 to 60% identity in amino acid sequence to the related proteases produced by Pseudomonas aeruginosa and Erwinia chrysanthemi. Two conserved sequence domains possibly associated with Ca²⁺ and Zn²⁺ binding were identified. Immediately adjacent to the aprX structural gene, a gene (inh) encoding a putative protease inhibitor and three genes (aprD, aprE, and aprF), possibly required for the transport of AprX, were also identified. The organization of the gene cluster involved in the synthesis and secretion of AprX in P. fluorescens CY091 appears to be somewhat different from that previously demonstrated in P. aeruginosa and E. chrysanthemi.     

Ferrous Iron-Dependent Volatilization of Mercury by the Plasma Membrane of Thiobacillus ferrooxidans

Of 100 strains of iron-oxidizing bacteria isolated, Thiobacillus ferrooxidans SUG 2-2 was the most resistant to mercury toxicity and could grow in an Fe²⁺ medium (pH 2.5) supplemented with 6 μM Hg²⁺. In contrast, T. ferrooxidans AP19-3, a mercury-sensitive T. ferrooxidans strain, could not grow with 0.7 μM Hg²⁺. When incubated for 3 h in a salt solution (pH 2.5) with 0.7 μM Hg²⁺, resting cells of resistant and sensitive strains volatilized approximately 20 and 1.7%, respectively, of the total mercury added. The amount of mercury volatilized by resistant cells, but not by sensitive cells, increased to 62% when Fe²⁺ was added. The optimum pH and temperature for mercury volatilization activity were 2.3 and 30°C, respectively. Sodium cyanide, sodium molybdate, sodium tungstate, and silver nitrate strongly inhibited the Fe²⁺-dependent mercury volatilization activity of T. ferrooxidans. When incubated in a salt solution (pH 3.8) with 0.7 μM Hg²⁺ and 1 mM Fe²⁺, plasma membranes prepared from resistant cells volatilized 48% of the total mercury added after 5 days of incubation. However, the membrane did not have mercury reductase activity with NADPH as an electron donor. Fe²⁺-dependent mercury volatilization activity was not observed with plasma membranes pretreated with 2 mM sodium cyanide. Rusticyanin from resistant cells activated iron oxidation activity of the plasma membrane and activated the Fe²⁺-dependent mercury volatilization activity of the plasma membrane.