RNA-Seq Reveals an Integrated Immune Response in Nucleated Erythrocytes

Background: Throughout the primary literature and within textbooks, the erythrocyte has been tacitly accepted to have maintained a unique physiological role; namely gas transport and exchange. In non-mammalian vertebrates, nucleated erythrocytes are present in circulation throughout the life cycle and a fragmented series of observations in mammals support a potential role in non-respiratory biological processes. We hypothesised that nucleated erythrocytes could actively participate via ligand-induced transcriptional re-programming in the immune response. Methodology/Principal Findings: Nucleated erythrocytes from both fish and birds express and regulate specific pattern recognition receptor (PRR) mRNAs and, thus, are capable of specific pathogen associated molecular pattern (PAMP) detection that is central to the innate immune response. In vitro challenge with diverse PAMPs led to de novo specific mRNA synthesis of both receptors and response factors including interferon-alpha (IFNα) that exhibit a stimulus-specific polysomal shift supporting active translation. RNA-Seq analysis of the PAMP (Poly (I:C), polyinosinic:polycytidylic acid)-erythrocyte response uncovered diverse cohorts of differentially expressed mRNA transcripts related to multiple physiological systems including the endocrine, reproductive and immune. Moreover, erythrocyte-derived conditioned mediums induced a type-1 interferon response in macrophages thus supporting an integrative role for the erythrocytes in the immune response. Conclusions/Significance: We demonstrate that nucleated erythrocytes in non-mammalian vertebrates spanning significant phylogenetic distance participate in the immune response. RNA-Seq studies highlight a mRNA repertoire that suggests a previously unrecognized integrative role for the erythrocytes in other physiological systems

Comprehensive Analysis of MGMT Promoter Methylation: Correlation with MGMT Expression and Clinical Response in GBM

O6-methylguanine DNA-methyltransferase (MGMT) promoter methylation has been identified as a potential prognostic marker for glioblastoma patients. The relationship between the exact site of promoter methylation and its effect on gene silencing, and the patient's subsequent response to therapy, is still being defined. The aim of this study was to comprehensively characterize cytosine-guanine (CpG) dinucleotide methylation across the entire MGMT promoter and to correlate individual CpG site methylation patterns to mRNA expression, protein expression, and progression-free survival. To best identify the specific MGMT promoter region most predictive of gene silencing and response to therapy, we determined the methylation status of all 97 CpG sites in the MGMT promoter in tumor samples from 70 GBM patients using quantitative bisulfite sequencing. We next identified the CpG site specific and regional methylation patterns most predictive of gene silencing and improved progression-free survival. Using this data, we propose a new classification scheme utilizing methylation data from across the entire promoter and show that an analysis based on this approach, which we call 3R classification, is predictive of progression-free survival (HR  = 5.23, 95% CI [2.089–13.097], p<0.0001). To adapt this approach to the clinical setting, we used a methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) test based on the 3R classification and show that this test is both feasible in the clinical setting and predictive of progression free survival (HR  = 3.076, 95% CI [1.301–7.27], p = 0.007). We discuss the potential advantages of a test based on this promoter-wide analysis and compare it to the commonly used methylation-specific PCR test. Further prospective validation of these two methods in a large independent patient cohort will be needed to confirm the added value of promoter wide analysis of MGMT methylation in the clinical setting

Conducting robust ecological analyses with climate data

