Background

Following the outbreaks of 2009 pandemic H1N1 infection, rapid influenza diagnostic tests have been used to detect H1N1 infection. However, no meta-analysis has been undertaken to assess the diagnostic accuracy when this manuscript was drafted.

Methods

The literature was systematically searched to identify studies that reported the performance of rapid tests. Random effects meta-analyses were conducted to summarize the overall performance.

Results

Seventeen studies were selected with 1879 cases and 3477 non-cases. The overall sensitivity and specificity estimates of the rapid tests were 0.51 (95%CI: 0.41, 0.60) and 0.98 (95%CI: 0.94, 0.99). Studies reported heterogeneous sensitivity estimates, ranging from 0.11 to 0.88. If the prevalence was 30%, the overall positive and negative predictive values were 0.94 (95%CI: 0.85, 0.98) and 0.82 (95%CI: 0.79, 0.85). The overall specificities from different manufacturers were comparable, while there were some differences for the overall sensitivity estimates. BinaxNOW had a lower overall sensitivity of 0.39 (95%CI: 0.24, 0.57) compared to all the others (p-value < 0.001), whereas QuickVue had a higher overall sensitivity of 0.57 (95%CI: 0.50, 0.63) compared to all the others (p-value = 0.005).

Conclusions

Rapid tests have high specificity but low sensitivity and thus limited usefulness.

The RNA transcriptome varies in response to cellular differentiation as well as environmental factors, and can be characterized by the diversity and abundance of transcript isoforms. Differential transcription analysis, the detection of differences between the transcriptomes of different cells, may improve understanding of cell differentiation and development and enable the identification of biomarkers that classify disease types. The availability of high-throughput short-read RNA sequencing technologies provides in-depth sampling of the transcriptome, making it possible to accurately detect the differences between transcriptomes. In this article, we present a new method for the detection and visualization of differential transcription. Our approach does not depend on transcript or gene annotations. It also circumvents the need for full transcript inference and quantification, which is a challenging problem because of short read lengths, as well as various sampling biases. Instead, our method takes a divide-and-conquer approach to localize the difference between transcriptomes in the form of alternative splicing modules (ASMs), where transcript isoforms diverge. Our approach starts with the identification of ASMs from the splice graph, constructed directly from the exons and introns predicted from RNA-seq read alignments. The abundance of alternative splicing isoforms residing in each ASM is estimated for each sample and is compared across sample groups. A non-parametric statistical test is applied to each ASM to detect significant differential transcription with a controlled false discovery rate. The sensitivity and specificity of the method have been assessed using simulated data sets and compared with other state-of-the-art approaches. Experimental validation using qRT-PCR confirmed a selected set of genes that are differentially expressed in a lung differentiation study and a breast cancer data set, demonstrating the utility of the approach applied on experimental biological data sets. The software of DiffSplice is available at http://www.netlab.uky.edu/p/bioinfo/DiffSplice.

The response to stressful stimuli requires rapid, precise, and dynamic gene expression changes that must be coordinated across the genome. To gain insight into the temporal ordering of genome reorganization, we investigated dynamic relationships between changing nucleosome occupancy, transcription factor binding, and gene expression in Saccharomyces cerevisiae yeast responding to oxidative stress. We applied deep sequencing to nucleosomal DNA at six time points before and after hydrogen peroxide treatment and revealed many distinct dynamic patterns of nucleosome gain and loss. The timing of nucleosome repositioning was not predictive of the dynamics of downstream gene expression change but instead was linked to nucleosome position relative to transcription start sites and specific cis-regulatory elements. We measured genome-wide binding of the stress-activated transcription factor Msn2p over time and found that Msn2p binds different loci with different dynamics. Nucleosome eviction from Msn2p binding sites was common across the genome; however, we show that, contrary to expectation, nucleosome loss occurred after Msn2p binding and in fact required Msn2p. This negates the prevailing model that nucleosomes obscuring Msn2p sites regulate DNA access and must be lost before Msn2p can bind DNA. Together, these results highlight the complexities of stress-dependent chromatin changes and their effects on gene expression.

