Search

SUMMARY Next generation sequencing (NGS) technologies, which allow high speed automated DNA sequencing, DNA sequence analysis and comparisons using big data algorithms are being used more and more in medical diagnosis, to provide prognosis information and for choosing a treatment that best fits the patient. The French “Academie nationale de medecine ” (ANM) and “ Academie des technologies “ (NATF) jointly stress the scientific and medical importance of these technologies and call for bringing together strengths of the national computer science industry and the medical community in a public-private consortium in order to build a Demonstrator and Center network filling up the growing gap between France and the most advanced countries in this field with great potential. They underline also common pitfalls which need to be addressed in technic, legal, economic, education and ethics issues, in order to make this complex project a success.

Balanced chromosomal rearrangements (BCR) are associated with abnormal phenotypes in approximately 6% of balanced translocations and 9.4% of balanced inversions. Abnormal phenotypes can be caused by disruption of genes at the breakpoints, deletions, or positional effects. Conventional cytogenetic techniques have a limited resolution and do not enable a thorough genetic investigation. Molecular techniques applied to BCR carriers can contribute to the characterization of this type of chromosomal rearrangement and to the phenotype-genotype correlation. Fifteen individuals among 35 with abnormal phenotypes and BCR were selected for further investigation by molecular techniques. Chromosomal rearrangements involved 11 reciprocal translocations, 3 inversions, and 1 balanced insertion. Array genomic hybridization (AGH) was performed and genomic imbalances were detected in 20% of the cases, 1 at a rearrangement breakpoint and 2 further breakpoints in other chromosomes. Alterations were further confirmed by FISH and associated with the phenotype of the carriers. In the analyzed cases not showing genomic imbalances by AGH, next-generation sequencing (NGS), using whole genome libraries, prepared following the Illumina TruSeq DNA PCR-Free protocol (Illumina®) and then sequenced on an Illumina HiSEQ 2000 as 150-bp paired-end reads, was done. The NGS results suggested breakpoints in 7 cases that were similar or near those estimated by karyotyping. The genes overlapping 6 breakpoint regions were analyzed. Follow-up of BCR carriers would improve the knowledge about these chromosomal rearrangements and their consequences.

MetA and MetX are phylogenetically unrelated families of acyl-L-homoserine transferases. Experimental assignation of function and structural modeling of these families correct widespread misannotation, reveal convergence of function and uncover new functions in a subclass of MetX. Experimental validation of enzyme function is crucial for genome interpretation, but it remains challenging because it cannot be scaled up to accommodate the constant accumulation of genome sequences. We tackled this issue for the MetA and MetX enzyme families, phylogenetically unrelated families of acyl-L-homoserine transferases involved in L-methionine biosynthesis. Members of these families are prone to incorrect annotation because MetX and MetA enzymes are assumed to always use acetyl-CoA and succinyl-CoA, respectively. We determined the enzymatic activities of 100 enzymes from diverse species, and interpreted the results by structural classification of active sites based on protein structure modeling. We predict that >60% of the 10,000 sequences from these families currently present in databases are incorrectly annotated, and suggest that acetyl-CoA was originally the sole substrate of these isofunctional enzymes, which evolved to use exclusively succinyl-CoA in the most recent bacteria. We also uncovered a divergent subgroup of MetX enzymes in fungi that participate only in L-cysteine biosynthesis as O-succinyl-L-serine transferases.

The exponential increase in genome sequencing output has led to the accumulation of thousands of predicted genes lacking a proper functional annotation. Among this mass of hypothetical proteins, enzymes catalyzing new reactions or using novel ways to catalyze already known reactions might still wait to be identified. Here, we provide a structural and biochemical characterization of the 3-keto-5-aminohexanoate cleavage enzyme (Kce), an enzymatic activity long known as being involved in the anaerobic fermentation of lysine but whose catalytic mechanism has remained elusive so far. Although the enzyme shows the ubiquitous triose phosphate isomerase (TIM) barrel fold and a Zn2+ cation reminiscent of metal-dependent class II aldolases, our results based on a combination of x-ray snapshots and molecular modeling point to an unprecedented mechanism that proceeds through deprotonation of the 3-keto-5-aminohexanoate substrate, nucleophilic addition onto an incoming acetyl-CoA, intramolecular transfer of the CoA moiety, and final retro-Claisen reaction leading to acetoacetate and 3-aminobutyryl-CoA. This model also accounts for earlier observations showing the origin of carbon atoms in the products, as well as the absence of detection of any covalent acyl-enzyme intermediate. Kce is the first representative of a large family of prokaryotic hypothetical proteins, currently annotated as the “domain of unknown function” DUF849.

