Single-molecule DNA-mapping and whole-genome sequencing of individual cells.

To elucidate cellular diversity and clonal evolution in tissues and tumors, one must resolve genomic heterogeneity in single cells. To this end, we have developed low-cost, mass-producible micro-/nanofluidic chips for DNA extraction from individual cells. These chips have modules that collect genomic DNA for sequencing or map genomic structure directly, on-chip, with denaturation-renaturation (D-R) optical mapping [Marie R, et al. (2013) Proc Natl Acad Sci USA 110:4893-4898]. Processing of single cells from the LS174T colorectal cancer cell line showed that D-R mapping of single molecules can reveal structural variation (SV) in the genome of single cells. In one experiment, we processed 17 fragments covering 19.8 Mb of the cell's genome. One megabase-large fragment aligned well to chromosome 19 with half its length, while the other half showed variable alignment. Paired-end single-cell sequencing supported this finding, revealing a region of complexity and a 50-kb deletion. Sequencing struggled, however, to detect a 20-kb gap that D-R mapping showed clearly in a megabase fragment that otherwise mapped well to the reference at the pericentromeric region of chromosome 4. Pericentromeric regions are complex and show substantial sequence homology between different chromosomes, making mapping of sequence reads ambiguous. Thus, D-R mapping directly, from a single molecule, revealed characteristics of the single-cell genome that were challenging for short-read sequencing.

Ultra-deep sequencing of foraminiferal microbarcodes unveils hidden richness of early monothalamous lineages in deep-sea sediments.

Deep-sea floors represent one of the largest and most complex ecosystems on Earth but remain essentially unexplored. The vastness and remoteness of this ecosystem make deep-sea sampling difficult, hampering traditional taxonomic observations and diversity assessment. This problem is particularly true in the case of the deep-sea meiofauna, which largely comprises small-sized, fragile, and difficult-to-identify metazoans and protists. Here, we introduce an ultra-deep sequencing-based metagenetic approach to examine the richness of benthic foraminifera, a principal component of deep-sea meiofauna. We used Illumina sequencing technology to assess foraminiferal richness in 31 unsieved deep-sea sediment samples from five distinct oceanic regions. We sequenced an extremely short fragment (36 bases) of the small subunit ribosomal DNA hypervariable region 37f, which has been shown to accurately distinguish foraminiferal species. In total, we obtained 495,978 unique sequences that were grouped into 1,643 operational taxonomic units, of which about half (841) could be reliably assigned to foraminifera. The vast majority of the operational taxonomic units (nearly 90%) were either assigned to early (ancient) lineages of soft-walled, single-chambered (monothalamous) foraminifera or remained undetermined and yet possibly belong to unknown early lineages. Contrasting with the classical view of multichambered taxa dominating foraminiferal assemblages, our work reflects an unexpected diversity of monothalamous lineages that are as yet unknown using conventional micropaleontological observations. Although we can only speculate about their morphology, the immense richness of deep-sea phylotypes revealed by this study suggests that ultra-deep sequencing can improve understanding of deep-sea benthic diversity considered until now as unknowable based on a traditional taxonomic approach.

Analysis of the Elodea nuttallii transcriptome in response to mercury and cadmium pollution: development of sensitive tools for rapid ecotoxicological testing.

Toxic metals polluting aquatic ecosystems are taken up by inhabitants and accumulate in the food web, affecting species at all trophic levels. It is therefore important to have good tools to assess the level of risk represented by toxic metals in the environment. Macrophytes are potential organisms for the identification of metal-responsive biomarkers but are still underrepresented in ecotoxicology. In the present study, we used next-generation sequencing to investigate the transcriptomic response of Elodea nuttallii exposed to enhanced concentrations of Hg and Cd. We de novo assembled more than 60 000 contigs, of which we found 170 to be regulated dose-dependently by Hg and 212 by Cd. Functional analysis showed that these genes were notably related to energy and metal homeostasis. Expression analysis using nCounter of a subset of genes showed that the gene expression pattern was able to assess toxic metal exposure in complex environmental samples and was more sensitive than other end points (e.g., bioaccumulation, photosynthesis, etc.). In conclusion, we demonstrate the feasibility of using gene expression signatures for the assessment of environmental contamination, using an organism without previous genetic information. This is of interest to ecotoxicology in a wider sense given the possibility to develop specific and sensitive bioassays.