Summary

DNA methylation has emerged as an important hallmark of epigenetics. Numerous platforms including tiling arrays and next generation sequencing, and experimental protocols are available for profiling DNA methylation. Similar to other tiling array data, DNA methylation data shares the characteristics of inherent correlation structure among nearby probes. However, unlike gene expression or protein DNA binding data, the varying CpG density which gives rise to CpG island, shore and shelf definition provides exogenous information in detecting differential methylation. This paper aims to introduce a robust testing and probe ranking procedure based on a non-homogeneous hidden Markov model that incorporates the above-mentioned features for detecting differential methylation. We revisit the seminal work of Sun and Cai (2009, J. R. Stat. Soc. B. 71, 393-424) and propose modeling the non-null using a non-parametric symmetric distribution in two-sided hypothesis testing. We show that this model improves probe ranking and is robust to model misspecification based on extensive simulation studies. We further illustrate that our proposed framework achieves good operating characteristics as compared to commonly used methods in real DNA methylation data that aims to detect differential methylation sites.

Obesity has reached epidemic proportions worldwide. Several animal models of obesity exist, but studies are lacking that compare traditional lard-based high fat diets (HFD) to “Cafeteria diets" (CAF) consisting of nutrient poor human junk food. Our previous work demonstrated the rapid and severe obesogenic and inflammatory consequences of CAF compared to HFD including rapid weight gain, markers of Metabolic Syndrome, multi-tissue lipid accumulation, and dramatic inflammation. To identify potential mediators of CAF-induced obesity and Metabolic Syndrome, we used metabolomic analysis to profile serum, muscle, and white adipose from rats fed CAF, HFD, or standard control diets. Principle component analysis identified elevations in clusters of fatty acids and acylcarnitines. These increases in metabolites were associated with systemic mitochondrial dysfunction that paralleled weight gain, physiologic measures of Metabolic Syndrome, and tissue inflammation in CAF-fed rats. Spearman pairwise correlations between metabolites, physiologic, and histologic findings revealed strong correlations between elevated markers of inflammation in CAF-fed animals, measured as crown like structures in adipose, and specifically the pro-inflammatory saturated fatty acids and oxidation intermediates laurate and lauroyl carnitine. Treatment of bone marrow-derived macrophages with lauroyl carnitine polarized macrophages towards the M1 pro-inflammatory phenotype through downregulation of AMPK and secretion of pro-inflammatory cytokines. Results presented herein demonstrate that compared to a traditional HFD model, the CAF diet provides a robust model for diet-induced human obesity, which models Metabolic Syndrome-related mitochondrial dysfunction in serum, muscle, and adipose, along with pro-inflammatory metabolite alterations. These data also suggest that modifying the availability or metabolism of saturated fatty acids may limit the inflammation associated with obesity leading to Metabolic Syndrome.

Summary

DNA methylation, an epigenetic alteration typically occurring early in cancer development, could aid in the molecular diagnosis of melanoma. We determined technical feasibility for high-throughput DNA-methylation array-based profiling using formalin-fixed paraffin-embedded tissues for selection of candidate DNA-methylation differences between melanomas and nevi. Promoter methylation was evaluated in 27 common benign nevi and 22 primary invasive melanomas using a 1505 CpG-site microarray. Unsupervised hierarchical clustering distinguished melanomas from nevi; and 26 CpG sites in 22 genes were identified with significantly different methylation levels between melanomas and nevi after adjustment for age, sex, and multiple comparisons and with β-value differences of ≥ 0.2. Prediction Analysis for Microarrays identified 12 CpG loci that were highly predictive of melanoma, with area under the receiver operating characteristic curves of greater than 0.95. Of our panel of 22 genes, 14 were statistically significant in an independent sample set of 29 nevi (including dysplastic nevi) and 25 primary invasive melanomas after adjustment for age, sex, and multiple comparisons. This first report of a DNA-methylation signature discriminating melanomas from nevi indicates that DNA methylation appears promising as an additional tool for enhancing melanoma diagnosis.