Prochlorococcus marinus , the dominant photosynthetic organism in the ocean, is found in two main ecological forms: high-light-adapted genotypes in the upper part of the water column and low-light-adapted genotypes at the bottom of the illuminated layer. P. marinus SS120, the complete genome sequence reported here, is an extremely low-light-adapted form. The genome of P. marinus SS120 is composed of a single circular chromosome of 1,751,080 bp with an average G+C content of 36.4%. It contains 1,884 predicted protein-coding genes with an average size of 825 bp, a single rRNA operon, and 40 tRNA genes. Together with the 1.66-Mbp genome of P. marinus MED4, the genome of P. marinus SS120 is one of the two smallest genomes of a photosynthetic organism known to date. It lacks many genes that are involved in photosynthesis, DNA repair, solute uptake, intermediary metabolism, motility, phototaxis, and other functions that are conserved among other cyanobacteria. Systems of signal transduction and environmental stress response show a particularly drastic reduction in the number of components, even taking into account the small size of the SS120 genome. In contrast, housekeeping genes, which encode enzymes of amino acid, nucleotide, cofactor, and cell wall biosynthesis, are all present. Because of its remarkable compactness, the genome of P. marinus SS120 might approximate the minimal gene complement of a photosynthetic organism.

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.

Aldosterone and vasopressin are responsible for the final adjustment of sodium and water reabsorption in the kidney. In principal cells of the kidney cortical collecting duct (CCD), the integral response to aldosterone and the long-term functional effects of vasopressin depend on transcription. In this study, we analyzed the transcriptome of a highly differentiated mouse clonal CCD principal cell line (mpkCCD cl4 ) and the changes in the transcriptome induced by aldosterone and vasopressin. Serial analysis of gene expression (SAGE) was performed on untreated cells and on cells treated with either aldosterone or vasopressin for 4 h. The transcriptomes in these three experimental conditions were determined by sequencing 169,721 transcript tags from the corresponding SAGE libraries. Limiting the analysis to tags that occurred twice or more in the data set, 14,654 different transcripts were identified, 3,642 of which do not match known mouse sequences. Statistical comparison (at P < 0.05 level) of the three SAGE libraries revealed 34 AITs (aldosterone-induced transcripts), 29 ARTs (aldosterone-repressed transcripts), 48 VITs (vasopressin-induced transcripts) and 11 VRTs (vasopressin-repressed transcripts). A selection of the differentially-expressed, hormone-specific transcripts (5 VITs, 2 AITs and 1 ART) has been validated in the mpkCCD cl4 cell line either by Northern blot hybridization or reverse transcription–PCR . The hepatocyte nuclear transcription factor HNF-3-α (VIT39), the receptor activity modifying protein RAMP3 (VIT48), and the glucocorticoid-induced leucine zipper protein (GILZ) (AIT28) are candidate proteins playing a role in physiological responses of this cell line to vasopressin and aldosterone.

This paper reports the genomic analysis of strain CBS732 of Zygosaccharomyces rouxii, a homothallic diploid yeast. We explored the sequences of 4934 random sequencing tags of about 1 kb in size and compared them to the Saccharomyces cerevisiae gene products. Approximately 2250 nuclear genes, 57 tRNAs, the rDNA locus, the endogenous pSR1 plasmid and 15 mitochondrial genes were identified. According to 18S and 25S rRNA cladograms and to synteny analysis, Z. rouxii could be placed among the S. cerevisiae sensu lato yeasts.

As part of the comparative genomics project ‘GENOLEVURES’, we studied the Kluyveromyces marxianus var. marxianus strain CBS712 using a partial random sequencing strategy. With a 0.2×genome equivalent coverage, we identified ca. 1300 novel genes encoding proteins, some containing spliceosomal introns with consensus splice sites identical to those of Saccharomyces cerevisiae, 28 tRNA genes, the whole rDNA repeat, and retrotransposons of the Ty1/2 family of S. cerevisiae with diverged Long Terminal Repeats. Functional classification of the K. marxianus genes, as well as the analysis of the paralogous gene families revealed few differences with respect to S. cerevisiae. Only 42 K. marxianus identified genes are without detectable homolog in the baker’s yeast. However, we identified several genetic rearrangements between these two yeast species.

