Number Of Biological Replicates Research Articles

Epistasis, aneuploidy, and functional mutations underlie evolution of resistance to induced microtubule depolymerization.

Cells with blocked microtubule polymerization are delayed in mitosis, but eventually manage to proliferate despite substantial chromosome missegregation. While several studies have analyzed the first cell division after microtubule depolymerization, we have asked how cells cope long-term with microtubule impairment. We allowed 24 clonal populations of yeast cells with beta-tubulin mutations preventing proper microtubule polymerization, to evolve for ˜150 generations. At the end of the laboratory evolution experiment, cells had regained the ability to form microtubules and were less sensitive to microtubule-depolymerizing drugs. Whole-genome sequencing identified recurrently mutated genes, in particular for tubulins and kinesins, as well as pervasive duplication of chromosome VIII. Recreating these mutations and chromosome VIII disomy prior to evolution confirmed that they allow cells to compensate for the original mutation in beta-tubulin. Most of the identified mutations did not abolish function, but rather restored microtubule functionality. Analysis of the temporal order of resistance development in independent populations repeatedly revealed the same series of events: disomy of chromosome VIII followed by a single additional adaptive mutation in either tubulins or kinesins. Since tubulins are highly conserved among eukaryotes, our results have implications for understanding resistance to microtubule-targeting drugs widely used in cancer therapy.

Beyond fish eDNA metabarcoding: Field replicates disproportionately improve the detection of stream-associated vertebrate species

Fast, reliable, and comprehensive biodiversity monitoring data are needed for environmental decision making and management. Recent work on fish environmental DNA (eDNA) metabarcoding shows that aquatic diversity can be captured fast, reliably, and non-invasively at moderate costs. Because freshwater ecosystems act as sinks in the landscape, they also collect traces of terrestrial species via surface runoff or when specimens come into direct contact with water (e.g., for drinking purposes). Thus, fish eDNA metabarcoding data can provide information on fish but also on other, even terrestrial vertebrate species that live in riparian habitats. This data become available and may offer a much more comprehensive approach for assessing vertebrate diversity at no additional costs. Studies on how the sampling strategy affects species detection especially of stream-associated communities, however, are scarce. We therefore performed an analysis on the effects of biological replication on both fish as well as (semi-)terrestrial species detection. Along a 2-km stretch of the river Mulde (Germany), we collected 18 1L water samples and analyzed the relation of detected species richness and quantity of biological replicates taken. We detected 58 vertebrate species, of which 25 were fish and lamprey, 18 mammals, and 15 birds, which account for 50%, 24% and 7% of all native species to the German federal state of Saxony-Anhalt. However, while increasing the number of biological replicates resulted in only 25 % more detected fish and lamprey species, mammal and bird species richness increased disproportionately by 69 % and 84 %, respectively. Contrary, PCR replicates showed little stochasticity. We thus emphasize to increase the number of biological replicates when the aim is to improve general species detections. This holds especially true, when the focus is on rare aquatic taxa or on (semi-)terrestrial species, the so-called ‘bycatch’. As a clear advantage, this information can be obtained without any additional sampling or laboratory effort when the sampling strategy regarding biological replication is chosen carefully. With the consideration of frequent eDNA metabarcoding as part of national biomonitoring programs, the additional information provided by the bycatch can be used to further investigate the state of the environment and its biodiversity on a much broader scale.

High heterogeneity undermines generalization of differential expression results in RNA-Seq analysis

BackgroundRNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible.ResultsOur findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis.ConclusionsHigh heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.

An experimental design tool to optimize inference precision in data-driven mathematical models of bacterial infections in vivo.

