Using reference-free compressed data structures to analyze sequencing reads from thousands of human genomes.

We are rapidly approaching the point where we have sequenced millions of human genomes. There is a pressing need for new data structures to store raw sequencing data and efficient algorithms for population scale analysis. Current reference-based data formats do not fully exploit the redundancy in population sequencing nor take advantage of shared genetic variation. In recent years, the Burrows-Wheeler transform (BWT) and FM-index have been widely employed as a full-text searchable index for read alignment and de novo assembly. We introduce the concept of a population BWT and use it to store and index the sequencing reads of 2705 samples from the 1000 Genomes Project. A key feature is that, as more genomes are added, identical read sequences are increasingly observed, and compression becomes more efficient. We assess the support in the 1000 Genomes read data for every base position of two human reference assembly versions, identifying that 3.2 Mbp with population support was lost in the transition from GRCh37 with 13.7 Mbp added to GRCh38. We show that the vast majority of variant alleles can be uniquely described by overlapping 31-mers and show how rapid and accurate SNP and indel genotyping can be carried out across the genomes in the population BWT. We use the population BWT to carry out nonreference queries to search for the presence of all known viral genomes and discover human T-lymphotropic virus 1 integrations in six samples in a recognized epidemiological distribution

The Light Responsive Transcriptome of the Zebrafish: Function and Regulation

Most organisms possess circadian clocks that are able to anticipate the day/night cycle and are reset or “entrained” by the ambient light. In the zebrafish, many organs and even cultured cell lines are directly light responsive, allowing for direct entrainment of the clock by light. Here, we have characterized light induced gene transcription in the zebrafish at several organizational levels. Larvae, heart organ cultures and cell cultures were exposed to 1- or 3-hour light pulses, and changes in gene expression were compared with controls kept in the dark. We identified 117 light regulated genes, with the majority being induced and some repressed by light. Cluster analysis groups the genes into five major classes that show regulation at all levels of organization or in different subset combinations. The regulated genes cover a variety of functions, and the analysis of gene ontology categories reveals an enrichment of genes involved in circadian rhythms, stress response and DNA repair, consistent with the exposure to visible wavelengths of light priming cells for UV-induced damage repair. Promoter analysis of the induced genes shows an enrichment of various short sequence motifs, including E- and D-box enhancers that have previously been implicated in light regulation of the zebrafish period2 gene. Heterologous reporter constructs with sequences matching these motifs reveal light regulation of D-box elements in both cells and larvae. Morpholino-mediated knock-down studies of two homologues of the D-box binding factor Tef indicate that these are differentially involved in the cell autonomous light induction in a gene-specific manner. These findings suggest that the mechanisms involved in period2 regulation might represent a more general pathway leading to light induced gene expression