Motivation: The Illumina BeadArray is a popular platform for profiling DNA methylation, an important epigenetic event associated with gene silencing and chromosomal instability. However, current approaches rely on an arbitrary detection P-value cutoff for excluding probes and samples from subsequent analysis as a quality control step, which results in missing observations and information loss. It is desirable to have an approach that incorporates the whole data, but accounts for the different quality of individual observations.

Results: We first investigate and propose a statistical framework for removing the source of biases in Illumina Methylation BeadArray based on several positive control samples. We then introduce a weighted model-based clustering called LumiWCluster for Illumina BeadArray that weights each observation according to the detection P-values systematically and avoids discarding subsets of the data. LumiWCluster allows for discovery of distinct methylation patterns and automatic selection of informative CpG loci. We demonstrate the advantages of LumiWCluster on two publicly available Illumina GoldenGate Methylation datasets (ovarian cancer and hepatocellular carcinoma).

Availability: R package LumiWCluster can be downloaded from http://www.unc.edu/~pfkuan/LumiWCluster

Contact: pfkuan@bios.unc.edu

Supplementary information: Supplementary data are available at Bioinformatics online.

Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is rapidly replacing chromatin immunoprecipitation combined with genome-wide tiling array analysis (ChIP-chip) as the preferred approach for mapping transcription-factor binding sites and chromatin modifications. The state of the art for analyzing ChIP-seq data relies on using only reads that map uniquely to a relevant reference genome (uni-reads). This can lead to the omission of up to 30% of alignable reads. We describe a general approach for utilizing reads that map to multiple locations on the reference genome (multi-reads). Our approach is based on allocating multi-reads as fractional counts using a weighted alignment scheme. Using human STAT1 and mouse GATA1 ChIP-seq datasets, we illustrate that incorporation of multi-reads significantly increases sequencing depths, leads to detection of novel peaks that are not otherwise identifiable with uni-reads, and improves detection of peaks in mappable regions. We investigate various genome-wide characteristics of peaks detected only by utilization of multi-reads via computational experiments. Overall, peaks from multi-read analysis have similar characteristics to peaks that are identified by uni-reads except that the majority of them reside in segmental duplications. We further validate a number of GATA1 multi-read only peaks by independent quantitative real-time ChIP analysis and identify novel target genes of GATA1. These computational and experimental results establish that multi-reads can be of critical importance for studying transcription factor binding in highly repetitive regions of genomes with ChIP-seq experiments.

Author Summary

Annotating repetitive regions of genomes experimentally is a challenging task. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) provides valuable data for characterizing repetitive regions of genomes in terms of transcription factor binding. Although ChIP-seq technology has been maturing, available ChIP-seq analysis methods and software rely on discarding sequence reads that map to multiple locations on the reference genome (multi-reads), thereby generating a missed opportunity for assessing transcription factor binding to highly repetitive regions of genomes. We develop a computational algorithm that takes multi-reads into account in ChIP-seq analysis. We show with computational experiments that multi-reads lead to significant increase in sequencing depths and identification of binding regions that are otherwise not identifiable when only reads that uniquely map to the reference genome (uni-reads) are used. In particular, we show that the number of binding regions identified can increase up to 36%. We support our computational predictions with independent quantitative real-time ChIP validation of binding regions identified only when multi-reads are incorporated in the analysis of a mouse GATA1 ChIP-seq experiment.