The management of bacterial diseases calls for a detailed knowledge about the dynamic changes in host-bacteria interactions. Biological insights are gained by integrating experimental data with mechanistic mathematical models to infer experimentally unobservable quantities. This inter-disciplinary field would benefit from experiments with maximal information content yielding high-precision inference. Here, we present a computationally efficient tool for optimizing experimental design in terms of parameter inference in studies using isogenic-tagged strains. We study the effect of three experimental design factors: number of biological replicates, sampling timepoint selection and number of copies per tagged strain. We conduct a simulation study to establish the relationship between our optimality criterion and the size of parameter estimate confidence intervals, and showcase its application in a range of biological scenarios reflecting different dynamics patterns observed in experimental infections. We show that in low-variance systems with low killing and replication rates, predicting high-precision experimental designs is consistently achieved; higher replicate sizes and strategic timepoint selection yield more precise estimates. Finally, we address the question of resource allocation under constraints; given a fixed number of host animals and a constraint on total inoculum size per host, infections with fewer strains at higher copies per strain lead to higher-precision inference.

Markov chain Monte Carlo for active module identification problem

BackgroundIntegrative network methods are commonly used for interpretation of high-throughput experimental biological data: transcriptomics, proteomics, metabolomics and others. One of the common approaches is finding a connected subnetwork of a global interaction network that best encompasses significant individual changes in the data and represents a so-called active module. Usually methods implementing this approach find a single subnetwork and thus solve a hard classification problem for vertices. This subnetwork inherently contains erroneous vertices, while no instrument is provided to estimate the confidence level of any particular vertex inclusion. To address this issue, in the current study we consider the active module problem as a soft classification problem.ResultsWe propose a method to estimate probabilities of each vertex to belong to the active module based on Markov chain Monte Carlo (MCMC) subnetwork sampling. As an example of the performance of our method on real data, we run it on two gene expression datasets. For the first many-replicate expression dataset we show that the proposed approach is consistent with an existing resampling-based method. On the second dataset the jackknife resampling method is inapplicable due to the small number of biological replicates, but the MCMC method can be run and shows high classification performance.ConclusionsThe proposed method allows to estimate the probability that an individual vertex belongs to the active module as well as the false discovery rate (FDR) for a given set of vertices. Given the estimated probabilities, it becomes possible to provide a connected subgraph in a consistent manner for any given FDR level: no vertex can disappear when the FDR level is relaxed. We show, on both simulated and real datasets, that the proposed method has good computational performance and high classification accuracy.

Exploration of the sputum methylome and omics deconvolution by quadratic programming in molecular profiling of asthma and COPD: the road to sputum omics 2.0

BackgroundTo date, most studies involving high-throughput analyses of sputum in asthma and COPD have focused on identifying transcriptomic signatures of disease. No whole-genome methylation analysis of sputum cells has been performed yet. In this context, the highly variable cellular composition of sputum has potential to confound the molecular analyses.MethodsWhole-genome transcription (Agilent Human 4 × 44 k array) and methylation (Illumina 450 k BeadChip) analyses were performed on sputum samples of 9 asthmatics, 10 healthy and 10 COPD subjects. RNA integrity was checked by capillary electrophoresis and used to correct in silico for bias conferred by RNA degradation during biobank sample storage. Estimates of cell type-specific molecular profiles were derived via regression by quadratic programming based on sputum differential cell counts. All analyses were conducted using the open-source R/Bioconductor software framework.ResultsA linear regression step was found to perform well in removing RNA degradation-related bias among the main principal components of the gene expression data, increasing the number of genes detectable as differentially expressed in asthma and COPD sputa (compared to controls). We observed a strong influence of the cellular composition on the results of mixed-cell sputum analyses. Exemplarily, upregulated genes derived from mixed-cell data in asthma were dominated by genes predominantly expressed in eosinophils after deconvolution. The deconvolution, however, allowed to perform differential expression and methylation analyses on the level of individual cell types and, though we only analyzed a limited number of biological replicates, was found to provide good estimates compared to previously published data about gene expression in lung eosinophils in asthma. Analysis of the sputum methylome indicated presence of differential methylation in genomic regions of interest, e.g. mapping to a number of human leukocyte antigen (HLA) genes related to both major histocompatibility complex (MHC) class I and II molecules in asthma and COPD macrophages. Furthermore, we found the SMAD3 (SMAD family member 3) gene, among others, to lie within differentially methylated regions which has been previously reported in the context of asthma.ConclusionsIn this methodology-oriented study, we show that methylation profiling can be easily integrated into sputum analysis workflows and exhibits a strong potential to contribute to the profiling and understanding of pulmonary inflammation. Wherever RNA degradation is of concern, in silico correction can be effective in improving both sensitivity and specificity of downstream analyses. We suggest that deconvolution methods should be integrated in sputum omics analysis workflows whenever possible in order to facilitate the unbiased discovery and interpretation of molecular patterns of inflammation.