We have analyzed the evolution of chromosome maps of Hemiascomycetes by comparing gene order and orientation of the 13 yeast species partially sequenced in this program with the genome map of Saccharomyces cerevisiae. From the analysis of nearly 8000 situations in which two distinct genes having homologs in S. cerevisiae could be identified on the sequenced inserts of another yeast species, we have quantified the loss of synteny, the frequency of single gene deletion and the occurrence of gene inversion. Traces of ancestral duplications in the genome of S. cerevisiae could be identified from the comparison with the other species that do not entirely coincide with those identified from the comparison of S. cerevisiae with itself. From such duplications and from the correlation observed between gene inversion and loss of synteny, a model is proposed for the molecular evolution of Hemiascomycetes. This model, which can possibly be extended to other eukaryotes, is based on the reiteration of events of duplication of chromosome segments, creating transient merodiploids that are subsequently resolved by single gene deletion events.

The generation of sequencing data for the hemiascomycetous yeast random sequence tag project was performed using the procedures established at GENOSCOPE. These procedures include a series of protocols for the sequencing reactions, using infra-red labelled primers, performed on both ends of the plasmid inserts in the same reaction tube, and their analysis on automated DNA sequencers. They also include a package of computer programs aimed at detecting potential assignation errors, selecting good quality sequences and estimating their useful length.

The identification of molecular evolutionary mechanisms in eukaryotes is approached by a comparative genomics study of a homogeneous group of species classified as Hemiascomycetes. This group includes Saccharomyces cerevisiae, the first eukaryotic genome entirely sequenced, back in 1996. A random sequencing analysis has been performed on 13 different species sharing a small genome size and a low frequency of introns. Detailed information is provided in the 20 following papers. Additional tables available on websites describe the ca. 20 000 newly identified genes. This wealth of data, so far unique among eukaryotes, allowed us to examine the conservation of chromosome maps, to identify the ‘yeast-specific’ genes, and to review the distribution of gene families into functional classes. This project conducted by a network of seven French laboratories has been designated ‘Genolevures’.

Random sequencing of the Kluyveromyces lactis genome allowed the identification of 2235–2601 open reading frames (ORFs) homologous to S. cerevisiae ORFs, 51 ORFs which were homologous to genes from other species, 64 tRNAs, the complete rDNA repeat, and a few Ty1- and Ty2-like sequences. In addition, the complete sequence of plasmid pKD1 and a large coverage of the mitochondrial genome were obtained. The global distribution into general functional categories found in Saccharomyces cerevisiae and as defined by MIPS is well conserved in K. lactis. However, detailed examination of certain subcategories revealed a small excess of genes involved in amino acid metabolism in K. lactis. The sequences are deposited at EMBL under the accession numbers AL424881–AL430960.

Autosomal dominant hereditary spastic paraplegia (AD-HSP) is a genetically heterogeneous neurodegenerative disorder characterized by progressive spasticity of the lower limbs. Among the four loci causing AD-HSP identified so far, the SPG4 locus at chromosome 2p2-1p22 has been shown to account for 40-50% of all AD-HSP families. Using a positional cloning strategy based on obtaining sequence of the entire SPG4 interval, we identified a candidate gene encoding a new member of the AAA protein family, which we named spastin. Sequence analysis of this gene in seven SPG4-linked pedigrees revealed several DNA modifications, including missense, nonsense and splice-site mutations. Both SPG4 and its mouse orthologue were shown to be expressed early and ubiquitously in fetal and adult tissues. The sequence homologies and putative subcellular localization of spastin suggest that this ATPase is involved in the assembly or function of nuclear protein complexes.

We performed high-throughput sequencing of DNA from fossilized faeces to evaluate this material as a source of information on the genome and diet of Pleistocene carnivores. We analysed coprolites derived from the extinct cave hyena ( Crocuta crocuta spelaea ), and sequenced 90 million DNA fragments from two specimens. The DNA reads enabled a reconstruction of the cave hyena mitochondrial genome with up to a 158-fold coverage. This genome, and those sequenced from extant spotted ( Crocuta crocuta ) and striped ( Hyaena hyaena ) hyena specimens, allows for the establishment of a robust phylogeny that supports a close relationship between the cave and the spotted hyena. We also demonstrate that high-throughput sequencing yields data for cave hyena multi-copy and single-copy nuclear genes, and that about 50 per cent of the coprolite DNA can be ascribed to this species. Analysing the data for additional species to indicate the cave hyena diet, we retrieved abundant sequences for the red deer ( Cervus elaphus ), and characterized its mitochondrial genome with up to a 3.8-fold coverage. In conclusion, we have demonstrated the presence of abundant ancient DNA in the coprolites surveyed. Shotgun sequencing of this material yielded a wealth of DNA sequences for a Pleistocene carnivore and allowed unbiased identification of diet.