Previous studies described functional roles for Rho GDP dissociation inhibitor 2 (RhoGDI2) in bladder, gastric and breast cancers. However, only limited expression and no functional analyses have been done for RhoGDI2 in ovarian cancer. We determined RhoGDI2 protein expression and function in ovarian cancer. First, protein gel blot analysis was performed to determine the expression levels of RhoGDI2 in ovarian cells lines. RhoGDI2 but not RhoGDI1 protein expression levels varied widely in ovarian carcinoma cell lines, with elevated levels seen in Ras-transformed ovarian epithelial cells. Next, immunohistochemistry was performed to detect RhoGDI2 expression in patient samples of ovarian cysts and ovarian cancer with known histological subtype, stage, grade and outcome. RhoGDI2 protein was significantly overexpressed in high-grade compared with low-grade ovarian cancers, correlated with histological subtype, and did not correlate with stage of ovarian cancer nor between carcinomas and benign cysts. Unexpectedly, stable suppression of RhoGDI2 protein expression in HeyA8 ovarian cancer cells increased anchorage-independent growth and Matrigel invasion in vitro and in tail-vein lung colony metastatic growth in vivo. Finally, we found that RhoGDI2 stably-associated preferentially with Rac1 and suppression of RhoGDI2 expression resulted in decreased Rac1 activity and Rac-associated JNK and p38 mitogenactivated protein kinase signaling. RhoGDI2 antagonizes the invasive and metastatic phenotype of HeyA8 ovarian cancer cells. In summary, our results suggest significant cell context differences in RhoGDI2 function in cancer cell growth.

The ability to map individual nucleosomes accurately across genomes enables the study of relationships between dynamic changes in nucleosome positioning/occupancy and gene regulation. However, the highly heterogeneous nature of nucleosome densities across genomes and short linker regions pose challenges in mapping nucleosome positions based on high-throughput microarray data of micrococcal nuclease (MNase) digested DNA. Previous works rely on additional detrending and careful visual examination to detect low-signal nucleosomes, which may exist in a subpopulation of cells. We propose a non-homogeneous hidden-state model based on first order differences of experimental data along genomic coordinates that bypasses the need for local detrending and can automatically detect nucleosome positions of various occupancy levels. Our proposed approach is applicable to both low and high resolution MNase-Chip and MNase-Seq (high throughput sequencing) data, and is able to map nucleosome-linker boundaries accurately. This automated algorithm is also computationally efficient and only requires a simple preprocessing step. We provide several examples illustrating the pitfalls of existing methods, the difficulties of detrending the observed hybridization signals and demonstrate the advantages of utilizing first order differences in detecting nucleosome occupancies via simulations and case studies involving MNase-Chip and MNase-Seq data of nucleosome occupancy in yeast S. cerevisiae.

Whole genome tiling arrays at a user specified resolution are becoming a versatile tool in genomics. Chromatin immunoprecipitation on microarrays (ChIP-chip) is a powerful application of these arrays. Although there is an increasing number of methods for analyzing ChIP-chip data, perhaps the most simple and commonly used one, due to its computational efficiency, is testing with a moving average statistic. Current moving average methods assume exchangeability of the measurements within an array. They are not tailored to deal with the issues due to array designs such as overlapping probes that result in correlated measurements. We investigate the correlation structure of data from such arrays and propose an extension of the moving average testing via a robust and rapid method called CMARRT. We illustrate the pitfalls of ignoring the correlation structure in simulations and a case study. Our approach is implemented as an R package called CMARRT and can be used with any tiling array platform.