Reporting of methodologies used for clonogenic assays to determine radiosensitivity.

Radiotherapy treatment strategies should be personalized based on the radiosensitivity of individual tumors. Clonogenic assays are the gold standard method for in vitro assessment of radiosensitivity. Reproducibility is the critical factor for scientific rigor; however, this is reduced by insufficient reporting of methodologies. In reality, the reporting standards of methodologies pertaining to clonogenic assays remain unclear. To address this, we performed a literature search and qualitative analysis of the reporting of methodologies pertaining to clonogenic assays. A comprehensive literature review identified 1672 papers that report the radiosensitivity of human cancer cells based on clonogenic assays. From the identified papers, important experimental parameters (i.e. number of biological replicates, technical replicates, radiation source and dose rate) were recorded and analyzed. We found that, among the studies, (i) 30.5% did not report biological or technical replicates; (ii) 47.0% did not use biological or technical replicates; (iii) 3.8% did not report the radiation source; and (iv) 32.3% did not report the dose rate. These data suggest that reporting of methodologies pertaining to clonogenic assays in a considerable number of previously published studies is insufficient, thereby threatening reproducibility. This highlights the need to raise awareness of standardization of the methodologies used to conduct clonogenic assays.

Genetic interactions derived from high-throughput phenotyping of 6589 yeast cell cycle mutants

Over the last 30 years, computational biologists have developed increasingly realistic mathematical models of the regulatory networks controlling the division of eukaryotic cells. These models capture data resulting from two complementary experimental approaches: low-throughput experiments aimed at extensively characterizing the functions of small numbers of genes, and large-scale genetic interaction screens that provide a systems-level perspective on the cell division process. The former is insufficient to capture the interconnectivity of the genetic control network, while the latter is fraught with irreproducibility issues. Here, we describe a hybrid approach in which the 630 genetic interactions between 36 cell-cycle genes are quantitatively estimated by high-throughput phenotyping with an unprecedented number of biological replicates. Using this approach, we identify a subset of high-confidence genetic interactions, which we use to refine a previously published mathematical model of the cell cycle. We also present a quantitative dataset of the growth rate of these mutants under six different media conditions in order to inform future cell cycle models.

Discriminant analysis and feature selection in mass spectrometry imaging using constrained repeated random sampling - Cross validation (CORRS-CV)

The identification of biomarkers through Mass spectrometry imaging (MSI) is gaining popularity in the clinical field. However, considering the complexity of spectral and spatial variables faced, data mining of the hyperspectral images can be troublesome. The discovery of markers generally depends on the creation of classification models which should be validated to ensure the statistical significance of the discriminants m/z detected. Internal validation using resampling methods such as cross validation (CV) are widely used for model selection, the estimation of its generalization performance and biomarker discovery when sample sizes are limited and an independent test set is not available. Here, we introduce for first time the use of Constrained Repeated Random Subsampling CV (CORRS-CV) on multi-images for the validation of classification models on MSI. Although several aspects must be taken into account (e.g. image size, CORRS-CV∂value, the similarity across spatially close pixels, the total computation time), CORRS-CV provides more accurate estimates of the model performance than k-fold CV using of biological replicates to define the data split when the number of biological replicates is scarce and holding images back for testing is a waste of valuable information. Besides, the combined use of CORRS-CV and rank products increases the robustness of the selection of discriminant features as candidate biomarkers which is an important issue due to the increased biological, environmental and technical variabilities when analysing multiple images, especially from human tissues collected in clinical studies.