Although the number of studies discerning the impact of climate change on ecological systems continues to increase, there has been relatively little sharing of the lessons learnt when accumulating this evidence. At a recent workshop entitled ‘Using climate data in ecological research’ held at the UK Met Office, ecologists and climate scientists came together to discuss the robust analysis of climate data in ecology. The discussions identified three common pitfalls encountered by ecologists: 1) selection of inappropriate spatial resolutions for analysis; 2) improper use of publically available data or code; and 3) insufficient representation of the uncertainties behind the adopted approach. Here, we discuss how these pitfalls can be avoided, before suggesting ways that both ecology and climate science can move forward. Our main recommendation is that ecologists and climate scientists collaborate more closely, on grant proposals and scientific publications, and informally through online media and workshops. More sharing of data and code (e.g. via online repositories), lessons and guidance would help to reconcile differing approaches to the robust handling of data. We call on ecologists to think critically about which aspects of the climate are relevant to their study system, and to acknowledge and actively explore uncertainty in all types of climate data. And we call on climate scientists to make simple estimates of uncertainty available to the wider research community. Through steps such as these, we will improve our ability to robustly attribute observed ecological changes to climate or other factors, while providing the sort of influential, comprehensive analyses that efforts to mitigate and adapt to climate change so urgently require

Novel Avian Influenza H7N3 Strain Outbreak, British Columbia

Genome sequences of chicken (low pathogenic avian influenza [LPAI] and highly pathogenic avian influenza [HPAI]) and human isolates from a 2004 outbreak of H7N3 avian influenza in Canada showed a novel insertion in the HA0 cleavage site of the human and HPAI isolate. This insertion likely occurred by recombination between the hemagglutination and matrix genes in the LPAI virus

Meraculous: De Novo Genome Assembly with Short Paired-End Reads

We describe a new algorithm, meraculous, for whole genome assembly of deep paired-end short reads, and apply it to the assembly of a dataset of paired 75-bp Illumina reads derived from the 15.4 megabase genome of the haploid yeast Pichia stipitis. More than 95% of the genome is recovered, with no errors; half the assembled sequence is in contigs longer than 101 kilobases and in scaffolds longer than 269 kilobases. Incorporating fosmid ends recovers entire chromosomes. Meraculous relies on an efficient and conservative traversal of the subgraph of the k-mer (deBruijn) graph of oligonucleotides with unique high quality extensions in the dataset, avoiding an explicit error correction step as used in other short-read assemblers. A novel memory-efficient hashing scheme is introduced. The resulting contigs are ordered and oriented using paired reads separated by ∼280 bp or ∼3.2 kbp, and many gaps between contigs can be closed using paired-end placements. Practical issues with the dataset are described, and prospects for assembling larger genomes are discussed

Big Genomes Facilitate the Comparative Identification of Regulatory Elements

The identification of regulatory sequences in animal genomes remains a significant challenge. Comparative genomic methods that use patterns of evolutionary conservation to identify non-coding sequences with regulatory function have yielded many new vertebrate enhancers. However, these methods have not contributed significantly to the identification of regulatory sequences in sequenced invertebrate taxa. We demonstrate here that this differential success, which is often attributed to fundamental differences in the nature of vertebrate and invertebrate regulatory sequences, is instead primarily a product of the relatively small size of sequenced invertebrate genomes. We sequenced and compared loci involved in early embryonic patterning from four species of true fruit flies (family Tephritidae) that have genomes four to six times larger than those of Drosophila melanogaster. Unlike in Drosophila, where virtually all non-coding DNA is highly conserved, blocks of conserved non-coding sequence in tephritids are flanked by large stretches of poorly conserved sequence, similar to what is observed in vertebrate genomes. We tested the activities of nine conserved non-coding sequences flanking the even-skipped gene of the teprhitid Ceratis capitata in transgenic D. melanogaster embryos, six of which drove patterns that recapitulate those of known D. melanogaster enhancers. In contrast, none of the three non-conserved tephritid non-coding sequences that we tested drove expression in D. melanogaster embryos. Based on the landscape of non-coding conservation in tephritids, and our initial success in using conservation in tephritids to identify D. melanogaster regulatory sequences, we suggest that comparison of tephritid genomes may provide a systematic means to annotate the non-coding portion of the D. melanogaster genome. We also propose that large genomes be given more consideration in the selection of species for comparative genomics projects, to provide increased power to detect functional non-coding DNAs and to provide a less biased view of the evolution and function of animal genomes