Massive production of small RNAs from a non-coding region of Cauliflower mosaic virus in plant defense and viral counter-defense.

To successfully infect plants, viruses must counteract small RNA-based host defense responses. During infection of Arabidopsis, Cauliflower mosaic pararetrovirus (CaMV) is transcribed into pregenomic 35S and subgenomic 19S RNAs. The 35S RNA is both reverse transcribed and also used as an mRNA with highly structured 600 nt leader. We found that this leader region is transcribed into long sense- and antisense-RNAs and spawns a massive quantity of 21, 22 and 24 nt viral small RNAs (vsRNAs), comparable to the entire complement of host-encoded small-interfering RNAs and microRNAs. Leader-derived vsRNAs were detected bound to the Argonaute 1 (AGO1) effector protein, unlike vsRNAs from other viral regions. Only negligible amounts of leader-derived vsRNAs were bound to AGO4. Genetic evidence showed that all four Dicer-like (DCL) proteins mediate vsRNA biogenesis, whereas the RNA polymerases Pol IV, Pol V, RDR1, RDR2 and RDR6 are not required for this process. Surprisingly, CaMV titers were not increased in dcl1/2/3/4 quadruple mutants that accumulate only residual amounts of vsRNAs. Ectopic expression of CaMV leader vsRNAs from an attenuated geminivirus led to increased accumulation of this chimeric virus. Thus, massive production of leader-derived vsRNAs does not restrict viral replication but may serve as a decoy diverting the silencing machinery from viral promoter and coding regions.

Spliced leader trapping reveals widespread alternative splicing patterns in the highly dynamic transcriptome of Trypanosoma brucei.

Trans-splicing of leader sequences onto the 5'ends of mRNAs is a widespread phenomenon in protozoa, nematodes and some chordates. Using parallel sequencing we have developed a method to simultaneously map 5'splice sites and analyze the corresponding gene expression profile, that we term spliced leader trapping (SLT). The method can be applied to any organism with a sequenced genome and trans-splicing of a conserved leader sequence. We analyzed the expression profiles and splicing patterns of bloodstream and insect forms of the parasite Trypanosoma brucei. We detected the 5' splice sites of 85% of the annotated protein-coding genes and, contrary to previous reports, found up to 40% of transcripts to be differentially expressed. Furthermore, we discovered more than 2500 alternative splicing events, many of which appear to be stage-regulated. Based on our findings we hypothesize that alternatively spliced transcripts present a new means of regulating gene expression and could potentially contribute to protein diversity in the parasite. The entire dataset can be accessed online at TriTrypDB or through: http://splicer.unibe.ch/.

Efficient mutation identification in zebrafish by microarray capturing and next generation sequencing.

Fish models like medaka, stickleback or zebrafish provide a valuable resource to study vertebrate genes. However, finding genetic variants e.g. mutations in the genome is still arduous. Here we used a combination of microarray capturing and next generation sequencing to identify the affected gene in the mozartkugelp11cv (mzlp11cv) mutant zebrafish. We discovered a 31-bp deletion in macf1 demonstrating the potential of this technique to efficiently isolate mutations in a vertebrate genome.

Sequencing of human genomes extracted from single cancer cells isolated in a valveless microfluidic device.

Sequencing the genomes of individual cells enables the direct determination of genetic heterogeneity amongst cells within a population. We have developed an injection-moulded valveless microfluidic device in which single cells from colorectal cancer derived cell lines (LS174T, LS180 and RKO) and fresh colorectal tumors have been individually trapped, their genomes extracted and prepared for sequencing using multiple displacement amplification (MDA). Ninety nine percent of the DNA sequences obtained mapped to a reference human genome, indicating that there was effectively no contamination of these samples from non-human sources. In addition, most of the reads are correctly paired, with a low percentage of singletons (0.17 ± 0.06%) and we obtain genome coverages approaching 90%. To achieve this high quality, our device design and process shows that amplification can be conducted in microliter volumes as long as the lysis is in sub-nanoliter volumes. Our data thus demonstrates that high quality whole genome sequencing of single cells can be achieved using a relatively simple, inexpensive and scalable device. Detection of genetic heterogeneity at the single cell level, as we have demonstrated for freshly obtained single cancer cells, could soon become available as a clinical tool to precisely match treatment with the properties of a patient's own tumor.