Bayesian estimation of genetic regulatory effects in high-throughput reporter assays.

High-throughput reporter assays dramatically improve our ability to assign function to noncoding genetic variants, by measuring allelic effects on gene expression in the controlled setting of a reporter gene. Unlike genetic association tests, such assays are not confounded by linkage disequilibrium when loci are independently assayed. These methods can thus improve the identification of causal disease mutations. While work continues on improving experimental aspects of these assays, less effort has gone into developing methods for assessing the statistical significance of assay results, particularly in the case of rare variants captured from patient DNA. We describe a Bayesian hierarchical model, called Bayesian Inference of Regulatory Differences, which integrates prior information and explicitly accounts for variability between experimental replicates. The model produces substantially more accurate predictions than existing methods when allele frequencies are low, which is of clear advantage in the search for disease-causing variants in DNA captured from patient cohorts. Using the model, we demonstrate a clear tradeoff between variant sequencing coverage and numbers of biological replicates, and we show that the use of additional biological replicates decreases variance in estimates of effect size, due to the properties of the Poisson-binomial distribution. We also provide a power and sample size calculator, which facilitates decision making in experimental design parameters. The software is freely available from www.geneprediction.org/bird. The experimental design web tool can be accessed at http://67.159.92.22:8080. Supplementary data are available at Bioinformatics online.

IntEREst: intron-exon retention estimator

BackgroundIn-depth study of the intron retention levels of transcripts provide insights on the mechanisms regulating pre-mRNA splicing efficiency. Additionally, detailed analysis of retained introns can link these introns to post-transcriptional regulation or identify aberrant splicing events in human diseases.ResultsWe present IntEREst, Intron–Exon Retention Estimator, an R package that supports rigorous analysis of non-annotated intron retention events (in addition to the ones annotated by RefSeq or similar databases), and support intra-sample in addition to inter-sample comparisons. It accepts binary sequence alignment/map (.bam) files as input and determines genome-wide estimates of intron retention or exon-exon junction levels. Moreover, it includes functions for comparing subsets of user-defined introns (e.g. U12-type vs U2-type) and its plotting functions allow visualization of the distribution of the retention levels of the introns. Statistical methods are adapted from the DESeq2, edgeR and DEXSeq R packages to extract the significantly more or less retained introns. Analyses can be performed either sequentially (on single core) or in parallel (on multiple cores). We used IntEREst to investigate the U12- and U2-type intron retention in human and plant RNAseq dataset with defects in the U12-dependent spliceosome due to mutations in the ZRSR2 component of this spliceosome. Additionally, we compared the retained introns discovered by IntEREst with that of other methods and studies.ConclusionIntEREst is an R package for Intron retention and exon-exon junction levels analysis of RNA-seq data. Both the human and plant analyses show that the U12-type introns are retained at higher level compared to the U2-type introns already in the control samples, but the retention is exacerbated in patient or plant samples carrying a mutated ZRSR2 gene. Intron retention events caused by ZRSR2 mutation that we discovered using IntEREst (DESeq2 based function) show considerable overlap with the retained introns discovered by other methods (e.g. IRFinder and edgeR based function of IntEREst). Our results indicate that increase in both the number of biological replicates and the depth of sequencing library promote the discovery of retained introns, but the effect of library size gradually decreases with more than 35 million reads mapped to the introns.

Identification of differentially expressed peptides in high-throughput proteomics data.

With the advent of high-throughput proteomics, the type and amount of data pose a significant challenge to statistical approaches used to validate current quantitative analysis. Whereas many studies focus on the analysis at the protein level, the analysis of peptide-level data provides insight into changes at the sub-protein level, including splice variants, isoforms and a range of post-translational modifications. Statistical evaluation of liquid chromatography-mass spectrometry/mass spectrometry peptide-based label-free differential data is most commonly performed using a t-test or analysis of variance, often after the application of data imputation to reduce the number of missing values. In high-throughput proteomics, statistical analysis methods and imputation techniques are difficult to evaluate, given the lack of gold standard data sets. Here, we use experimental and resampled data to evaluate the performance of four statistical analysis methods and the added value of imputation, for different numbers of biological replicates. We find that three or four replicates are the minimum requirement for high-throughput data analysis and confident assignment of significant changes. Data imputation does increase sensitivity in some cases, but leads to a much higher actual false discovery rate. Additionally, we find that empirical Bayes method (limma) achieves the highest sensitivity, and we thus recommend its use for performing differential expression analysis at the peptide level.

