Jacques Corbeil

Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species

BackgroundThe process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly.ResultsIn Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies.ConclusionsMany current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

Ray: Simultaneous Assembly of Reads from a Mix of High-Throughput Sequencing Technologies

An accurate genome sequence of a desired species is now a pre-requisite for genome research. An important step in obtaining a high-quality genome sequence is to correctly assemble short reads into longer sequences accurately representing contiguous genomic regions. Current sequencing technologies continue to offer increases in throughput, and corresponding reductions in cost and time. Unfortunately, the benefit of obtaining a large number of reads is complicated by sequencing errors, with different biases being observed with each platform. Although software are available to assemble reads for each individual system, no procedure has been proposed for high-quality simultaneous assembly based on reads from a mix of different technologies. In this paper, we describe a parallel short-read assembler, called Ray, which has been developed to assemble reads obtained from a combination of sequencing platforms. We compared its performance to other assemblers on simulated and real datasets. We used a combination of Roche/454 and Illumina reads to assemble three different genomes. We showed that mixing sequencing technologies systematically reduces the number of contigs and the number of errors. Because of its open nature, this new tool will hopefully serve as a basis to develop an assembler that can be of universal utilization (availability: http://deNovoAssembler.sf.Net/). For online Supplementary Material , see www.liebertonline.com.

Ray Meta: scalable de novo metagenome assembly and profiling.

AbstractaVoluminous parallel sequencing datasets, especially metagenomic experiments, require distributed computing for de novo assembly and taxonomic profiling. Ray Meta is a massively distributed metagenome assembler that is coupled with Ray Communities, which profiles microbiomes based on uniquely-colored k-mers. It can accurately assemble and profile a three billion read metagenomic experiment representing 1,000 bacterial genomes of uneven proportions in 15 hours with 1,024 processor cores, using only 1.5 GB per core. The software will facilitate the processing of large and complex datasets, and will help in generating biological insights for specific environments. Ray Meta is open source and available at http://denovoassembler.sf.net.

Bendamustine (Treanda) Displays a Distinct Pattern of Cytotoxicity and Unique Mechanistic Features Compared with Other Alkylating Agents

Purpose: Bendamustine has shown clinical activity in patients with disease refractory to conventional alkylator chemotherapy. The purpose of this study was to characterize the mechanisms of action of bendamustine and to compare it with structurally related compounds. Experimental Design: Bendamustine was profiled in the National Cancer Institute in vitro antitumor screen. Microarray-based gene expression profiling, real-time PCR, immunoblot, cell cycle, and functional DNA damage repair analyses were used to characterize response to bendamustine and compare it with chlorambucil and phosphoramide mustard. Results: Bendamustine displays a distinct pattern of activity unrelated to other DNA-alkylating agents. Its mechanisms of action include activation of DNA-damage stress response and apoptosis, inhibition of mitotic checkpoints, and induction of mitotic catastrophe. In addition, unlike other alkylators, bendamustine activates a base excision DNA repair pathway rather than an alkyltransferase DNA repair mechanism. Conclusion: These results suggest that bendamustine possesses mechanistic features that differentiate it from other alkylating agents and may contribute to its distinct clinical efficacy profile.

Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia

Microarray techniques hold great promise for identifying risk factors for schizophrenia (SZ) but have not yet generated widely reproducible results due to methodological differences between studies and the high risk of type I inferential errors. Here we established a protocol for conservative analysis and interpretation of gene expression data from the dorsolateral prefrontal cortex of SZ patients using statistical and bioinformatic methods that limit false positives. We also compared brain gene expression profiles with those from peripheral blood cells of a separate sample of SZ patients to identify disease-associated genes that generalize across tissues and populations and further substantiate the use of gene expression profiling of blood for detecting valid SZ biomarkers. Implementing this systematic approach, we: (i) discovered 177 putative SZ risk genes in brain, 28 of which map to linked chromosomal loci; (ii) delineated six biological processes and 12 molecular functions that may be particularly disrupted in the illness; (iii) identified 123 putative SZ biomarkers in blood, 6 of which (BTG1, GSK3A, HLA-DRB1, HNRPA3, SELENBP1, and SFRS1) had corresponding differential expression in brain; (iv) verified the differential expression of the strongest candidate SZ biomarker (SELENBP1) in blood; and (v) demonstrated neuronal and glial expression of SELENBP1 protein in brain. The continued application of this approach in other brain regions and populations should facilitate the discovery of highly reliable and reproducible candidate risk genes and biomarkers for SZ. The identification of valid peripheral biomarkers for SZ may ultimately facilitate early identification, intervention, and prevention efforts as well.

A new perspective on V3 phenotype prediction.

The particular coreceptor used by a strain of HIV-1 to enter a host cell is highly indicative of its pathology. HIV-1 coreceptor usage is primarily determined by the amino add sequences of the V3 loop region of the viral envelope glycoprotein. The canonical approach to sequence-based prediction of coreceptor usage was derived via statistical analysis of a less reliable and significantly smaller data set than is presently available. We aimed to produce a superior phenotypic classifier by applying modern machine learning (ML) techniques to the current database of V3 loop sequences with known phenotype. The trained classifiers along with the sequence data are available for public use at the supplementary website: http://genomiac2.ucsd.edu:8080/wetcat/v3.html and http://www.cs.waikato.ac.nz/ml/weka[corrected].

Modulation of gene expression in drug resistant Leishmania is associated with gene amplification, gene deletion and chromosome aneuploidy

BackgroundDrug resistance can be complex, and several mutations responsible for it can co-exist in a resistant cell. Transcriptional profiling is ideally suited for studying complex resistance genotypes and has the potential to lead to novel discoveries. We generated full genome 70-mer oligonucleotide microarrays for all protein coding genes of the human protozoan parasites Leishmania major and Leishmania infantum. These arrays were used to monitor gene expression in methotrexate resistant parasites.ResultsLeishmania is a eukaryotic organism with minimal control at the level of transcription initiation and few genes were differentially expressed without concomitant changes in DNA copy number. One exception was found in Leishmania major, where the expression of whole chromosomes was down-regulated. The microarrays highlighted several mechanisms by which the copy number of genes involved in resistance was altered; these include gene deletion, formation of extrachromosomal circular or linear amplicons, and the presence of supernumerary chromosomes. In the case of gene deletion or gene amplification, the rearrangements have occurred at the sites of repeated (direct or inverted) sequences. These repeats appear highly conserved in both species to facilitate the amplification of key genes during environmental changes. When direct or inverted repeats are absent in the vicinity of a gene conferring a selective advantage, Leishmania will resort to supernumerary chromosomes to increase the levels of a gene product.ConclusionAneuploidy has been suggested as an important cause of drug resistance in several organisms and additional studies should reveal the potential importance of this phenomenon in drug resistance in Leishmania.

HIV induces lymphocyte apoptosis by a p53-initiated, mitochondrial-mediated mechanism

HIV‐1 induces apoptosis and leads to CD4+ T‐lymphocyte depletion in humans. It is still unclear whether HIV‐1 kills infected cells directly or indirectly. To elucidate the mechanisms of HIV‐1–induced apoptosis, we infected human CD4+ T cells with HIV‐1. Enzymatic analysis with fluorometric substrates showed that caspase 2, 3, and 9 were activated in CD4+ T cells with peak levels 48 h after infection. Immunoblotting analysis confirmed the cleavage of pro‐caspase 3 and 9, and of specific caspase substrates. Release of cytochrome c and apoptosis‐inducing factor (AIF) from mitochondria was observed in HIV‐infected cells. The cytochrome c and AIF release preceded the reduction of the mitochondrial transmembrane potential and nuclear chromatin condensation. HIV infection led to phosphorylation of p53 at the Ser15 residue, detectable as early as 24 h after infection. The p53 phosphorylation was followed by increased mRNA and protein expression of p21, Bax, HDM2, and p53. Up‐regulation of surface FasL expression, accompanied by a down‐regulation of Fas‐associated proteins (FADD, DAXX, and RIP), was observed 72 h after infection. Our results suggest that HIV activates the p53 pathway, leading to cytochrome c and AIF release with ensuing caspase activation.

Gene expression modulation is associated with gene amplification, supernumerary chromosomes and chromosome loss in antimony-resistant Leishmania infantum.

Antimonials remain the first line drug against the protozoan parasite Leishmania but their efficacy is threatened by resistance. We carried out a RNA expression profiling analysis comparing an antimony-sensitive and -resistant (Sb2000.1) strain of Leishmania infantum using whole-genome 70-mer oligonucleotide microarrays. Several genes were differentially expressed between the two strains, several of which were found to be physically linked in the genome. MRPA, an ATP-binding cassette (ABC) gene known to be involved in antimony resistance, was overexpressed in the antimony-resistant mutant along with three other tandemly linked genes on chromosome 23. This four gene locus was flanked by 1.4 kb repeated sequences from which an extrachromosomal circular amplicon was generated in the resistant cells. Interestingly, gene expression modulation of entire chromosomes occurred in the antimony-resistant mutant. Southern blots analyses and comparative genomic hybridizations revealed that this was either due to the presence of supernumerary chromosomes or to the loss of one chromosome. Leishmania parasites with haploid chromosomes were viable. Changes in copy number for some of these chromosomes were confirmed in another antimony-resistant strain. Selection of a partial revertant line correlated antimomy resistance levels and the copy number of aneuploid chromosomes, suggesting a putative link between aneuploidy and drug resistance in Leishmania.

Metagenomic Analysis of Stress Genes in Microbial Mat Communities from Antarctica and the High Arctic

ABSTRACT Polar and alpine microbial communities experience a variety of environmental stresses, including perennial cold and freezing; however, knowledge of genomic responses to such conditions is still rudimentary. We analyzed the metagenomes of cyanobacterial mats from Arctic and Antarctic ice shelves, using high-throughput pyrosequencing to test the hypotheses that consortia from these extreme polar habitats were similar in terms of major phyla and subphyla and consequently in their potential responses to environmental stresses. Statistical comparisons of the protein-coding genes showed similarities between the mats from the two poles, with the majority of genes derived from Proteobacteria and Cyanobacteria; however, the relative proportions differed, with cyanobacterial genes more prevalent in the Antarctic mat metagenome. Other differences included a higher representation of Actinobacteria and Alphaproteobacteria in the Arctic metagenomes, which may reflect the greater access to diasporas from both adjacent ice-free lands and the open ocean. Genes coding for functional responses to environmental stress (exopolysaccharides, cold shock proteins, and membrane modifications) were found in all of the metagenomes. However, in keeping with the greater exposure of the Arctic to long-range pollutants, sequences assigned to copper homeostasis genes were statistically (30%) more abundant in the Arctic samples. In contrast, more reads matching the sigma B genes were identified in the Antarctic mat, likely reflecting the more severe osmotic stress during freeze-up of the Antarctic ponds. This study underscores the presence of diverse mechanisms of adaptation to cold and other stresses in polar mats, consistent with the proportional representation of major bacterial groups.