Continuous Culture Adaptation of Methylobacterium extorquens AM1 and TK 0001 to Very High Methanol Concentrations.

The bio-economy relies on microbial strains optimized for efficient large scale production of chemicals and fuels from inexpensive and renewable feedstocks under industrial conditions. The reduced one carbon compound methanol, whose production does not involve carbohydrates needed for the feed and food sector, can be used as sole carbon and energy source by methylotrophic bacteria like Methylobacterium extorquens AM1. This strain has already been engineered to produce various commodity and high value chemicals from methanol. The toxic effect of methanol limits its concentration as feedstock to 1% v/v. We obtained M. extorquens chassis strains tolerant to high methanol via adaptive directed evolution using the GM3 technology of automated continuous culture. Turbidostat and conditional medium swap regimes were employed for the parallel evolution of the recently characterized strain TK 0001 and the reference strain AM1 and enabled the isolation of derivatives of both strains capable of stable growth with 10% methanol. The isolates produced more biomass at 1% methanol than the ancestor strains. Genome sequencing identified the gene metY coding for an O-acetyl-L-homoserine sulfhydrylase as common target of mutation. We showed that the wildtype enzyme uses methanol as substrate at elevated concentrations. This side reaction produces methoxine, a toxic homolog of methionine incorporated in polypeptides during translation. All mutated metY alleles isolated from the evolved populations coded for inactive enzymes, designating O-acetyl-L-homoserine sulfhydrylase as a major vector of methanol toxicity. A whole cell transcriptomic analysis revealed that genes coding for chaperones and proteases were upregulated in the evolved cells as compared with the wildtype, suggesting that the cells had to cope with aberrant proteins formed during the adaptation to increasing methanol exposure. In addition, the expression of ribosomal proteins and enzymes related to energy production from methanol like formate dehydrogenases and ATP synthases was boosted in the evolved cells upon a short-term methanol stress. D-lactate production from methanol by adapted cells overexpressing the native D-lactate dehydrogenase was quantified. A significant higher lactate yield was obtained compared with control cells, indicating an enhanced capacity of the cells resistant to high methanol to assimilate this one carbon feedstock more efficiently.

Genome sequences of Escherichia coli B strains REL606 and BL21(DE3).

Escherichia coli K-12 and B have been the subjects of classical experiments from which much of our understanding of molecular genetics has emerged. We present here complete genome sequences of two E. coli B strains, REL606, used in a long-term evolution experiment, and BL21(DE3), widely used to express recombinant proteins. The two genomes differ in length by 72,304 bp and have 426 single base pair differences, a seemingly large difference for laboratory strains having a common ancestor within the last 67 years. Transpositions by IS1 and IS150 have occurred in both lineages. Integration of the DE3 prophage in BL21(DE3) apparently displaced a defective prophage in the lambda attachment site of B. As might have been anticipated from the many genetic and biochemical experiments comparing B and K-12 over the years, the B genomes are similar in size and organization to the genome of E. coli K-12 MG1655 and have >99% sequence identity over approximately 92% of their genomes. E. coli B and K-12 differ considerably in distribution of IS elements and in location and composition of larger mobile elements. An unexpected difference is the absence of a large cluster of flagella genes in B, due to a 41 kbp IS1-mediated deletion. Gene clusters that specify the LPS core, O antigen, and restriction enzymes differ substantially, presumably because of horizontal transfer. Comparative analysis of 32 independently isolated E. coli and Shigella genomes, both commensals and pathogenic strains, identifies a minimal set of genes in common plus many strain-specific genes that constitute a large E. coli pan-genome.

Contrasted microcolinearity and gene evolution within a homoeologous region of wheat and barley species.

We study here the evolution of genes located in the same physical locus using the recently sequenced Ha locus in seven wheat genomes in diploid, tetraploid, and hexaploid species and compared them with barley and rice orthologous regions. We investigated both the conservation of microcolinearity and the molecular evolution of genes, including coding and noncoding sequences. Microcolinearity is restricted to two groups of genes (Unknown gene-2, VAMP, BGGP, Gsp-1, and Unknown gene-8 surrounded by several copies of ATPase), almost conserved in rice and barley, but in a different relative position. Highly conserved genes between wheat and rice run along with genes harboring different copy numbers and highly variable sequences between close wheat genomes. The coding sequence evolution appeared to be submitted to heterogeneous selective pressure and intronic sequences analysis revealed that the molecular clock hypothesis is violated in most cases.

Genome structure impacts molecular evolution at the AvrLm1 avirulence locus of the plant pathogen Leptosphaeria maculans.