Identification of marginal causal relationships in gene networks from observational and interventional expression data.

Causal network inference is an important methodological challenge in biology as well as other areas of application. Although several causal network inference methods have been proposed in recent years, they are typically applicable for only a small number of genes, due to the large number of parameters to be estimated and the limited number of biological replicates available. In this work, we consider the specific case of transcriptomic studies made up of both observational and interventional data in which a single gene of biological interest is knocked out. We focus on a marginal causal estimation approach, based on the framework of Gaussian directed acyclic graphs, to infer causal relationships between the knocked-out gene and a large set of other genes. In a simulation study, we found that our proposed method accurately differentiates between downstream causal relationships and those that are upstream or simply associative. It also enables an estimation of the total causal effects between the gene of interest and the remaining genes. Our method performed very similarly to a classical differential analysis for experiments with a relatively large number of biological replicates, but has the advantage of providing a formal causal interpretation. Our proposed marginal causal approach is computationally efficient and may be applied to several thousands of genes simultaneously. In addition, it may help highlight subsets of genes of interest for a more thorough subsequent causal network inference. The method is implemented in an R package called MarginalCausality (available on GitHub).

The design of experiments for the transcriptome studies by high-throughput sequencing methods

The common questions in the design of the highthroughput sequencing experiments using RNA-Seq or Ribo-Seq methods are reviewed. The ENCODE guidelines (2011) as well as the recently published advances in the design of the studies of mammalian, animal and plant transcriptomes are also summarized in this review. The optimal limit of the sequencing depth does exist for the identification of almost all actively transcribed genes. This limit depends on the transcriptome size in the biological object studied. Additional sequencing does not provide any substantial additional information about the transcriptome complexity. For mammals, the optimal limit of the sequencing depth for the identification of the actively transcribed genes is equal to ~ 2 × 109 bp per biological sample. For other species, the optimal limit of the sequencing depth per biological sample is determined similarly for mammals; however, the transcriptome size and the mean RNA content in the studied object should be taken into account, in comparison to the mammalian transcriptomes. The discovery of differentially expressed genes, as well as the identification of splicing sites in the mRNA could be enhanced by increasing the number of biological samples analyzed per each experimental group. The minimal number of biological replicates per experimental group is equal to 2. However, the optimal number of biological replicates per experimental group is equal to 5–8 (similar to the experiments quantifying the expression of single genes by qRT-PCR). For the transcriptome studies, it is recommended to use the sequencing technologies that have the accuracy of sequencing ≥ 0.999 per bp. For RNASeq, it is also recommended to use the technologies that are able to produce reads equal to or larger than 75 bp, to minimize the cost of the effective identification of the sequences. The relative cost for the sequencing of the control samples could be reduced by increasing the number of experimental groups in the experiment or by combining several independent experiments with similar control groups. The present notes could be utilized during the design step in the experimental studies devoted to the research of transcriptomes.

Analysis of Differential Gene Expression Based on Bayesian Estimation of Variance