The ability to map individual nucleosomes accurately across genomes enables the study of relationships between dynamic changes in nucleosome positioning/occupancy and gene regulation. However, the highly heterogeneous nature of nucleosome densities across genomes and short linker regions pose challenges in mapping nucleosome positions based on high-throughput microarray data of micrococcal nuclease (MNase) digested DNA. Previous works rely on additional detrending and careful visual examination to detect low-signal nucleosomes, which may exist in a subpopulation of cells. We propose a non-homogeneous hidden-state model based on first order differences of experimental data along genomic coordinates that bypasses the need for local detrending and can automatically detect nucleosome positions of various occupancy levels. Our proposed approach is applicable to both low and high resolution MNase-Chip and MNase-Seq (high throughput sequencing) data, and is able to map nucleosome-linker boundaries accurately. This automated algorithm is also computationally efficient and only requires a simple preprocessing step. We provide several examples illustrating the pitfalls of existing methods, the difficulties of detrending the observed hybridization signals and demonstrate the advantages of utilizing first order differences in detecting nucleosome occupancies via simulations and case studies involving MNase-Chip and MNase-Seq data of nucleosome occupancy in yeast S. cerevisiae.

Background

Chromatin immunoprecipitation combined with DNA microarrays (ChIP-chip) is an assay used for investigating DNA-protein-binding or post-translational chromatin/histone modifications. As with all high-throughput technologies, it requires thorough bioinformatic processing of the data for which there is no standard yet. The primary goal is to reliably identify and localize genomic regions that bind a specific protein. Further investigation compares binding profiles of functionally related proteins, or binding profiles of the same proteins in different genetic backgrounds or experimental conditions. Ultimately, the goal is to gain a mechanistic understanding of the effects of DNA binding events on gene expression.

Results

We present a free, open-source R/Bioconductor package Starr that facilitates comparative analysis of ChIP-chip data across experiments and across different microarray platforms. The package provides functions for data import, quality assessment, data visualization and exploration. Starr includes high-level analysis tools such as the alignment of ChIP signals along annotated features, correlation analysis of ChIP signals with complementary genomic data, peak-finding and comparative display of multiple clusters of binding profiles. It uses standard Bioconductor classes for maximum compatibility with other software. Moreover, Starr automatically updates microarray probe annotation files by a highly efficient remapping of microarray probe sequences to an arbitrary genome.

Conclusion

Starr is an R package that covers the complete ChIP-chip workflow from data processing to binding pattern detection. It focuses on the high-level data analysis, e.g., it provides methods for the integration and combined statistical analysis of binding profiles and complementary functional genomics data. Starr enables systematic assessment of binding behaviour for groups of genes that are alingned along arbitrary genomic features.

RNA interference (RNAi) is a modality in which small double-stranded RNA molecules (siRNAs) designed to lead to the degradation of specific mRNAs are introduced into cells or organisms. siRNA libraries have been developed in which siRNAs targeting virtually every gene in the human genome are designed, synthesized and are presented for introduction into cells by transfection in a microtiter plate array. These siRNAs can then be transfected into cells using high-throughput screening (HTS) methodologies. The goal of RNAi HTS is to identify a set of siRNAs that inhibit or activate defined cellular phenotypes. The commonly used analysis methods including median ± kMAD have issues about error rates in multiple hypothesis testing and plate-wise versus experiment-wise analysis. We propose a methodology based on a Bayesian framework to address these issues. Our approach allows for sharing of information across plates in a plate-wise analysis, which obviates the need for choosing either a plate-wise or experimental-wise analysis. The proposed approach incorporates information from reliable controls to achieve a higher power and a balance between the contribution from the samples and control wells. Our approach provides false discovery rate (FDR) control to address multiple testing issues and it is robust to outliers.

Well-defined relationships between oligonucleotide properties and hybridization signal intensities (HSI) can aid chip design, data normalization and true biological knowledge discovery. We clarify these relationships using the data from two microarray experiments containing over three million probes from 48 high-density chips. We find that melting temperature (Tm) has the most significant effect on HSI while length for the long oligonucleotides studied has very little effect. Analysis of positional effect using a linear model provides evidence that the protruding ends of probes contribute more than tethered ends to HSI, which is further validated by specifically designed match fragment sliding and extension experiments. The impact of sequence similarity (SeqS) on HSI is not significant in comparison with other oligonucleotide properties. Using regression and regression tree analysis, we prioritize these oligonucleotide properties based on their effects on HSI. The implications of our discoveries for the design of unbiased oligonucleotides are discussed. We propose that isothermal probes designed by varying the length is a viable strategy to reduce sequence bias, though imposing selection constraints on other oligonucleotide properties is also essential.