Leptosphaeria maculans, a dothideomycete fungus causing stem canker on oilseed rape, develops gene-for-gene interactions with its host plants. It has the ability to rapidly adapt to selection pressure exerted by cultivars harbouring novel resistance genes as exemplified recently by the 3-year evolution towards virulence at the AvrLm1 locus in French populations. The AvrLm1 avirulence gene was recently cloned and shown to be a solo gene within a 269 kb non-coding, heterochromatin-like region. Here we describe the sequencing of the AvrLm1 genomic region in one avirulent and two virulent isolates to investigate the molecular basis of evolution towards virulence at the AvrLm1 locus. For these virulent isolates, the gain of virulence was linked to a 260 kb deletion of a chromosomal segment spanning AvrLm1 and deletion breakpoints were identical or similar. Among the 460 isolates analysed from France, Australia and Mexico, a similar large deletion was apparent in > 90% of the virulent isolates. Deletion breakpoints were also strongly conserved in most of the virulent isolates, which led to the hypothesis that a unique deletion event leading to the avrLm1 virulence has diffused in pathogen populations. These data finally suggest that retrotransposons are key drivers in genome evolution and adaptation to novel selection pressure in L. maculans.

Complete genome sequence of the entomopathogenic and metabolically versatile soil bacterium Pseudomonas entomophila.

Pseudomonas entomophila is an entomopathogenic bacterium that, upon ingestion, kills Drosophila melanogaster as well as insects from different orders. The complete sequence of the 5.9-Mb genome was determined and compared to the sequenced genomes of four Pseudomonas species. P. entomophila possesses most of the catabolic genes of the closely related strain P. putida KT2440, revealing its metabolically versatile properties and its soil lifestyle. Several features that probably contribute to its entomopathogenic properties were disclosed. Unexpectedly for an animal pathogen, P. entomophila is devoid of a type III secretion system and associated toxins but rather relies on a number of potential virulence factors such as insecticidal toxins, proteases, putative hemolysins, hydrogen cyanide and novel secondary metabolites to infect and kill insects. Genome-wide random mutagenesis revealed the major role of the two-component system GacS/GacA that regulates most of the potential virulence factors identified.

Lost in the middle of nowhere: the AvrLm1 avirulence gene of the Dothideomycete Leptosphaeria maculans.

Leptosphaeria maculans, a Dothideomycete causing stem canker on oilseed rape (Brassica napus), develops gene-for-gene interactions with its host plants. To date, nine resistance genes (Rlm1-9) have been identified in Brassica spp. The corresponding nine avirulence genes (AvrLm1-9) in L. maculans have been mapped at four independent loci, thereby revealing two clusters of three and four linked avirulence genes. Here, we report the completion of map-based cloning of AvrLm1. AvrLm1 was genetically delineated within a 7.3 centimorgan interval corresponding to a 439 kb BAC contig. AvrLm1 is a single copy gene isolated within a 269 kb non-coding, heterochromatin-like region. The region comprised a number of degenerated, nested copies of four long-terminal repeat (LTR) retrotransposons, including Pholy and three novel Gypsy-like retrotransposons. AvrLm1 restored the avirulent phenotype on Rlm1 cultivars following functional complementation of virulent isolates. AvrLm1 homologues were not detected in other Leptosphaeria species or in known fungal genomes including the closely related species Stagonospora nodorum. The predicted AvrLm1 protein is composed of 205 amino acids, of which only one is a cysteine residue. It contains a peptide signal suggesting extracellular localization. Unlike most other fungal avirulence genes, AvrLm1 is constitutively expressed, with a probable increased level of expression upon plant infection, suggesting the absence of tight regulation of AvrLm1 expression.

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation.

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.

Molecular basis of evolutionary events that shaped the hardness locus in diploid and polyploid wheat species (Triticum and Aegilops).