Gene expression is arguably the most important indicator of biological function. Thus identifying differentially expressed genes is one of the main aims of high throughout studies that use microarray and RNAseq platforms to study deregulated cellular pathways. There are many tools for analysing differentia gene expression from transciptomic datasets. The major challenge of this topic is to estimate gene expression variance due to the high amount of ‘background noise’ that is generated from biological equipment and the lack of biological replicates. Bayesian inference has been widely used in the bioinformatics field. In this work, we reveal that the prior knowledge employed in the Bayesian framework also helps to improve the accuracy of differential gene expression analysis when using a small number of replicates. We have developed a differential analysis tool that uses Bayesian estimation of the variance of gene expression for use with small numbers of biological replicates. Our method is more consistent when compared to the widely used cyber-t tool that successfully introduced the Bayesian framework to differential analysis. We also provide a user-friendly web based Graphic User Interface for biologists to use with microarray and RNAseq data. Bayesian inference can compensate for the instability of variance caused when using a small number of biological replicates by using pseudo replicates as prior knowledge. We also show that our new strategy to select pseudo replicates will improve the performance of the analysis. - See more at: http://www.eurekaselect.com/node/138761/article#sthash.VeK9xl5k.dpuf

Comparison of normalization and differential expression analyses using RNA-Seq data from 726 individual Drosophila melanogaster.

BackgroundA generally accepted approach to the analysis of RNA-Seq read count data does not yet exist. We sequenced the mRNA of 726 individuals from the Drosophila Genetic Reference Panel in order to quantify differences in gene expression among single flies. One of our experimental goals was to identify the optimal analysis approach for the detection of differential gene expression among the factors we varied in the experiment: genotype, environment, sex, and their interactions. Here we evaluate three different filtering strategies, eight normalization methods, and two statistical approaches using our data set. We assessed differential gene expression among factors and performed a statistical power analysis using the eight biological replicates per genotype, environment, and sex in our data set.ResultsWe found that the most critical considerations for the analysis of RNA-Seq read count data were the normalization method, underlying data distribution assumption, and numbers of biological replicates, an observation consistent with previous RNA-Seq and microarray analysis comparisons. Some common normalization methods, such as Total Count, Quantile, and RPKM normalization, did not align the data across samples. Furthermore, analyses using the Median, Quantile, and Trimmed Mean of M-values normalization methods were sensitive to the removal of low-expressed genes from the data set. Although it is robust in many types of analysis, the normal data distribution assumption produced results vastly different than the negative binomial distribution. In addition, at least three biological replicates per condition were required in order to have sufficient statistical power to detect expression differences among the three-way interaction of genotype, environment, and sex.ConclusionsThe best analysis approach to our data was to normalize the read counts using the DESeq method and apply a generalized linear model assuming a negative binomial distribution using either edgeR or DESeq software. Genes having very low read counts were removed after normalizing the data and fitting it to the negative binomial distribution. We describe the results of this evaluation and include recommended analysis strategies for RNA-Seq read count data.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2353-z) contains supplementary material, which is available to authorized users.

Quantification of co-transcriptional splicing from RNA-Seq data

During gene expression, protein-coding transcripts are shaped by multiple processing events: 5′ end capping, pre-mRNA splicing, RNA editing, and 3′ end cleavage and polyadenylation. These events are required to produce mature mRNA, which can be subsequently translated. Nearly all of these RNA processing steps occur during transcription, while the nascent RNA is still attached to the DNA template by RNA polymerase II (i.e. co-transcriptionally). Polyadenylation occurs after 3′ end cleavage or post-transcriptionally. Pre-mRNA splicing – the removal of introns and ligation of exons – can be initiated and concluded co-transcriptionally, although this is not strictly required. Recently, a number of studies using global methods have shown that the majority of splicing is co-transcriptional, yet not all published studies agree in their conclusions. Short read sequencing of RNA (RNA-Seq) is the prevailing approach to measuring splicing levels in nascent RNA, mRNA or total RNA. Here, we compare four different strategies for analyzing and quantifying co-transcriptional splicing. To do so, we reanalyze two nascent RNA-Seq datasets of the same species, but different cell type and RNA isolation procedure. Average co-transcriptional splicing values calculated on a per intron basis are similar, independent of the strategy used. We emphasize the technical requirements for identifying co-transcriptional splicing events with high confidence, e.g. how to calculate co-transcriptional splicing from nascent RNA- versus mRNA-Seq data, the number of biological replicates needed, depletion of polyA+RNA, and appropriate normalization. Finally, we present guidelines for planning a nascent RNA-Seq experiment.