Motivation: Current computational approaches to function prediction are mostly based on protein sequence classification and transfer of annotation from known proteins to their closest homologous sequences relying on the orthology concept of function conservation. This approach suffers a major weakness: annotation reliability depends on global sequence similarity to known proteins and is poorly efficient for enzyme superfamilies that catalyze different reactions. Structural biology offers a different strategy to overcome the problem of annotation by adding information about protein 3D structures. This information can be used to identify amino acids located in active sites, focusing on detection of functional polymorphisms residues in an enzyme superfamily. Structural genomics programs are providing more and more novel protein structures at a high-throughput rate. However, there is still a huge gap between the number of sequences and available structures. Computational methods, such as homology modeling provides reliable approaches to bridge this gap and could be a new precise tool to annotate protein functions. Results: Here, we present Active Sites Modeling and Clustering (ASMC) method, a novel unsupervised method to classify sequences using structural information of protein pockets. ASMC combines homology modeling of family members, structural alignment of modeled active sites and a subsequent hierarchical conceptual classification. Comparison of profiles obtained from computed clusters allows the identification of residues correlated to subfamily function divergence, called specificity determining positions. ASMC method has been validated on a benchmark of 42 Pfam families for which previous resolved holo-structures were available. ASMC was also applied to several families containing known protein structures and comprehensive functional annotations. We will discuss how ASMC improves annotation and understanding of protein families functions by giving some specific illustrative examples on nucleotidyl cyclases, protein kinases and serine proteases. Availability: http://www.genoscope.fr/ASMC/. Contact: raquelcm@dcc.ufmg.br; kbastard@genoscope.cns.fr; artigue@genoscope.cns.fr Supplementary information: Supplementary data are available at Bioinformatics online.

The Genoscope annotation workflow for eukaryote genomes relies on evidence from ab initio gene models predictions combined with homology searches, using collections of expressed sequences - full length cDNAs, ESTs or massive-scale mRNA sequences from the same or closely related organisms – proteins or other genomic sequences. Global analysis of these drafts or complete sequences are then combining both approaches in the form of gene prediction data integration using GAZE, capable to identify a majority of the existing gene features. Although of very good quality, gene-modelling remains still tentative at the end of the process. Even though computational predictors are useful on large scale annotation for global genomics analysis, there is no complete genome for which all gene structures, in terms of exons, introns and coding regions, have been experimentally confirmed.Finished genomes can provide exciting insights into the genome organization and evolution. Additional experimental data generated by genome sequencing projects give assistance to genome annotation aiming to a better understanding of the biology of the organism. Therefore, gene models and annotation can be improved by human curation to find errors or to resolve incongruous evidence on the automatic annotation of the genome. We now provide to collaborators carrying sequencing projects with a distributed annotation platform allowing expert evaluation of the annotation, in addition to our automated gene prediction pipeline.To ensure at most the participation of the scientific community, an annotation tool for revising annotations has been set up using components of the Generic Model Organism Database toolkit, which provides tools for managing organism databases. A CHADO database, linked to an Apollo graphical interface, permit users to correct gene structures and store them in a dedicated organism database, as we will show on a few examples. Such a tool would facilitate connecting and comparing predicted annotations with existing biological data, becoming the repository of complete annotated finished genome sequence.

Background Massively parallel DNA sequencing instruments are enabling the decoding of whole genomes at significantly lower cost and higher throughput than classical Sanger technology. Each of these technologies have been estimated to yield assemblies with more problematic features than the standard method. These problems are of a different nature depending on the techniques used. So, an appropriate mix of technologies may help resolve most difficulties, and eventually provide assemblies of high quality without requiring any Sanger-based input. Results We compared assemblies obtained using Sanger data with those from different inputs from New Sequencing Technologies. The assemblies were systematically compared with a reference finished sequence. We found that the 454 GSFLX can efficiently produce high continuity when used at high coverage. The potential to enhance continuity by scaffolding was tested using 454 sequences from circularized genomic fragments. Finally, we explore the use of Solexa-Illumina short reads to polish the genome draft by implementing a technique to correct 454 consensus errors. Conclusion High quality drafts can be produced for small genomes without any Sanger data input. We found that 454 GSFLX and Solexa/Illumina show great complementarity in producing large contigs and supercontigs with a low error rate.