DNA methylation, an epigenetic alteration typically occurring early in cancer development, could aid in the molecular diagnosis of melanoma. We determined technical feasibility for high-throughput DNA-methylation array-based profiling using formalin-fixed paraffin-embedded tissues for selection of candidate DNA-methylation differences between melanomas and nevi. Promoter methylation was evaluated in 27 common benign nevi and 22 primary invasive melanomas using a 1505 CpG site microarray. Unsupervised hierarchical clustering distinguished melanomas from nevi; 26 CpG sites in 22 genes were identified with significantly different methylation levels between melanomas and nevi after adjustment for age, sex, and multiple comparisons and with β-value differences of ≥0.2. Prediction analysis for microarrays identified 12 CpG loci that were highly predictive of melanoma, with area under the receiver operating characteristic curves of >0.95. Of our panel of 22 genes, 14 were statistically significant in an independent sample set of 29 nevi (including dysplastic nevi) and 25 primary invasive melanomas after adjustment for age, sex, and multiple comparisons. This first report of a DNA-methylation signature discriminating melanomas from nevi indicates that DNA methylation appears promising as an additional tool for enhancing melanoma diagnosis.

Background

EBV DNA is found within the malignant cells of 10% of gastric cancers. Modern molecular technology facilitates identification of virus-related biochemical effects that could assist in early diagnosis and disease management.

Methods

In this study, RNA expression profiling was performed on 326 macrodissected paraffin-embedded tissues including 204 cancers and, when available, adjacent non-malignant mucosa. Nanostring nCounter probes targeted 96 RNAs (20 viral, 73 human, and 3 spiked RNAs).

Results

In 182 tissues with adequate housekeeper RNAs, distinct profiles were found in infected versus uninfected cancers, and in malignant versus adjacent benign mucosa. EBV-infected gastric cancers expressed nearly all of the 18 latent and lytic EBV RNAs in the test panel. Levels of EBER1 and EBER2 RNA were highest and were proportional to the quantity of EBV genomes as measured by Q-PCR. Among protein coding EBV RNAs, EBNA1 from the Q promoter and BRLF1 were highly expressed while EBNA2 levels were low positive in only 6/14 infected cancers. Concomitant upregulation of cellular factors implies that virus is not an innocent bystander but rather is linked to NFKB signaling (FCER2, TRAF1) and immune response (TNFSF9, CXCL11, IFITM1, FCRL3, MS4A1 and PLUNC), with PPARG expression implicating altered cellular metabolism. Compared to adjacent non-malignant mucosa, gastric cancers consistently expressed INHBA, SPP1, THY1, SERPINH1, CXCL1, FSCN1, PTGS2 (COX2), BBC3, ICAM1, TNFSF9, SULF1, SLC2A1, TYMS, three collagens, the cell proliferation markers MYC and PCNA, and EBV BLLF1 while they lacked CDH1 (E-cadherin), CLDN18, PTEN, SDC1 (CD138), GAST (gastrin) and its downstream effector CHGA (chromogranin). Compared to lymphoepithelioma-like carcinoma of the uterine cervix, gastric cancers expressed CLDN18, EPCAM, REG4, BBC3, OLFM4, PPARG, and CDH17 while they had diminished levels of IFITM1 and HIF1A. The druggable targets ERBB2 (Her2), MET, and the HIF pathway, as well as several other potential pharmacogenetic indicators (including EBV infection itself, as well as SPARC, TYMS, FCGR2B and REG4) were identified in some tumor specimens.

Conclusion

This study shows how modern molecular technology applied to archival fixed tissues yields novel insights into viral oncogenesis that could be useful in managing affected patients.