A Bayesian hierarchical model to detect differentially methylated loci from single nucleotide resolution sequencing data

DNA methylation is an important epigenetic modification that has essential roles in cellular processes including gene regulation, development and disease and is widely dysregulated in most types of cancer. Recent advances in sequencing technology have enabled the measurement of DNA methylation at single nucleotide resolution through methods such as whole-genome bisulfite sequencing and reduced representation bisulfite sequencing. In DNA methylation studies, a key task is to identify differences under distinct biological contexts, for example, between tumor and normal tissue. A challenge in sequencing studies is that the number of biological replicates is often limited by the costs of sequencing. The small number of replicates leads to unstable variance estimation, which can reduce accuracy to detect differentially methylated loci (DML). Here we propose a novel statistical method to detect DML when comparing two treatment groups. The sequencing counts are described by a lognormal-beta-binomial hierarchical model, which provides a basis for information sharing across different CpG sites. A Wald test is developed for hypothesis testing at each CpG site. Simulation results show that the proposed method yields improved DML detection compared to existing methods, particularly when the number of replicates is low. The proposed method is implemented in the Bioconductor package DSS.

High proteome variation between ecotypes of Littorina saxatilis cannot be explained by tissue heterogeneity or a common-garden × ecotype effect

In small organisms where gene expression analyses are accomplished on whole specimens rather than individual tissues, the differences observed in gene expression levels between groups of samples are assumed to be caused by regulatory differences in gene expression within tissues. An alternative explanation is the lack of homogeneity distribution of the different tissues across groups of samples. In the case of the marine snail Littorina saxatilis, previous proteomic studies suggested a high differentiation in protein expression between the two (RB and SU) sympatric ecotypes existing on exposed rocky shores of the Galician coast (NW of Spain). As these ecotypes are known to differ in the proportion of muscular section (foot) contributing to the whole body, it remains to be checked which of the above explanations better explains this high proteome differentiation. Results from this new study suggest that tissue heterogeneity cannot explain the high proteomic differences observed between ecotypes of L. saxatilis. On the other hand, different estimates of proteome differentiation between ecotypes were obtained depending on the number of replicates used in each study (ranging from 7 to 30%). A reanalysis of previous published data shows a clear positive relationship between the degree of proteome differentiation observed and the number of biological replicates used in the analysis. This stresses the importance of investigating the effect of sample size on gene expression analyses. Finally, the results from the present study also discard the idea that observed proteome differences between ecotypes under similar laboratory conditions might be due to an interaction effect between common garden and the two different ecotypes. General conclusions drawn from the present study could be useful for setting up future gene expression studies in other species.

Metagenomic Predictions: From Microbiome to Complex Health and Environmental Phenotypes in Humans and Cattle

Mammals have a large cohort of endo- and ecto- symbiotic microorganisms (the microbiome) that potentially influence host phenotypes. There have been numerous exploratory studies of these symbiotic organisms in humans and other animals, often with the aim of relating the microbiome to a complex phenotype such as body mass index (BMI) or disease state. Here, we describe an efficient methodology for predicting complex traits from quantitative microbiome profiles. The method was demonstrated by predicting inflammatory bowel disease (IBD) status and BMI from human microbiome data, and enteric greenhouse gas production from dairy cattle rumen microbiome profiles. The method uses unassembled massively parallel sequencing (MPS) data to form metagenomic relationship matrices (analogous to genomic relationship matrices used in genomic predictions) to predict IBD, BMI and methane production phenotypes with useful accuracies (r = 0.423, 0.422 and 0.466 respectively). Our results show that microbiome profiles derived from MPS can be used to predict complex phenotypes of the host. Although the number of biological replicates used here limits the accuracy that can be achieved, preliminary results suggest this approach may surpass current prediction accuracies that are based on the host genome. This is especially likely for traits that are largely influenced by the gut microbiota, for example digestive tract disorders or metabolic functions such as enteric methane production in cattle.