The Hardness (Ha) locus controls grain hardness in hexaploid wheat (Triticum aestivum) and its relatives (Triticum and Aegilops species) and represents a classical example of a trait whose variation arose from gene loss after polyploidization. In this study, we investigated the molecular basis of the evolutionary events observed at this locus by comparing corresponding sequences of diploid, tertraploid, and hexaploid wheat species (Triticum and Aegilops). Genomic rearrangements, such as transposable element insertions, genomic deletions, duplications, and inversions, were shown to constitute the major differences when the same genomes (i.e., the A, B, or D genomes) were compared between species of different ploidy levels. The comparative analysis allowed us to determine the extent and sequences of the rearranged regions as well as rearrangement breakpoints and sequence motifs at their boundaries, which suggest rearrangement by illegitimate recombination. Among these genomic rearrangements, the previously reported Pina and Pinb genes loss from the Ha locus of polyploid wheat species was caused by a large genomic deletion that probably occurred independently in the A and B genomes. Moreover, the Ha locus in the D genome of hexaploid wheat (T. aestivum) is 29 kb smaller than in the D genome of its diploid progenitor Ae. tauschii, principally because of transposable element insertions and two large deletions caused by illegitimate recombination. Our data suggest that illegitimate DNA recombination, leading to various genomic rearrangements, constitutes one of the major evolutionary mechanisms in wheat species.

Truncated and RIP-degenerated copies of the LTR retrotransposon Pholy are clustered in a pericentromeric region of the Leptosphaeria maculans genome.

The LMR1 5.2 kb interspersed repeat of Leptosphaeria maculans was described by Taylor and Borgmann [Mol. Plant Microbe Interact. 7 (1994) 181] as an uncharacterized repeated element sharing homologies with both LINEs and SINEs. Here, we used the LMR1 sequence as a template to identify the full-length element within a 184-kb genomic sequence corresponding to the pericentromeric region of the 2.80 Mb chromosome of isolate v23.1.3. This region comprises (i) one 6980-bp full-sized Pholy element bordered by two 275- to 280-bp long terminal repeats (LTRs), (ii) five Pholy-related sequences, usually truncated at their 3' ends, and (iii) five solo-LTRs. Structural features strongly suggested that Pholy corresponds to an ancient copia-like retrotransposon, sharing strong homologies with the Elsa retrotransposon of Stagonospora nodorum. Pholy was also suggested to be specific to pericentromeric regions. Comparative analysis of the structure of the Pholy-like sequences occurring in the 184-kb contig and in other parts of the genome showed that this family of repeats is highly degenerated following extensive repeat induced point mutation (RIP).

Bracoviruses contain a large multigene family coding for protein tyrosine phosphatases.

The relationship between parasitic wasps and bracoviruses constitutes one of the few known mutualisms between viruses and eukaryotes. The virions produced in the wasp ovaries are injected into host lepidopteran larvae, where virus genes are expressed, allowing successful development of the parasite by inducing host immune suppression and developmental arrest. Bracovirus-bearing wasps have a common phylogenetic origin, and contemporary bracoviruses are hypothesized to have been inherited by chromosomal transmission from a virus that originally integrated into the genome of the common ancestor wasp living 73.7 +/- 10 million years ago. However, so far no conserved genes have been described among different braconid wasp subfamilies. Here we show that a gene family is present in bracoviruses of different braconid wasp subfamilies (Cotesia congregata, Microgastrinae, and Toxoneuron nigriceps, Cardiochilinae) which likely corresponds to an ancient component of the bracovirus genome that might have been present in the ancestral virus. The genes encode proteins belonging to the protein tyrosine phosphatase family, known to play a key role in the control of signal transduction pathways. Bracovirus protein tyrosine phosphatase genes were shown to be expressed in different tissues of parasitized hosts, and two protein tyrosine phosphatases were produced with recombinant baculoviruses and tested for their biochemical activity. One protein tyrosine phosphatase is a functional phosphatase. These results strengthen the hypothesis that protein tyrosine phosphatases are involved in virally induced alterations of host physiology during parasitism.

Genome sequence of a polydnavirus: insights into symbiotic virus evolution.

Little is known of the fate of viruses involved in long-term obligatory associations with eukaryotes. For example, many species of parasitoid wasps have symbiotic viruses to manipulate host defenses and to allow development of parasitoid larvae. The complete nucleotide sequence of the DNA enclosed in the virus particles injected by a parasitoid wasp revealed a complex organization, resembling a eukaryote genomic region more than a viral genome. Although endocellular symbiont genomes have undergone a dramatic loss of genes, the evolution of symbiotic viruses appears to be characterized by extensive duplication of virulence genes coding for truncated versions of cellular proteins.

Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype.

Tetraodon nigroviridis is a freshwater puffer fish with the smallest known vertebrate genome. Here, we report a draft genome sequence with long-range linkage and substantial anchoring to the 21 Tetraodon chromosomes. Genome analysis provides a greatly improved fish gene catalogue, including identifying key genes previously thought to be absent in fish. Comparison with other vertebrates and a urochordate indicates that fish proteins have diverged markedly faster than their mammalian homologues. Comparison with the human genome suggests approximately 900 previously unannotated human genes. Analysis of the Tetraodon and human genomes shows that whole-genome duplication occurred in the teleost fish lineage, subsequent to its divergence from mammals. The analysis also makes it possible to infer the basic structure of the ancestral bony vertebrate genome, which was composed of 12 chromosomes, and to reconstruct much of the evolutionary history of ancient and recent chromosome rearrangements leading to the modern human karyotype.