In most bilaterian organisms so far studied, Hox genes are organized in genomic clusters and determine development along the anteroposterior axis. It has been suggested that this clustering, together with spatial and temporal colinearity of gene expression, represents the ancestral condition. However, in organisms with derived modes of embryogenesis and lineage-dependent mechanisms for the determination of cell fate, temporal colinearity of expression can be lost and Hox cluster organization disrupted, as is the case for the ecdysozoans Drosophila melanogaster and Caenorhabditis elegans and the urochordates Ciona intestinalis and Oikopleura dioica. We sought to determine whether a lophotrochozoan, the platyhelminth parasite Schistosoma mansoni, possesses a conserved or disrupted Hox cluster. Using a polymerase chain reaction (PCR)-based strategy, we have cloned and characterized three novel S. mansoni genes encoding orthologues of Drosophila labial (SmHox 1), deformed (SmHox4), and abdominal A (SmHox8), as well as the full-length coding sequence of the previously described Smox1, which we identify as an orthologue of fushi tarazu. Quantitative reverse transcriptase-PCR showed that the four genes were expressed at all life-cycle stages but that levels of expression were differentially regulated. Phylogenetic analysis and the conservation of "parapeptide'' sequences ('"-terminal to the homeodomains of SmHox8 and Smox1 support the grouping of platyhelminths within the lophotrochozoan clade. However, Bacterial Artificial Chromosome (BAC) library screening followed by genome walking failed to reconstitute a cluster. The BAC clones containing Hox genes were sequenced. and in no case were other Hox genes found on the same clone. Moreover, the SmHox4 and SmHox8 genes contained single very large introns (>40 kbp) further indicating that the schistosome Hox cluster is highly extended. Localization of the Hox genes to chromosomes using fluorescence in situ hybridization showed that SinHox4 and SmHox8 are on the long arm of chromosome 4, whereas SmHox1 and Smox1 are on chromosome 3. In silico screening of the available genome sequences corroborated results of Southern blotting and BAC library screening that indicate that there are no paralogues of SmHox1, SmHox4, or SmHox8. The schistosome Hox cluster is therefore not duplicated, but is both dispersed and disintegrated in the genome.

We report the construction of a tiling path of around 650 clones covering more than 99% of human chromosome 14. Clone overlap information to assemble the map was derived by comparing fully sequenced clones with a database of clone end sequences (sequence tag connector strategy). We selected homogeneously distributed seed points using an auxiliary high-resolution radiation hybrid map comprising 1,895 distinct positions. The high long-range continuity and low redundancy of the tiling path indicates that the sequence tag connector approach compares favourably with alternative mapping strategies.

Saccharomyces bayanus var. uvarum investigated here is the species closest to Saccharomyces cerevisiae. Random sequence tags (RSTs) allowed us to identify homologues to 2789 open reading frames (ORFs) in S. cerevisiae, ORFs duplicated in S. uvarum but not in S. cerevisiae, centromeres, tRNAs, homologues of Ty1/2 and Ty4 retrotransposons, and a complete rDNA repeat. Only 13 RSTs seem to be homologous to sequences in other organisms but not in S. cerevisiae. As the synteny between the two species is very high, cases in which synteny is lost suggest special mechanisms of genome evolution. The corresponding RSTs revealed that S. uvarum can exist without any S. cerevisiae DNA introgression. Accession numbers are from AL397139 to AL402278 in the EMBL databank.

A genomic exploration of Kluyveromyces thermotolerans was performed by random sequence tag (RST) analysis. We sequenced 2653 RSTs corresponding to inserts sequenced from both ends. We performed a systematic comparison with a complete set of proteins from Saccharomyces cerevisiae, other completely sequenced genomes and SwissProt. We identified six mitochondrial genes and 1358–1496 nuclear genes by comparison with S. cerevisiae. In addition, 25 genes were identified by comparison with other organisms. This corresponds to about 24% of the estimated gene content of this organism. A lower level of conservation is observed with orthologues to genes of S. cerevisiae previously classified as orphans. Gene order was found to be conserved between S. cerevisiae and K. thermotolerans in 56.5% of studied cases.

A method is described for mutagenesis of Pseudomonas fluorescens strains by electroporation with the transposon delivery vector pUT/mini-Tn5 Km. The transposition process was shown to be optimal at 12.5 kV cm-1 for a pulse time (Bowen and Koslak, 1992) of about 4 ms. The Pseudomonas fluorescens L6.5 target strain exhibited maximal electrocompetence when harvested at the middle of the exponential growth phase. As many as 7.7 10(5) mutants per picomole of delivery vector (7.5 kb) could be obtained, and these kanamycin-resistant mutants were shown to have lost the pUT plasmid. By external calibration with plasmids of increasing size (from 11.5 to 60.1 kb), the efficiency of the transformation process was evaluated to be approximately 1.31 x 10(8) transformants per picomole of delivery vector. Efficiency of the transposition process was 0.58%. This rapid method was used to tag for the cloning three independent chromosomal loci responsible for the Alk+ phenotype of Pseudomonas fluorescens L6.5 strain.