MISIS-2: A bioinformatics tool for in-depth analysis of small RNAs and representation of consensus master genome in viral quasispecies.

In most eukaryotes, small RNA (sRNA) molecules such as miRNAs, siRNAs and piRNAs regulate gene expression and repress transposons and viruses. AGO/PIWI family proteins sort functional sRNAs based on size, 5'-nucleotide and other sequence features. In plants and some animals, viral sRNAs are extremely diverse and cover the entire viral genome sequences, which allows for de novo reconstruction of a complete viral genome by deep sequencing and bioinformatics analysis of viral sRNAs. Previously, we have developed a tool MISIS to view and analyze sRNA maps of viruses and cellular genome regions which spawn multiple sRNAs. Here we describe a new release of MISIS, MISIS-2, which enables to determine and visualize a consensus sequence and count sRNAs of any chosen sizes and 5'-terminal nucleotide identities. Furthermore we demonstrate the utility of MISIS-2 for identification of single nucleotide polymorphisms (SNPs) at each position of a reference sequence and reconstruction of a consensus master genome in evolving viral quasispecies. MISIS-2 is a Java standalone program. It is freely available along with the source code at the website http://www.fasteris.com/apps.

Microbial Communities of Conducting and Respiratory Zones of Lung-Transplanted Patients.

Background: Lung transplantation (LT) is a recognized treatment for end-stage pulmonary disease. Bacteria from the recipient nasopharynx seed the new lungs leading to infections and allograft damage. Understanding the characteristics and topological variations of the microbiota may be important to apprehend the pathophysiology of allograft dysfunction. Objectives: To examine the characteristics and relationship of bacterial compositions between conducting and respiratory zones of the allograft. Methods: We performed 16S rRNA gene sequencing on bronchial aspirates (BAs) and bronchoalveolar lavages (BALs) collected in pairs in 19 patients at several time-points post-LT. Results: The respiratory zone was characterized independently of the time post-LT by a higher bacterial richness than the conducting zone (p = 0.041). The phyla Firmicutes and Proteobacteria dominated both sampling zones, with an inverse correlation between these two phyla (Spearman r = -0.830). Samples of the same pair, as well as pairs from the same individual clustered together (Pseudo-F = 3.8652, p < 0.01). Microbiota of BA and BAL were more closely related in samples from the same patient than each sample type across different patients, with variation in community structure being mainly inter-individual (p < 0.01). Both number of antibiotics administered (p < 0.01) and time interval post-LT (p < 0.01) contributed to the variation in global microbiota structure. Longitudinal analysis of BA-BAL pairs of two patients showed dynamic wave like fluctuations of the microbiota. Conclusions: Our results show that post-transplant respiratory zones harbor higher bacterial richness, but overall similar bacterial profiles as compared to conductive zones. They further support an individual microbial signature following LT.

MISIS: a bioinformatics tool to view and analyze maps of small RNAs derived from viruses and genomic loci generating multiple small RNAs.

In eukaryotes, diverse small RNA (sRNA) populations including miRNAs, siRNAs and piRNAs regulate gene expression and repress transposons, transgenes and viruses. Functional sRNAs are associated with effector proteins based on their size and nucleotide composition. The sRNA populations are currently analyzed by deep sequencing that generates millions of reads which are then mapped to a reference sequence or database. Here we developed a tool called MISIS to view and analyze sRNA maps of genomic loci and viruses which spawn multiple sRNAs. MISIS displays sRNA reads as a histogram where the x-axis indicates positions of the 5'- or 3'-terminal nucleotide of sense and antisense sRNAs, respectively, along a given reference sequence or its selected region and the y-axis the number of reads starting (for sense sRNA) or ending (for antisense sRNA) at each position. Size-classes of sRNAs can be visualized and compared separately or in combination. Thus, MISIS gives an overview of sRNA distribution along the reference sequence as well as detailed information on single sRNA species of different size-classes and abundances. MISIS reads standard BAM/SAM files outputted by mapping tools and generates table files containing counts of sRNA reads at each position of the reference sequence forward and reverse strand and for each of the chosen size-classes of sRNAs. These table files can be used by other tools such as Excel for further quantitative analysis and visualization. MISIS is a Java standalone program. It is freely available along with the source code at the following website: http://www.fasteris.com/apps.