Genome evolution in yeasts.

Identifying the mechanisms of eukaryotic genome evolution by comparative genomics is often complicated by the multiplicity of events that have taken place throughout the history of individual lineages, leaving only distorted and superimposed traces in the genome of each living organism. The hemiascomycete yeasts, with their compact genomes, similar lifestyle and distinct sexual and physiological properties, provide a unique opportunity to explore such mechanisms. We present here the complete, assembled genome sequences of four yeast species, selected to represent a broad evolutionary range within a single eukaryotic phylum, that after analysis proved to be molecularly as diverse as the entire phylum of chordates. A total of approximately 24,200 novel genes were identified, the translation products of which were classified together with Saccharomyces cerevisiae proteins into about 4,700 families, forming the basis for interspecific comparisons. Analysis of chromosome maps and genome redundancies reveal that the different yeast lineages have evolved through a marked interplay between several distinct molecular mechanisms, including tandem gene repeat formation, segmental duplication, a massive genome duplication and extensive gene loss.

Characterization of the genomic organization of the region bordering the centromere of chromosome V of Podospora anserina by direct sequencing.

A Podospora anserina BAC library of 4800 clones has been constructed in the vector pBHYG allowing direct selection in fungi. Screening of the BAC collection for centromeric sequences of chromosome V allowed the recovery of clones localized on either sides of the centromere, but no BAC clone was found to contain the centromere. Seven BAC clones containing 322,195 and 156,244bp from either sides of the centromeric region were sequenced and annotated. One 5S rRNA gene, 5 tRNA genes, and 163 putative coding sequences (CDS) were identified. Among these, only six CDS seem specific to P. anserina. The gene density in the centromeric region is approximately one gene every 2.8kb. Extrapolation of this gene density to the whole genome of P. anserina suggests that the genome contains about 11,000 genes. Synteny analyses between P. anserina and Neurospora crassa show that co-linearity extends at the most to a few genes, suggesting rapid genome rearrangements between these two species.

Identification of the fertility restoration locus, Rfo, in radish, as a member of the pentatricopeptide-repeat protein family.

Ogura cytoplasmic male sterility (CMS) in radish (Raphanus sativus) is caused by an aberrant mitochondrial gene, Orf138, that prevents the production of functional pollen without affecting female fertility. Rfo, a nuclear gene that restores male fertility, alters the expression of Orf138 at the post-transcriptional level. The Ogura CMS/Rfo two-component system is a useful model for investigating nuclear-cytoplasmic interactions, as well as the physiological basis of fertility restoration. Using a combination of positional cloning and microsynteny analysis of Arabidopsis thaliana and radish, we genetically and physically delimited the Rfo locus to a 15-kb DNA segment. Analysis of this segment shows that Rfo is a member of the pentatricopeptide repeat (PPR) family. In Arabidopsis, this family contains more than 450 members of unknown function, although most of them are predicted to be targeted to mitochondria and chloroplasts and are thought to have roles in organellar gene expression.

The DNA sequence and analysis of human chromosome 14.

Chromosome 14 is one of five acrocentric chromosomes in the human genome. These chromosomes are characterized by a heterochromatic short arm that contains essentially ribosomal RNA genes, and a euchromatic long arm in which most, if not all, of the protein-coding genes are located. The finished sequence of human chromosome 14 comprises 87,410,661 base pairs, representing 100% of its euchromatic portion, in a single continuous segment covering the entire long arm with no gaps. Two loci of crucial importance for the immune system, as well as more than 60 disease genes, have been localized so far on chromosome 14. We identified 1,050 genes and gene fragments, and 393 pseudogenes. On the basis of comparisons with other vertebrate genomes, we estimate that more than 96% of the chromosome 14 genes have been annotated. From an analysis of the CpG island occurrences, we estimate that 70% of these annotated genes are complete at their 5' end.

Analysis of molecular markers genetically linked to the Leptosphaeria maculans avirulence gene AvrLm1 in field populations indicates a highly conserved event leading to virulence on Rlm1 genotypes.

Map-based cloning of the avirulence gene AvrLm1 of Leptosphaeria maculans was initiated utilizing a genetic map of the fungus and a BAC library constructed from an AvrLm1 isolate. Seven polymorphic DNA markers closely linked to AvrLm1 were identified. Of these, two were shown to border the locus on its 5' end and were present, with size polymorphism, in both the virulent and the avirulent isolates. In contrast, three markers, J19-1.1, J53-1.3 (in coupling phase with avirulence), and Vir1 (in repulsion phase with avirulence), cosegregated with AvrLm1 in 312 progeny from five in vitro crosses. J19-1.1 and J53-1.3 were never amplified in the virulent parents or progeny, whereas Vir1 was never amplified in the avirulent parents or progeny. J19-1.1 and J53-1.3 were shown to be separated by 40 kb within a 184-kb BAC contig. In addition, the 1.6-cM genetic distance between J53-1.3 and the nearest recombinant marker corresponded to a 121-kb physical distance. When analyzing a European Union-wide collection of 192 isolates, J53-1.3, J19-1.1, and Vir1 were found to be closely associated with the AvrLm1 locus. The results of polymerase chain reaction amplification with primers for the three markers were in accordance with the interaction phenotype for 92.2% (J53-1.3), 90.6% (J19-1.1), and 88.0% (Vir1) of the isolates. In addition, genome organization of the AvrLm1 region was highly conserved in field isolates, because 89.1% of the avirulent isolates and 79.0% of the virulent isolates showed the same association of markers as that of the parents of in vitro crosses. The large-scale analysis of field isolates with markers originating from the genetic map therefore confirms (i) the physical proximity between the markers and the target locus and (ii) that AvrLm1 is located in (or close to) a recombination-deficient genome region. As a consequence, map-based markers provided us with high-quality markers for an overview of the occurrence of race "AvrLm1" at the field scale. These data were used to propose hypotheses on evolution towards virulence in field isolates.

Comparative sequence analysis of the X-inactivation center region in mouse, human, and bovine.

We have sequenced to high levels of accuracy 714-kb and 233-kb regions of the mouse and bovine X-inactivation centers (Xic), respectively, centered on the Xist gene. This has provided the basis for a fully annotated comparative analysis of the mouse Xic with the 2.3-Mb orthologous region in human and has allowed a three-way species comparison of the core central region, including the Xist gene. These comparisons have revealed conserved genes, both coding and noncoding, conserved CpG islands and, more surprisingly, conserved pseudogenes. The distribution of repeated elements, especially LINE repeats, in the mouse Xic region when compared to the rest of the genome does not support the hypothesis of a role for these repeat elements in the spreading of X inactivation. Interestingly, an asymmetric distribution of LINE elements on the two DNA strands was observed in the three species, not only within introns but also in intergenic regions. This feature is suggestive of important transcriptional activity within these intergenic regions. In silico prediction followed by experimental analysis has allowed four new genes, Cnbp2, Ftx, Jpx, and Ppnx, to be identified and novel, widespread, complex, and apparently noncoding transcriptional activity to be characterized in a region 5' of Xist that was recently shown to attract histone modification early after the onset of X inactivation.

Molecular analysis of nonrandom 8q12 deletions in acute lymphoblastic leukemia: identification of two candidate genes.

Acute lymphoblastic leukemia is the most common malignancy in childhood. High-resolution allelotyping performed in our laboratory showed new chromosomal sites of nonrandom deletions. We have focused our work on 8q12 deletions, which we have found in about 4% of patients (eight of 205 informative cases). These deletions were of small size (less than 1 Mb) in all but one patient, and the deleted region common to all patients was delineated between two microsatellite markers (D8S1113 and D8S1763). This region was sequenced entirely from two overlapping bacterial artificial chromosomes. The common deleted region (120 kb) had a low GC content (37%), was composed more than 50% of LINE sequences, and contained only two candidate genes. The centromeric deletion borders were clustered within an interval of 33 kb between two microsatellite markers. This interval contains the first exon of an HMG-1-related gene (KIAA0808) and a putative gene, DL8q12, predicted to encode a protein with 231 amino acid residues with no homolog in protein databases. Analysis of the available mRNA from lymphoblastic cells of two patients with 8q12 deletions using common polymorphisms in the 3' UTR of KIAA0808 showed monoallelic expression of this gene. Identification of a biallelic polymorphism in the first exon of DL8q12 showed that this gene was deleted in two of four informative cases. Sequencing of the exons of both genes from all patients with 8q12 deletions did not show any mutation, which suggests that neither of these genes behaves as a classic tumor suppressor gene.