Crumble: reference free lossy compression of sequence quality values

Motivation: The bulk of space taken up by NGS sequencing CRAM files consists of per-base quality values. Most of these are unnecessary for variant calling, offering an opportunity for space saving. Results: On the Syndip test set, a 17 fold reduction in the quality storage portion of a CRAM file can be achieved while maintaining variant calling accuracy. The size reduction of an entire CRAM file varied from 2.2 to 7.4 fold, depending on the non-quality content of the original file (see Supplementary Material S6 for details). Availability and implementation: Crumble is OpenSource and can be obtained from https://github.com/jkbonfield/crumble. Supplementary information: Supplementary data are available at Bioinformatics online.This work was funded by the Wellcome Trust [WT098051]

Mapping the Human Memory B Cell and Serum Neutralizing Antibody Responses to Dengue Virus Serotype 4 Infection and Vaccination

ABSTRACT The four dengue virus (DENV) serotypes are mosquito-borne flaviviruses responsible for dengue fever and dengue hemorrhagic fever. People exposed to DENV develop antibodies (Abs) that strongly neutralize the serotype responsible for infection. Historically, infection with DENV serotype 4 (DENV4) has been less common and less studied than infections with the other three serotypes. However, DENV4 has been responsible for recent large and sustained epidemics in Asia and Latin America. The neutralizing antibody responses and the epitopes targeted against DENV4 have not been characterized in human infection. In this study, we mapped and characterized epitopes on DENV4 recognized by neutralizing antibodies in people previously exposed to DENV4 infections or to a live attenuated DENV4 vaccine. To study the fine specificity of DENV4 neutralizing human antibodies, B cells from two people exposed to DENV4 were immortalized and screened to identify DENV-specific clones. Two human monoclonal antibodies (MAbs) that neutralized DENV4 were isolated, and their epitopes were finely mapped using recombinant viruses and alanine scan mutation array techniques. Both antibodies bound to quaternary structure epitopes near the hinge region between envelope protein domain I (EDI) and EDII. In parallel, to characterize the serum neutralizing antibody responses, convalescence-phase serum samples from people previously exposed to primary DENV4 natural infections or a monovalent DENV4 vaccine were analyzed. Natural infection and vaccination also induced serum-neutralizing antibodies that targeted similar epitope domains at the EDI/II hinge region. These studies defined a target of neutralizing antigenic site on DENV4 targeted by human antibodies following natural infection or vaccination. IMPORTANCE The four serotypes of dengue virus are the causative agents of dengue fever and dengue hemorrhagic fever. People exposed to primary DENV infections develop long-term neutralizing antibody responses, but these principally recognize only the infecting serotype. An effective vaccine against dengue should elicit long-lasting protective antibody responses to all four serotypes simultaneously. We and others have defined antigenic sites on the envelope (E) protein of viruses of dengue virus serotypes 1, 2, and 3 targeted by human neutralizing antibodies. The epitopes on DENV4 E protein targeted by the human neutralizing antibodies and the mechanisms of serotype 4 neutralization are poorly understood. Here, we report the properties of human antibodies that neutralize dengue virus serotype 4. People exposed to serotype 4 infections or a live attenuated serotype 4 vaccine developed neutralizing antibodies that bound to similar sites on the viral E protein. These studies have provided a foundation for developing and evaluating DENV4 vaccines

Integrative Genomics Viewer

Author Manuscript 2012 May 07.To the Editor: Rapid improvements in sequencing and array-based platforms are resulting in a flood of diverse genome-wide data, including data from exome and whole-genome sequencing, epigenetic surveys, expression profiling of coding and noncoding RNAs, single nucleotide polymorphism (SNP) and copy number profiling, and functional assays. Analysis of these large, diverse data sets holds the promise of a more comprehensive understanding of the genome and its relation to human disease. Experienced and knowledgeable human review is an essential component of this process, complementing computational approaches. This calls for efficient and intuitive visualization tools able to scale to very large data sets and to flexibly integrate multiple data types, including clinical data. However, the sheer volume and scope of data pose a significant challenge to the development of such tools.National Institute of General Medical Sciences (U.S.) (R01GM074024)National Cancer Institute (U.S.) (R21CA135827)National Human Genome Research Institute (U.S.) (U54HG003067

Non-perturbative dynamics of hot non-Abelian gauge fields: beyond leading log approximation

Many aspects of high-temperature gauge theories, such as the electroweak baryon number violation rate, color conductivity, and the hard gluon damping rate, have previously been understood only at leading logarithmic order (that is, neglecting effects suppressed only by an inverse logarithm of the gauge coupling). We discuss how to systematically go beyond leading logarithmic order in the analysis of physical quantities. Specifically, we extend to next-to-leading-log order (NLLO) the simple leading-log effective theory due to Bodeker that describes non-perturbative color physics in hot non-Abelian plasmas. A suitable scaling analysis is used to show that no new operators enter the effective theory at next-to-leading-log order. However, a NLLO calculation of the color conductivity is required, and we report the resulting value. Our NLLO result for the color conductivity can be trivially combined with previous numerical work by G. Moore to yield a NLLO result for the hot electroweak baryon number violation rate.Comment: 20 pages, 1 figur

Novel Role of Prostate Apoptosis Response-4 Tumor Suppressor in B-Cell Chronic Lymphocytic Leukemia

Prostate apoptosis response-4 (Par-4), a proapoptotic tumor suppressor protein, is downregulated in many cancers including renal cell carcinoma, glioblastoma, endometrial, and breast cancer. Par-4 induces apoptosis selectively in various types of cancer cells but not normal cells. We found that chronic lymphocytic leukemia (CLL) cells from human patients and from Eµ-Tcl1 mice constitutively express Par-4 in greater amounts than normal B-1 or B-2 cells. Interestingly, knockdown of Par-4 in human CLL-derived Mec-1 cells results in a robust increase in p21/WAF1 expression and decreased growth due to delayed G1-to-S cell-cycle transition. Lack of Par-4 also increased the expression of p21 and delayed CLL growth in Eμ-Tcl1 mice. Par-4 expression in CLL cells required constitutively active B-cell receptor (BCR) signaling, as inhibition of BCR signaling with US Food and Drug Administration (FDA)–approved drugs caused a decrease in Par-4 messenger RNA and protein, and an increase in apoptosis. In particular, activities of Lyn, a Src family kinase, spleen tyrosine kinase, and Bruton tyrosine kinase are required for Par-4 expression in CLL cells, suggesting a novel regulation of Par-4 through BCR signaling. Together, these results suggest that Par-4 may play a novel progrowth rather than proapoptotic role in CLL and could be targeted to enhance the therapeutic effects of BCR-signaling inhibitors

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software

Extent, causes, and consequences of small RNA expression variation in human adipose tissue.

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population

Bos taurus genome assembly

We present here the assembly of the bovine genome. The assembly method combines the BAC plus WGS local assembly used for the rat and sea urchin with the whole genome shotgun (WGS) only assembly used for many other animal genomes including the rhesus macaque. The assembly process consisted of multiple phases: First, BACs were assembled with BAC generated sequence, then subsequently in combination with the individual overlapping WGS reads. Different assembly parameters were tested to separately optimize the performance for each BAC assembly of the BAC and WGS reads. In parallel, a second assembly was produced using only the WGS sequences and a global whole genome assembly method. The two assemblies were combined to create a more complete genome representation that retained the high quality BAC-based local assembly information, but with gaps between BACs filled in with the WGS-only assembly. Finally, the entire assembly was placed on chromosomes using the available map information. Over 90% of the assembly is now placed on chromosomes. The estimated genome size is 2.87 Gb which represents a high degree of completeness, with 95% of the available EST sequences found in assembled contigs. The quality of the assembly was evaluated by comparison to 73 finished BACs, where the draft assembly covers between 92.5 and 100% (average 98.5%) of the finished BACs. The assembly contigs and scaffolds align linearly to the finished BACs, suggesting that misassemblies are rare. Genotyping and genetic mapping of 17,482 SNPs revealed that more than 99.2% were correctly positioned within the Btau_4.0 assembly, confirming the accuracy of the assembly. The biological analysis of this bovine genome assembly is being published, and the sequence data is available to support future bovine research

Genomewide Association Study for Determinants of HIV-1 Acquisition and Viral Set Point in HIV-1 Serodiscordant Couples with Quantified Virus Exposure

Host genetic factors may be important determinants of HIV-1 sexual acquisition. We performed a genome-wide association study (GWAS) for host genetic variants modifying HIV-1 acquisition and viral control in the context of a cohort of African HIV-1 serodiscordant heterosexual couples. To minimize misclassification of HIV-1 risk, we quantified HIV-1 exposure, using data including plasma HIV-1 concentrations, gender, and condom use.We matched couples without HIV-1 seroconversion to those with seroconversion by quantified HIV-1 exposure risk. Logistic regression of single nucleotide polymorphisms (SNPs) for 798 samples from 496 HIV-1 infected and 302 HIV-1 exposed, uninfected individuals was performed to identify factors associated with HIV-1 acquisition. In addition, a linear regression analysis was performed using SNP data from a subset (n = 403) of HIV-1 infected individuals to identify factors predicting plasma HIV-1 concentrations.After correcting for multiple comparisons, no SNPs were significantly associated with HIV-1 infection status or plasma HIV-1 concentrations.This GWAS controlling for HIV-1 exposure did not identify common host genotypes influencing HIV-1 acquisition. Alternative strategies, such as large-scale sequencing to identify low frequency variation, should be considered for identifying novel host genetic predictors of HIV-1 acquisition

The Genome Sequence of Caenorhabditis briggsae: A Platform for Comparative Genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes